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Chapter  I 

Introduction 


Purpose  and  Objective  of  This  Study 

Silane,  SiH4,  is  a  gas  which  has  important  application  in  the 
manufacture  of  thin  films  of  hydrogenated  amorphous  silicon  (a-Si:H)  for 
solar  cells,  xerographic  photoreceptors,  and  other  imaging  detectors.  It  is 
also  used  extensively  for  the  production  of  thin  films  of  silicon  nitride 
and  silicon  oxide  employed  in  the  fabrication  of  microelectronic  devices. 
One  of  the  most  widely  used  methods  for  the  production  of  these  films  is 
that  of  Plasma  Enhanced  Chemical  Vapor  Deposition  (PECVD).  In  this 
production  technique,  a  low  pressure  electric  discharge  is  applied  to  a  gas 
mixture  containing  SiH4  and  a  thin  film  of  the  desired  material  (depending 
on  the  other  gases  employed)  is  grown  on  a  substrate.  For  example,  in  the 
case  of  amorphous  silicon  film  production,  the  SiH4  is  diluted  in  a  rare  gas 
(usually  argon)1  and  for  the  production  of  silicon  nitride  films,  the  SiH4  is 
mixed  with  either  nitrogen  or  ammonia2.  In  all  of  these  PECVD  systems, 
one  of  the  key  steps  in  the  production  process  must  be  the  dissociation  of 
SiH4  since  in  all  cases  the  films  contain  proportionally  less  hydrogen  than 
is  present  in  SiH4.  The  rate  of  SiH4  dissociation  under  PECVD  conditions  is 
thus  an  important  parameter  governing  film  growth.  The  dissociation  rate 
under  a  given  set  of  conditions  is  strongly  influenced  by  the  mechanism  by 
which  SiH4  is  dissociated.  These  include  electron  impact,  ion-molecule 
reactions,  and  reactions  involving  neutral  heavy  particles.  Interest  in  the 
work  presented  here  began  in  this  laboratory  in  1983  while  studying 
electron  impact  dissociation  of  SiH4  in  a  direct  current  discharge3.  While 
most  PECVD  reactors  involve  radio  frequency  discharges,  a  dc  discharge 
was  chosen  to  begin  the  investigation  because  of  its  more  well  understood 
characteristics.  The  details  of  this  work  will  be  discussed  later  in  this 
chapter,  but  to  summarize,  it  was  found  that  the  rate  of  SiH4  dissociation 
when  diluted  in  nitrogen  was  far  greater  than  expected  from  electron 
impact.  Some  other  mechanism  had  to  be  responsible  for  the  observed 
dissociation.  In  an  effort  to  determine  what  the  driver  for  the  dissocia- 
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work  on  SiFU  reactions  in  active  nitrogen  is  given  at  the  end  of  the 
chapter. 

In  chapter  II,  the  experimental  apparatus  is  described  along  with 
calibration  procedures  where  appropriate.  The  chapter  begins  with  a 
description  of  the  flowing  afterglow  apparatus  along  with  its  associated 
flow  and  pressure  controls.  This  is  followed  by  a  description  of  the 
Tunable  Diode  Laser  (TDL)  system  used  to  measure  SiFU  density  and  the 
calibration  procedure  used  in  the  measurement.  To  measure  S1H4  dissocia¬ 
tion  in  the  afterglow  region,  the  density  was  measured  from  the  in¬ 
tegrated  absorbance  of  an  isolated  SiH4  absorption  line.  The  integrated 
absorbance  was  determined  by  least-squares  fitting  a  Voigt  profile  to  the 
measured  absorption  data  and  then  comparing  these  n.easurements  with 
data  from  a  known  gas  mixture.  For  calibration,  the  line  intensity  of  the 
Q(3),A2,1  line  of  the  v4  band  of  SiH4  was  measured  absolutely  using  a 
known  gas  mixture  of  SiH4  in  argon.  No  detailed  line  intensity  measure¬ 
ments  are  available  in  the  literature  for  SiFU.  The  ultraviolet/visible 
monochromator  system  with  associated  photomultiplier  and  signal 
processing  system  are  described  along  with  the  relative  intensity  calibra¬ 
tion  procedure  used.  Finally,  a  description  is  given  of  the  Fourier  Trans¬ 
form  Infrared  (FTIR)  spectrometer  system  with  long  path  absorption  cell 
used  for  making  infrared  absorption  measurements  of  the  gases  leaving 
the  reaction  zone. 

In  chapter  III  the  results  of  the  measurements  and  the  analysis  of 
the  data  are  presented.  A  general  description  is  given  of  the  behavior  of 
the  observed  chemiluminescence  as  it  was  affected  bv  pressure,  SiFU 
flow,  and  discharge  current.  These  qualitative  observations  were  very 
important  in  the  overall  formulation  of  the  reaction  model.  Specifically, 
the  movement  and  shape  of  the  reaction  flame  and  the  observations  of 
pulsing  of  the  reaction  flame  were  key  indicators  of  the  presence  of  a 
chain  reaction  in  the  overall  scheme.  Following  the  general  observations 
emission  spectra  of  the  chemiluminescence  are  given  and  key  emitting 
species  are  identified.  The  primary  spectral  features  in  the  reaction  flame 
at  pressures  above  2  torr  are  the  B-X  band  system  of  SiN  and  a  number  of 
lines  due  to  atomic  silicon.  In  addition,  there  is  an  unidentified  broadband 
continuum  emitter  in  the  range  of  450-870nm.  No  NHX,  Fl2,  FI,  or  SiFI 
emission  could  be  identified.  Spectra  taken  over  the  pressure  range  0.5-8 
torr  are  presented.  At  pressures  below  2  torr  the  dominant  spectral 
feature  is  the  aforementioned  broadband  continuum  while  at  pressures 
above  2  torr,  the  domain  features  are  the  SiN  and  Si  emission.  The 
dominance  of  both  the  Si  and  SiN  emissions  at  higher  pressures  indicates 
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that  these  emissions  increase  faster  than  linear  with  pressure.  Difficul¬ 
ties  in  defining  the  emitting  volume  prevented  a  precise  determination  of 
the  pressure  dependence.  Measurements  of  the  SiN  B-state  vibrational 
populations  observed  in  the  higher  pressure  chemiluminescence  are  also 
presented.  At  these  higher  pressures,  the  measured  B-state  vibrational 
populations  show  a  peak  in  the  populations  at  vibrational  level  v=5.  This 
type  of  distribution  has  been  seen  previously  in  experiments  where  SiCI4 
and  SiBr4  were  added  to  active  nitrogen.  Following  this,  dissociation 
measurements  of  SiH4  are  presented  and  compared  with  the  first  positive 
emission  from  the  afterglow.  The  measurements  show  that  the  amount  of 
SiH4  dissociated  at  the  "zero  remaining  limit"  (when  the  initial  SiH4 
density  is  such  that  exactly  this  amount  can  be  dissociated)  is  shown  to 
be  proportional  to  the  square  root  of  the  intensity  of  the  N2  first  positive 
emission  in  the  afterglow.  This  indicates  that  the  dissociation  is  propor¬ 
tional  to  the  N-atom  density  in  the  afterglow.  Absorption  spectra  taken 
during  the  FTIR/long  path  cell  investigation  are  presented.  The  only  absorp¬ 
tion  features  found  in  the  spectra  are  due  to  SiH4  and  an  unidentified 
species  with  bands  at  900  and  2200  cm-1.  The  very  broadband  nature  of 
these  later  features  and  the  positions  of  the  band  centers  indicate  a 
complex  molecule  containing  Si  and  N  and  perhaps  some  H.  Finally,  the 
results  of  an  experiment  involving  a  modification  of  the  original  flowing 
afterglow  apparatus  to  look  at  SiH4  dissociation  induced  by  collisions 
with  N2(A3X)  molecules  are  presented.  In  this  experiment,  argon  metas¬ 
tables  are  created  in  a  dc  discharge  and  allowed  to  collide  with  N2 
molecules  downstream  of  the  discharge.  In  the  resulting  collisions  some 
of  the  argon  metastables  transfer  their  energy  to  N2  to  form  N2(A3£) 
metastables  which  can  then  react  further  downstream  with  SiH4.  If  the 
SiH4  is  dissociated  by  the  collisions  with  the  N2  metastables,  the  result¬ 
ing  radical  fragments  react  with  SiH4  to  form  higher  silanes  such  as 
Si2H6.  The  experiment  showed  conclusively  that  Si2H6  formation  was 
initiated  by  the  reaction  of  N2{A3I)  metastables  with  SiH4. 

Chapter  IV  presents  development  of  the  reaction  mechanism  and  the 
justification  for  the  proposed  scheme.  It  begins  with  a  general  discussion 
of  the  rate-equation  approach  and  discusses  the  bounds  on  the  general  bi- 
molecular  reaction  rate.  Following  this,  a  discussion  is  given  of  the 
general  characteristics  the  reaction  scheme  must  possess  in  order  to 
describe  the  observed  results.  A  general  cataloging  of  the  most  abundant 
species  found  in  the  afterglow  is  also  given.  From  the  fact  that  N-atoms 
are  the  most  abundant  reactive  species  present  in  the  afterglow  and  from 
the  general  behavior  of  the  reaction  scheme  presented  in  chapter  III  it  is 
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concluded  that  the  direct  reaction  of  N-atoms  with  SiH4,  if  finite,  must  be 
very  slow.  This  is  in  spite  of  the  results  of  chapter  III  that  indicate  that 
SiH4  dissociation  is  proportional  to  the  N-atom  density  in  the  afterglow. 
This  key  observation  indicates  that  in  any  proposed  reaction  scheme,  N- 
atoms  must  be  the  driver  for  SiH4  dissociation  but  cannot  react  directly 
with  the  molecule.  From  these  observations  and  the  overall  characteris¬ 
tics  the  reaction  scheme  must  possess,  a  general  model  of  the  reaction 
mechanism  is  developed.  This  general  model  attempts  to  describe  the  bulk 
of  the  observations  available  about  the  behavior  of  the  chemiluminescence 
and  the  SiH4  dissociation.  Arguments  are  then  presented  to  specifically 
define  the  species  involved  in  each  step  cf  the  general  model.  The  argu¬ 
ments  are  based  on  how  the  measured  data  can  be  explained  while  satisfy¬ 
ing  thermodynamic  constraints.  At  the  same  time,  consideration  is  also 
given  to  realistic  densities  of  species  available  in  the  afterglow,  and 
perturbations  that  the  proposed  scheme  would  have  on  these  densities.  The 
steps  in  the  general  model  include  the  initial  (slow)  dissociation  of  SiH4, 
reaction  of  N-atoms  with  SiHx  radicals,  formation  of  electronically 
excited  SiN  and  Si,  and  production  of  a  species  to  drive  the  later  (fast) 
SiH4  dissociation.  Vibrationally  excited,  ground  electronic  state  N2  is 
proposed  as  the  primary  driver  for  the  later  stages  of  the  SiH4  dis¬ 
sociation.  The  specific  set  of  chemical  reactions  proposed  to  describe  the 
reactions  of  SiH4  in  active  nitrogen  is  summarized  at  the  conclusion  of 
the  chapter. 

Details  of  the  coupled  rate-equation  model  are  presented  in  chapter 
V  and  the  results  of  the  model  compared  with  experimental  data.  The 
overall  model  includes  a  simplified  set  of  reactions  to  describe  the 
interaction  of  the  various  afterglow  species,  the  reactions  between  the 
afterglow  species  and  the  silicon  bearing  species,  and  a  complete  model 
of  the  interaction  of  the  vibrational  levels  of  ^(X1!).  Formulation  of 
these  subsets  of  the  overall  model,  the  selection  of  rate  constants  for  the 
reactions,  and  the  initial  conditions  used  as  input  for  the  model  are 
discussed.  When  available,  literature  values  for  the  rate  constants  are 
used.  When  not  available  in  the  literature,  rate  constants  are  constrained 
by  gas  kinetic  values.  The  vibrational  state  densities  are  treated  as  a  set 
of  coupled  equations  which  interact  through  the  exchange  of  vibrational 
quanta  (v-v)  and  vibration-to-translation  exchange  (v-T).  The  rates  for 
these  processes  are  taken  from  the  literature.  Following  this,  the  general 
behavior  of  the  model  is  presented  including  a  discussion  of  the  SiH4 
dissociation  and  the  behavior  of  the  chemiluminescence.  The  model  is 
shown  to  describe  the  general  behavior  of  the  chemiluminescence  ob- 
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served  in  this  work  and  the  time  dependence  of  the  N-atom  and  SiH4 
concentrations  measured  by  mass  spectrometer  in  reference  30.  The  model 
results  of  the  time  dependence  of  the  vibrational  populations  in  the  N2 
electronic  ground  state  are  discussed  and  the  chapter  concludes  with  a 
discussion  of  the  likelihood  for  the  formation  of  higher  nitrides  of  silicon 
in  the  gas  phase.  Chapter  VI  presents  some  conclusions  and  suggestions 
for  further  work. 

Background  and  Motivation  for  this  Work 

The  SiH4  molecule  possesses  tetrahedral  symmetry  with  an  Si-H 
bond  distance  of  1.46A.  Due  to  symmetry,  the  nine  vibrational  modes  of 
the  molecule  are  reduced  to  four  fundamental  modes.  There  is  one  non¬ 
degenerate  stretching  mode  Vi  at  2187cm-1,  one  doubly  degenerate  bending 
mode  v2  at  978cm-1,  and  two  triply  degenerate  bending  and  stretching 
modes  V3  at  2183cm-1  and  v4  at  910cm-1.  Formally,  only  the  V3  and  v4 
modes  are  infrared  active  but  due  to  the  small  differences  in  frequency 
between  the  vi  and  V3  modes  and  the  \>2  and  v4  modes,  vibration-rotation 
interaction  between  these  modes  results  in  all  four  modes  being  active.5 
The  electronic  states  of  SiH4  are  not  well  known.  The  molecule  possesses 
no  strong  absorption  bands  in  the  visible  or  ultraviolet.  Some  extremely 
weak  absorption  bands  seen  in  the  visible  have  been  attributed  to  Av=6-9 
overtones  of  the  ground  state  vibrational  levels.6  Absorption  spectra  of 
the  vacuum  ultraviolet  show  an  onset  at  =160nm  (7.8  eV)  followed  by 
continuous  absorption  to  shorter  wavelengths.7  Gordon8  has  performed  ab 
initio  calculations  of  the  energies  of  the  vertical  excited  states  of  SiH4 
and  these  calculations  place  the  first  excited  state  at  =9  eV  above  the 
ground  state.  Gordon,  however,  states  that  the  calculations  tend  to 
overestimate  the  excitation  energies  by  =1  eV. 

As  mentioned  at  the  beginning  of  this  chapter,  interest  in  SiH4  began 
in  this  laboratory  in  1983  while  studying  the  dissociation  of  SiH4  in  a  dc 
discharge  as  part  of  an  effort  to  study  the  plasma  chemical  reactions 
associated  with  PECVD  of  amorphous  silicon3.  In  that  work,  measurements 
of  the  rate  of  dissociation  of  SiH4  diluted  in  either  argon,  helium,  or 
nitrogen  was  measured  using  the  TDL  system,  mentioned  earlier,  to 
monitor  the  SiH4  density  downstream  of  the  discharge.  The  density  was 
determined  by  measuring  the  absorbance  of  an  isolated  rotation-vibration 
line  of  the  v4  band.  Premixed  gases  containing  1%-5%  SiH4  were  flowed 
through  the  discharge  under  controlled  mass  flow  conditions  so  the 
residence  time  in  the  discharge  was  well  known.  From  the  downstream 
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density  measurements,  the  fractional  dissociation  was  measured  for  a 
fixed  residence  time  and  a  fixed  discharge  current.  By  varying  the  SiH4 
residence  time  in  the  discharge  an  effective  "rate"  of  dissociation  could 
be  measured.  These  measurements  were  then  compared  with  calculated 
rates  due  to  electron  impact.  The  rates  were  calculated  by  solving  the 
Boltzmann  transport  equation  for  the  electrons  in  the  positive  column  to 
arrive  at  an  electron  energy  distribution  which  was  then  folded  with  a 
model  cross  section  for  SiH4  dissociation.  The  model  cross  section  was 
assumed  to  have  a  threshold  of  7.8  eV  corresponding  to  the  onset  of 
vacuum  u.v.  absorption.  For  He/SiH4  and  Ar/SiH4  mixtures,  the  measured 
rates  and  the  calculated  rates  agreed  within  a  factor  of  2.  Considering  the 
fact  that  no  attempt  was  made  in  the  calculations  to  include  the  effects 
of  the  sheath  regions  of  the  discharge,  this  was  considered  reasonable 
agreement.  From  these  results,  the  conclusion  was  drawn  that  for  SiH4 
diluted  in  these  rare  gases,  the  dissociation  rate  in  a  dc  discharge  was  set 
by  electron  impact9.  In  the  case  of  N2/SiH4  mixtures,  reasonable 
agreement  between  the  calculated  and  measured  rates  could  not  be 
achieved.  For  example,  in  the  case  of  1%  SiH4  in  N2,  the  measured  rate  was 
a  factor  of  70  larger  than  the  calculated  rate.  In  this  case  it  was  con¬ 
cluded  that  other  processes  besides  electron  impact  must  be  occurring  in 
the  discharge  to  lead  to  such  a  large  rate  of  dissociation10.  As  discussed 
in  the  beginning  of  this  chapter,  it  was  the  purpose  of  this  work  to 
determine  what  some  of  these  other  processes  might  be. 

Silane  Dissociation 

Probably  because  of  the  lack  of  detailed  knowledge  of  the  electronic 
states  of  SiH4,  the  mechanisms  by  which  the  molecule  is  dissociated  are 
not  well  known.  The  earliest  studied  dissociation  process  is  that  of 
pyrolysis.  These  early  studies  of  pyrolysis  lead  to  the  development  of  the 
first  method  of  amorphous  silicon  production  known  as  Chemical  Vapor 
Deposition  (CVD).  This  process  is  based  on  the  thermal  decomposition  of 
SiH4  at  atmospheric  pressure  on  a  hot  (typically  450-700°C)  substrate. 
This  is  in  contrast  with  PECVD  of  amorphous  silicon  which  is  done  at 
pressures  of  less  than  a  torr,  requires  a  (usually  RF)  gas  discharge,  and 
employs  substrate  temperatures  of  less  than  200°C.  There  have  been  a 
number  of  measurements  of  pyrolysis  of  SiH4  made  in  both  heated  vessels 
and  shock  tubes.  The  most  important  pyrolysis  work  using  a  heated  vessel 
is  probably  the  work  of  Purnell  and  Walsh11.  They  found  that  both  H2  and 
Si2H6  were  formed  at  a  pressure  dependent  rate  which  exhibited  an 
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activation  energy  of  Ea=55.9±0.5  Kcal/mole  (2.43  eV/molecule).  The 
pressure  dependence  found  by  Purnell  and  Walsh  was  given  as  P1-5.  Pres¬ 
sures  in  their  experiment  ranged  from  35-400  torr  and  temperatures 
ranged  from  650-700°C.  The  most  significant  shock  tube  measurements 
are  probably  those  of  Newman  et  al.12.  Under  the  shock  conditions  es¬ 
timated  at  PtotaMOOO  torr  and  T=1035-1 184°K  they  found  an  activation 
energy  of  Ea=52.7±1.4  Kcal/mole  (2.29  eV/molecule).  Recently  Coltrin  et 
al.13  have  published  a  unimolecular  decomposition  rate  for  SiH4,  based  on 
an  unpublished  RRKM  analysis  by  O'Neal  of  reference  12’s  original  data. 
Their  Arrhenius  expression  for  the  rate  is  given  as 

k(T.P)  *  (a+bP)ATPEXP(-Ea<RT)  (s-i)  (eq.  1-1) 

where  a=0.0504,  b=0.9496,  P  is  in  atmospheres,  A=2.54x1038,  p=-7.95,  and 
Ea=62.0  Kcal/mole  (2.69  eV/molecule).  Coltrin  et  al.13  claim  the  expres¬ 
sion  is  valid  over  the  pressure  range  of  5  torr  to  1  atm.  For  a  typical  CVD 
reactor  case  of  T=500°C  and  P=1  atm,  eq.  1-1  gives  a  unimolecular  decom¬ 
position  rate  of  k=8.46x10-3  s_1.  With  the  exception  of  a  1984  paper  by 
Robertson,  Hils,  and  Gallagher14,  who  claimed  that  SiH4  pyrolysis  was 
initiated  by  surface  decomposition,  the  pyrolysis  of  SiH4  is  considered  a 
unimolecular,  homogeneous  process.  Both  Purnell  and  Walsh11  and  New¬ 
man12  propose  the  decomposition  proceeds  via  the  reaction, 

SiH4  -  SiH2  +  H2  .  (eq.  1-2) 

From  the  heat  of  formation  data  given  in  appendix  B,  the  enthalpy  change 
in  the  above  reaction  is  AE=+59.9  Kcal/mole,  thus  the  activation  energy 
given  by  Coltrin  et  al.13  is  consistent  with  this  picture  of  thermal  dis¬ 
sociation. 

Dissociation  of  SiH4  in  discharges  can  occur  by  reactions  with  ions. 
For  example,  reactions  of  SiH4  with  rare-gas  ions  can  result  in  dissocia¬ 
tion  by  reactions  such  as 


SiH4  +  Ar  -*SiH3  +  H  +  Ar  (eq.  1-3) 

\ 

Since  Ar  is  often  used  as  a  buffer  gas  in  amorphous  silicon  deposition, 
reactions  such  as  eq.  1-3  could  be  important  in  governing  the  overall  SiH4 
dissociation  under  discharge  conditions.  Ion-molecule  reactions  involving 
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silane  and  ionized  fragments  from  the  break  up  of  silane  can  also  be 
important.  Reactions  such  as 


SiH4  +  Hn  -*  SiH3  +  H2  +  (eq.  1-4) 

and 

SiH4  +  SiHm -*  SiH3  +  SiH1  ,m  (eq.  1.5) 

where  n=2  or  3  and  m=0,1,2,  can  increase  the  rate  of  silane  dissociation 
under  discharge  conditions.  Haaland15  has  recently  discussed  a  number  of 
ion-molecule  reactions  that  can  occur  in  silane  discharges.  Haaland’s  work 
includes  measurements  of  reaction  rates  of  SiHU  with  the  noble-gas  ions 
He+,  Ne+,  Ar+,  Kr+,  and  Xe+  and  a  number  of  ion-molecule  reactions  that  can 
occur  among  dissociated  fragments.  Also,  Chatham,  et  al.16  measured 
reactions  rates  for  SiH4  with  Ar+,  He+  and  Ne+  while  Turbin,  et  al.17  and 
Haller18  have  looked  at  ion-molecule  reactions  among  silane  fragments.  It 
is  assumed  in  the  work  presented  here  that  ion  densities  in  the  afterglow 
are  sufficiently  small  that  ion-molecule  reactions  of  this  type  can  be 
neglected. 

Likewise  electron  impact  dissociation  of  SiH4  will  not  be  considered 
in  this  work,  but  under  discharge  conditions  this  process  is  significant. 
Measurements  of  electron  impact  dissociation  of  SiH4  have  been  made 
under  discharge  conditions  and  with  low  energy  electron  beams.  Nolet19 
measured  the  rate  of  disappearance  of  SH4  diluted  in  argon  in  a  flowing- 
gas  D.C.  discharge,  while  Perrin  and  Aarts20  measured  dissociative  excita¬ 
tion  of  electronically  excited  fragments  of  SiH4  using  a  low  energy 
electron  beam.  Perrin,  et  al.21  used  a  multipole  reactor  and  a  mass 
spectrometer  to  measure  the  rate  of  disappearance  of  SiH4  in  a  flowing- 
gas  system.  The  narrow  energy  spread  of  the  electrons  in  the  multipole 
reactor  allowed  Perrin  et  al.21  to  measure  a  total  dissociation  cross 
section  as  a  function  of  electron  energy  from  *17  eV  to  110  eV.  As 
mentioned  earlier,  DeJoseph,  et  al.3  used  a  method  similar  to  Nolet19  to 
measure  the  rate  of  silane  dissociation  under  D.C.  discharge  conditions  and 
used  a  modelling  approach  to  arrive  at  a  dissociation  cross  section.  The 
model  cross  section  was  assumed  to  have  a  threshold  coincident  with  the 
vacuum  ultraviolet  absorption  threshold.  This  threshold  is  also  very  close 
to  the  *8  eV  threshold  arrived  at  by  extrapolation  of  the  data  of  Perrin,  et 
al.21  It  thus  seems  reasonable  that  the  threshold  for  electron  impact 
dissociation  of  SiH4  is  very  close  to,  if  not  coincident  with,  the  threshold 
for  vacuum  ultraviolet  absorption  of  7.8  eV.  The  products  formed  from  the 
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electron  impact  dissociation  are  dependent  on  the  electron  energy,  but  in 
the  case  of  vacuum  u.v.  absorption  at  147  nm  (8.4  eV),  Perkins,  et  al.22 
suggest  the  products  are  primarily  SiH2  (83%)  and  SiH3  (17%).  However  the 
product  distribution  from  photodissociation  may  differ  from  electron 
impact  dissociation  since  the  electrons  can  induce  transitions  which  are 
not  necessarily  optically  allowed. 

Heavy  particle  dissociation  of  silane  can  be  divided  into  two  types 
of  collisions.  In  a  reactive  collision,  the  colliding  species  chemically 
bonds  to  one  of  the  atoms  of  the  silane  molecule.  The  hydrogen  abstraction 
reaction  of  atomic  hydrogen  with  silane  is  of  this  type: 

SiH4  +  H  -*  SiH3  +  H2.  (eq.  1-6) 

The  second  type  of  collision  to  consider  is  a  super  elastic  collision  in 
which  the  internal  energy  of  the  colliding  species  is  sufficient  to  dis¬ 
sociate  silane,  but  the  reacting  species  maintains  its  identity  (but  with 
less  internal  energy)  after  the  collision.  The  reaction  of  Hg(3Pi)  with 
silane  is  such  a  reaction: 

SiH4  +  Hg(3Pi)  -  SiH3  +  H  +  HgOSo).  (eq.  1-7) 

Reaction  1-6  has  been  studied  by  a  number  of  authors.  Kushner23,  has 
recently  reviewed  the  rate  measurements  and  gives  a  rate  constant  of 
2.7x1 0-12  cm3/s.  Reaction  1-7  is  known  from  experiments  with  mercury 
sensitization  of  silane.  In  these  experiments,  a  gas  mixture  of  SiH4  and  Hg 
is  irradiated  with  257  nm  (4.8  eV)  resonance  radiation  from  a  Hg  lamp 
resulting  in  SiH4  decomposition  and  hydrogen  evolution.  While  this  reac¬ 
tion  has  been  known  since  the  observations  of  Emeleus  and  Stewart24  in 
1936,  it  has  only  been  recently  that  details  of  the  reaction  have  been 
studied25.  In  particular,  it  has  not  been  clear  whether  the  reaction  in¬ 
volved  the  Hg(3Pi)  state,  which  is  excited  by  the  257  nm  radiation,  or 
whether  the  reaction  involved  the  metastable  Hg(3P0)  state  formed  in  the 
reaction 

Hg(3Pi)  +  M  -  Hg(3P0)  +  M  (eq.  1-8) 

where  M  is  any  gas  species.  The  two  states  differ  in  energy  by  only  0.2 
eV26,  thus  reaction  1-8  might  be  expected  to  be  reasonably  fast.  Perrin 
and  Allain25  measured  the  rate  of  quenching  of  both  the  metastable 
Hg(3P0)  state  and  the  Hg(3Pi)  state  by  silane  and  report  quenching  rates  of 
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7.3±1 .3x10-1°  and  8.3±0.7x10_1°  cm3/s  respectively  for  the  two  states. 
Furthermore,  they  suggest  that  the  quenching  of  both  states  is  primarily 
by  electronic-to-electronic  transfer  to  SiH4,  followed  by  dissociation  of 
SiH4.  In  particular,  they  rule  out  the  reactive  collision 

SiH4  +  Hg^P,)  -  SiH3  +  HgH  (eq.  1-9) 

and  suggest  that  for  M=SiH4  reaction  1-8  is  very  slow.  This  is  in  spite  of 
the  lack  of  any  known  electronic  states  of  SiH4  below  7.8  eV.  This  points 
up  a  difficulty  with  proposing  reaction  schemes  involving  silane.  From  the 
energetics  alone,  one  would  be  tempted  to  conclude  that  reaction  1-9 
would  be  more  likely  than  reaction  1-7.  If  one  were  then  to  consider  the 
lack  of  excited  states  of  SiH4  below  7.8  eV,  one  would  most  surely  rule 
out  reaction  1-7.  The  difficulty  with  this  reasoning  lies  in  the  details  of 
the  collisions.  Because  of  the  thermal  energies  involved,  the  colliding 
particles  can  spend  a  relatively  long  time  in  close  proximity  and  thus 
there  is  some  probability  that  a  collision  complex  of  the  colliding  pairs 
will  be  formed.  The  subsequent  breakup  of  this  unstable  complex  may  lead 
to  the  formation  of  species  that  might  not  be  expected  to  be  formed  in  a 
brief,  or  direct  collision27.  Reactions  of  heavy  particles  with  SiH4  must 
therefore  be  considered  a  difficult  area  to  predict  a  priori  rates  of 
reaction  and  details  of  the  reaction  channels. 


Active  nitrogen  has  been  studied  for  more  than  100  years  and  the 
reader  is  referred  to  the  book  by  Wright  and  Winkler  for  a  good  historical 
perspective  up  until  19684.  While  extensively  studied  (Wright  and  Winkler 
list  1529  references)  active  nitrogen  remains  for  the  most  part  a  poorly 
understood  subject.  Wright  and  Winkler  summarize  the  difficulties 
associated  with  this  enigma  in  the  following  quote28: 

"Active  nitrogen  and  its  behavior  stand  revealed  as  more 
complex  than  generally  recognized.  Moreover,  apparent 
contradictions  in  observations  and  interpretations  are 
indicated  as  reflecting  a  greater  sensitivity  of  the  systems 
to  experimental  conditions  than  has  generally  been 
appreciated.  As  a  consequence,  assessment  of  much  of  the 
available  data,  at  the  present  time,  is  likely  to  be  arbitrary" 


Unfortunately,  there  has  been  little  done  since  1968  to  change  the  above 
assessment. 

One  of  the  earliest  and  most  performed  experiments  with  active 
nitrogen  involved  the  introduction  of  various  compounds  into  the 
afterglow,  studying  the  resulting  changes  in  the  afterglow  (if  any), 
studying  the  spectroscopy  of  any  emission  from  the  reaction  zone,  and 
studying  any  new  compounds  formed  during  the  reaction.  Thus  while  the 
original  idea  of  introducing  the  SiH4  into  active  nitrogen  seemed 
straightforward,  it  fell  in  with  a  group  of  experiments  which  have  been 
extensively  studied  for  many  years  with  many  conflicting  results  and 
interpretations.  Thus  it  was  surprising  when  this  work  was  begun  in  1984 
that  only  one  paper,  by  Dewhurst  and  Cooper  in  I96029,  could  be  found  of 
any  study  of  SiH4  in  active  nitrogen.  Since  this  work  was  begun,  two  more 
studies  of  SiH4  in  active  nitrogen  have  appeared,  both  in  1986.  Horie,  et 
al.30  looked  at  SiH4  in  active  nitrogen  in  a  flowing  afterglow  similar  to 
the  one  used  in  this  work.  Piper  and  Caledonia31  also  used  a  flowing 
afterglow  experiment  but  with  higher  flow  speeds  than  used  in  this  work 
or  reference  30. 

Dewhurst  and  Cooper29  reacted  SiH4  in  a  glass  bulb  containing  active 
nitrogen  at  a  pressure  of  1.2  torr.  A  liquid  nitrogen  trap  was  positioned 
downstream  of  the  bulb  to  trap  condensables  from  the  reaction.  Experi¬ 
ments  were  performed  at  room  temperature  and  250°K.  They  described  a 
"lilac  colored  reaction  flame”  at  both  temperatures  and  reported  "an 
intense  blue  glow"  when  the  silane  flow  was  stopped  and  the  active 
nitrogen  allowed  to  reach  the  downstream  LN2  trap.  They  also  reported  "a 
preliminary  examination  of  the  emission  spectra"  in  which  they  identified 
atomic  silicon  emission.  Reaction  products  from  the  cold  trap  were  later 
analyzed  using  infrared  absorption  and  gas  chromatography.  No  ammonia, 
hydrazine,  or  other  condensable  products  could  be  identified  in  either  the 
room  temperature  or  250°K  data. 

Horie  et  al.30  used  a  small  flowing  afterglow  system  similar  to  the 
one  to  be  described  in  this  work.  Pressures  in  their  experiment  were  kept 
near  3  torr.  They  published  spectra  of  the  reaction  flame  over  the  range 
450-800nm  and  discussed  some  observed  features  in  the  range  of  200- 
450nm.  Spectra  were  not  corrected  for  system  response.  Emission  spectra 
from  Si*,  SiN*,  and  an  unidentified  continuum  from  450  to  810nm  were 
discussed.  Changes  in  the  continuum  when  deuterated  silane  (SiD4)  was 
substituted  were  also  discussed.  They  specifically  stated  that  no  emis¬ 
sions  from  SiH,  Si2,  NH,  or  NH2  could  be  identified.  In  addition  to  spectra 
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of  the  most  intense  region  of  the  flame,  spatially  resolved  emission  was 
presented  for  Si*,  SiN*,  N2(1st  positive),  and  at  two  wavelengths  in  the 
continuum.  Mass  spectroscopy  was  used  to  measure  the  density  of  N- 
atoms,  SiH4,  and  H2  as  a  function  of  time  in  the  reaction  zone.  The  mass 
spectrometer  was  also  used  to  search  for  reaction  products.  Only  H2  could 
be  found.  One  very  interesting  observation  by  Horie  et  al.  was  the 
measurement  of  chemiionization  in  the  reaction  flame.  Attempts  to 
determine  the  mass  of  the  ions  were  unsuccessful,  but  they  claimed  the 
magnitude  of  the  ion  current  was  comparable  to  that  due  to  NO+  when  NO 
was  introduced  into  the  afterglow.  Lastly,  Horie,  et  al.  proposed  a  sketchy 
reaction  scheme  to  explain  their  observed  results.  More  details  of  their 
proposed  scheme  will  be  discussed  later,  but  essentially  a  chain  reaction 
was  proposed  triggered  by  the  reaction  of  N2(A3L)  metastables  with  SiH4 
resulting  in  SiH4  dissociation.  In  their  proposed  scheme  N2(A3Z)  metas¬ 
tables  were  produced  by  exothermic  reactions  which  further  reacted  with 
SiH4.  They  also  invoked  N2(A3X)  to  explain  the  Si*  emission  by  collisions 
of  the  second  kind  with  Si. 

Piper  and  Caledonia31  also  measured  emission  spectra  of  the 
chemiluminescence  in  a  flowing  afterglow  system.  Pressure  in  their 
experiment  ranged  from  1.5  to  2  torr.  While  there  is  no  detailed  data  to 
support  this,  it  appears  from  the  spectra  presented  in  references  30  and 
31  that  Piper  and  Caledonia  had  higher  sensitivity  in  their  spectrometry. 
The  spectral  features  identified  by  Piper  and  Caledonia  differed  in  both 
species  identified  and  relative  intensity  from  those  of  Horie,  et  al.  Like 
Horie  et  al.,  Piper  and  Caledonia  identified  a  number  of  lines  from  Si*  and 
identified  some  SiN*  emission.  On  the  other  hand,  while  the  SiN*  emission 
was  a  dominant  feature  in  the  spectra  of  Horie  et  al.,  the  emission  was 
very  weak  in  the  spectra  of  Piper  and  Caledonia.  Unlike  Horie,  et  al.,  Piper 
and  Caledonia  observed  emission  from  both  SiH  and  NH,  the  NH  emission 
being  one  of  the  strongest  features  in  their  spectra.  Piper  and  Caledonia 
also  observed  emission  from  CN,  which  was  not  mentioned  by  Horie,  et  al. 
There  is  also  disagreement  between  the  two  on  the  observed  continuum 
emission.  Piper  and  Caledonia  reported  a  continuum  between  280  and 
430nm,  with  no  emission  above  430nm.  As  previously  discussed,  Horie,  et 
al.  observed  continuum  emission  from  450  to  810nm  with  no  emission 
seen  below  450nm.  Piper  and  Caledonia  also  discuss  spectra  taken  under 
different  afterglow  conditions.  By  changing  the  configuration  of  their 
experiment,  they  were  able  to  produce  active  nitrogen  that  only  contained 
N2(A3I),  and  also  varying  densities  of  N-atoms.  The  results  of  these 
measurements  were  that  no  emission  was  seen  with  only  N2(A3I)  present 
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and  the  intensity  of  the  emission  was  stronger  as  the  N-atom  density  was 
increased.  Finally,  Piper  and  Caledonia  reported  reaction  rates  for  a 
number  of  afterglow  species  to  react  with  SiH4.  These  rates  were  deter¬ 
mined  by  measuring  the  rate  of  loss  of  a  particular  species  due  to  reac¬ 
tions  with  SiH4.  The  table  below  summarizes  their  results. 


Table  1-1:  Rate  Coefficients  for  Various  Nitrogen  Species  with  Silane. 

Data  from  Reference  31. 


Nitrogen  Species  Rate  Coefficient  nO'^cm3/^ 


N2(A3I,  v=0) 
N2(A3L,  v=1) 
N2(a’1I,  v=0) 
N(2D) 

N(2P) 

N(4S) 


0.59±0.1 1 
0.78±0.1 7 
14.5±3.5 
4.0+0.9 
0.06+0.02 
<0.006 


The  reaction  rate  for  N(4S)  was  set  by  the  smallest  rate  they  could 
measure  in  their  experiment.  No  reaction  of  N(4S)  with  silane  could  be 
observed.  Besides  the  three  references  given  here,  one  might  be  tempted 
to  compare  the  reaction  of  SiH4  in  the  afterglow  to  that  of  CH4.  Wright  and 
Winkler4  give  a  number  of  references  for  CH4  in  active  nitrogen.  While  the 
results  of  the  investigations  (emission  spectra,  species  formed)  are 
similar,  interpretation  of  the  reaction  mechanism  is  an  open  question. 
Thus  one  is  in  much  surer  waters  by  staying  with  the  problem  at  hand.  In 
the  following  pages,  this  problem  will  be  presented. 
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Chapter  II 

Description  of  the  Experiment 


Most  flowing  afterglow  experiments  have  used  either  a  "condensed" 
dc  discharge  or  a  microwave  discharge  to  produce  the  afterglow  (see 
reference  4).  The  condensed  discharge  is  produced  by  placing  a  capacitor 
(formerly  a  condenser)  in  parallel  with  the  discharge  electrodes  to 
produce  an  oscillating  voltage  discharge  which  yields  a  high  concentration 
of  nitrogen  atoms  in  the  afterglow.  The  microwave  discharge  is  typically 
made  up  of  a  small  microwave  cavity  driven  at  a  few  GHz.  For  this  inves¬ 
tigation.  a  conventional  dc  discharge  with  current  regulation  was  used. 

The  reason  for  using  the  dc  discharge  was  twofold.  First,  since  it  was 
desired  to  measure  the  dissociation  of  SiH4  as  a  function  of  discharge 
parameters,  the  dc  discharge  offered  a  more  easily  measurable  and 
reproducible  set  of  parameters  (current  and  voltage)  than  the  microwave 
discharge.  Second,  the  positive  column  of  the  dc  discharge  appeared  to  be 
simpler  (but  not  simple)  system  for  modeling  the  afterglow  and  hence 
predicting  what  species  and  their  concentrations  would  be  present.  As 
will  be  discussed  later,  modeling  even  the  dc  discharge  proved  to  be  so 
difficult  that  this  latter  advantage  failed  to  materialize.  Figure  2-1  is  a 
schematic  of  the  experiment.  The  basic  experiment  consisted  of  a  flowing 
afterglow  with  associated  flow  and  pressure  controls,  an  absorption  cell 
and  Tunable  Diode  Laser  (TDL)  system  for  measuring  the  silane  density 
downstream  of  the  afterglow,  and  a  visible/ultraviolet  monochromator  for 
recording  the  spectrum  of  the  radiation  from  the  reaction  zone.  Details  of 
the  system  will  be  described  below. 

Description  of  the  Experiment-  Flowing  Afterglow  System 

Pure  nitrogen  (Matheson  ultrahigh  purity,  99.999%  min.)  was  flowed 
through  an  8  mm  i.d.  pyrex  discharge  tube  with  electrodes  out  of  the  flow 
as  shown  in  the  figure.  The  spacing  between  the  side  arms  housing  the 
electrodes  was  25  cm.  The  electrode  material  was  Kovar.  The  power 
supply  was  a  current  regulated  dc  supply  capable  of  6  kV  maximum 
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2-1:  Schematic  of  the  Experiment 


voltage.  A  Weston  model  911  analogue  milliammeter  was  used  to  measure 
discharge  current.  Discharge  voltage  measurements  used  a  Fluke  model 
80K-40  probe  (1000  Mohm  input  resistance)  and  a  Fluke  model  8050A  4 
1/2  digit  DMM  were  used.  The  nitrogen  mass  flow  was  controlled  by  an 
MKS  model  254  Mass  Flow  Controller.  This  unit  would  maintain  a  preset 
mass  flow  of  up  to  three  gases  by  measuring  the  mass  flow  of  each  gas 
;th  an  MKS  model  258  flow  meter  and  via  a  feedback  loop  control  the 
rlow  with  an  MKS  model  248A  flow  control  valve.  The  mass  flow  meters 
were  of  the  thermal  conductivity  type  and  the  model  254  allowed  gas- 
correction  factors  to  be  entered  via  front  panel  switches  so  that  true 
mass  Now  could  be  read  out  and  controlled  for  any  gas.  This  system  also 
controlled  the  silane/helium  mixture  which  was  introduced  into  the 
afterglow.  A  mixture  of  0.9±0.1%  silane  in  helium  was  used  throughout  the 
experiment  so  that  small  amounts  of  silane  could  be  titrated  using  the 
mass  flow  meters.  Flow  rates  of  pure  silane  would  have  been  too  small  to 
be  accurately  measured  and  controlled.  Since  its  lowest  energy  level  is  at 
19.8  eV,  helium  is  assumed  to  not  interact  with  the  afterglow  or  silane  in 
the  experiment.  Pressure  in  the  system  was  measured  at  two  points  using 
MKS  capacitance  manometers.  A  model  315BHS  (10  torr  head)  was  used  in 
the  afterglow  region  and  a  model  220BHS  (10  torr  head)  was  used  to 
measure  the  pressure  in  the  discharge  region.  The  pressure  signal  from  the 
discharge  region  was  fed  to  an  MKS  model  252A  Exhaust  Valve  Controller 
which  in  turn  controlled  an  MKS  model  253A  Throttling  Valve.  This  system 
thus  allowed  independent  control  and  good  stability  of  both  the  mass  flow 
and  the  pressure.  It  should  be  pointed  out  that  over  the  pressures  and 
flows  of  the  experiments  run,  the  pressures  recorded  by  the  two 
capacitance  manometers  never  differed  by  more  than  a  few  percent. 

Nitrogen  flowing  out  of  the  discharge  then  entered  a  36  mm  i.d. 
pyrex  section  where  the  straw-yellow  Lewis-Rayleigh  afterglow4  could  be 
seen.  This  larger  i.d.  section  extended  about  45  cm  from  the  end  of  the 
discharge  tube.  The  silane/helium  mixture  could  be  introduced  into  the 
afterglow  at  a  point  25  cm  from  the  end  of  the  discharge  tube.  During  the 
course  of  the  investigation,  the  silane  injection  tube  was  moved  to  a  point 
at  the  edge  of  the  discharge  tube,  but  most  of  the  data  was  recorded  with 
the  tube  at  the  25  cm  position.  A  2.5  cm  quartz  window  was  cemented  to 
the  side  of  the  36  mm  section  at  the  silane  injection  port  to  allow  recor¬ 
ding  of  emission  spectra  down  to  220  nm  in  the  u.v.  During  the  course  of 
the  experiment,  an  identical  flow  tube  was  constructed  entirely  of  quartz 
to  allow  less  restricted  viewing,  but  most  of  the  spectra  that  will  be 
presented  here  were  taken  with  the  quartz  window.  The  visible/u.v. 


spectrometer  will  be  described  later.  Gas  flowing  out  of  the  36  mm  i.d. 
section  then  entered  a  1.2  m  absorption  cell  fitted  with  a  barium  fluoride 
window  on  one  end  and  a  mirror  on  the  other.  This  allowed  the  gases 
leaving  the  reaction  zone  to  be  probed  with  the  Tunable  Diode  Laser  (TDL) 
system  which  is  described  below.  By  monitoring  the  infrared  absorption  of 
silane,  its  density  could  be  measured  downstream  of  the  reaction  zone. 

The  gas  then  left  the  absorption  cell  and  was  pumped  by  a  Welch  model 
1397  mechanical  pump  (425  liters/s  @  1  torr)  equipped  with  an  oil  trap. 
Due  to  the  low  concentrations  of  silane  used  in  the  experiment,  no  special 
precautions  were  taken  with  the  pump  except  to  vent  the  exhaust  to  the 
outside.  A  blower  with  an  enclosed  motor  was  used  to  draw  the  exhaust 
gases  from  the  pump. 

The  design  of  the  experiment  was  based  to  some  extent  on  ease  of 
fabrication  and,  unfortunately,  characterization  of  the  flow  velocity  was 
difficult.  The  flow  velocity  of  gas  in  a  cylindrical  tube  is  given  by32 

v-Q  .  - 2  7  Q(sccm)  (cm/sec) 

PA  P(torr)-A(cm2)  (eq.  2-1) 

where  Q  is  mass  flow,  P  is  pressure,  and  A  is  the  cross-sectional  area  of 
the  tube.  While  mass  flow  and  pressure  could  be  accurately  measured,  the 
abrupt  change  in  tube  diameter  from  the  8  mm  i.d.  discharge  section  to  the 
36  mm  i.d.  reactor  section  meant  that  the  flow  velocity  dropped  by  the 
ratios  of  the  areas  or  (8/36)2=1/20.  Thus  there  was  a  transition  region  in 
which  the  flow  velocity  changed  by  approximately  a  factor  of  20.  At  low 
pressures  (below  3  torr)  the  flow  velocity  was  limited  by  the  conductance 
of  the  lines  to  the  pump  so  that  at  1  torr  the  maximum  mass  flow  of  N2 
was  limited  to  40  seem.  At  higher  pressures,  the  mass  flow  measure¬ 
ments  were  constrained  to  the  limit  of  the  mass  flow  meter  which  was 
400  seem.  The  maximum  flow  velocity  in  the  small  diameter  section 
ranged  from  =1000  cm/s  at  1  torr  to  =3000  cm/s  at  3  torr  to  =1000  cm/s 
at  10  torr.  From  observations  of  the  straw-yellow  afterglow,  it  was 
apparent  that  at  high  pressure  (8-10  torr)  and  high  mass  flow,  the  transi¬ 
tion  region  extended  beyond  the  point  of  introduction  of  the  silane/helium 
mixture.  Thus  it  was  not  possible  to  accurately  know  the  flow  speed  in  the 
reaction  zone  at  all  values  of  mass  flow  and  pressure.  When  measurements 
were  to  be  compared  at  different  pressures,  the  quantity  Q/P  was  held 
constant. 
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A  schematic  of  the  TDL  system  is  shown  in  figure  2-2.  The  system  is 
a  commercial  system  built  by  Laser  Analytics,  Inc.  The  system  centers 
around  a  small  (=0.025  mm3),  low  power,  continuous  wave  solid-state 
semiconductor  laser  composed  of  lead  salts  such  as  PbS,  PbSe,  PbTe,  and 
PbSnTe33.  These  lasers  emit  in  the  wavelength  range  2-30  microns  in  the 
infrared  with  the  wavelength  range  of  a  particular  laser  fixed  by  its 
composition.  A  single  laser  may  only  cover  0.05  micron,  hence  many  lasers 
are  needed  to  cover  the  total  range  of  the  system.  The  advantages  of  these 
lasers  over  other  infrared  techniques  are  their  very  narrow  linewidth 
(=3x1 0-4  cm-1)  and  high  brightness  compared  to  a  standard  blackbody 
source  at  1500°K.  The  narrow  linewidth  allows  the  true  Doppler  profile  of 
an  absorption  line  at  low  pressure  to  be  seen  which  results  in  high  sen¬ 
sitivity  for  detecting  species  by  absorption.  In  addition,  the  narrow 
linewidth  allows  the  true  profile  of  the  line  to  be  measured  and  hence  no 
correction  for  instrument  resolution  is  needed.  The  high  brightness  of  the 
lasers  allow  rapid  data  acquisition  using  fast  solid  state  infrared  detec¬ 
tors.  A  major  disadvantage  is  that  the  lasers  operate  at  temperatures 
below  that  of  liquid  nitrogen,  hence  either  liquid  helium  or  a  helium 
refrigeration  system  must  be  used  to  cool  the  lasers.  This  system  used  a 
closed-cycle  helium  refrigeration  system  which  is  indicated  in  figure  2-2 
by  "Lines  to  Compressor".  Another  disadvantage  (although  minor)  is  that 
the  lasers  usually  run  multimode  so  that  some  form  of  mode  selection  is 
necessary  to  achieve  marrow  linewidth  operation.  This  is  accomplished  in 
this  system  by  passing  the  laser  through  the  half-meter  monochromator 
shown  in  figure  2-2.  The  wavelength  of  a  particular  composition  laser  is 
varied  by  changing  the  temperature  of  the  diode12.  In  practice,  this  is 
accomplished  in  two  ways.  First,  a  wavelength  range  of  interest  is 
selected  by  changing  the  temperature  of  the  heat-sink  on  which  the  diode 
is  mounted  using  a  small  resistive  heater  with  a  feedback  circuit  to 
control  the  heat-sink  temperature.  Once  the  heat-sink  temperature  has 
stabilized,  the  diode  temperature  is  varied  about  this  point  by  changing 
the  current  through  the  diode.  The  current  can  either  be  swept  rapidly  and 
the  laser  output  observed  with  an  oscilloscope  or  the  current  can  be  swept 
slowly  while  the  laser  beam  is  mechanically  chopped  and  phase  sensitive¬ 
ly  detected.  All  of  the  data  taken  in  this  work  used  this  latter  mode 
operation.  The  laser  beam  was  chopped  with  a  tuning  fork  chopper  at  400 
Hz  (figure  2-2)  and  the  beam  detected  with  a  liquid  nitrogen  cooled 
HgCdTe  detector  with  matched  preamp.  The  preamp  output  was  fed  to 
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a  PAR  model  HR-8  lockin  amplifier.  The  output  from  the  lockin  was  then 
digitized  and  sent  to  the  H.P.  9836C  computer  (see  figure  2-1).  The  diode 
current  was  swept  by  the  computer  using  a  digital-to-analogue  converter 
(DAC)  to  drive  an  external  current  drive  on  the  laser  system.  For  the 
experiments  here,  the  system  was  operated  in  the  'external  beam  mode' 
with  the  two  beam  turning  mirrors  removed  from  the  sample  compartment 
(see  figure  2-2).  Thus  the  beam  was  directed  out  of  the  instrument 
through  the  lower  sample  compartment  port  (figure  2-2),  through  the 
absorption  cell,  and  then  back  into  the  instrument  through  the  other 
sample  compartment  port.  To  measure  the  silane  density  in  the  absorption 
cell,  a  pair  of  PbSnTe  diodes  were  chosen  which  covered  the  wavelength 
range  of  10.5  to  11.7  microns.  These  lasers  were  mounted  in  the  system 
cold  head  and  either  one  could  be  selected  by  physically  shifting  the  cold 
head  in  its  mounts.  With  these  diodes,  the  v4  band  of  silane  could  be 
scanned.  This  band  was  chosen  for  monitoring  the  silane  density  primarily 
because  of  the  availability  of  stable  diodes  in  this  wavelength  range 
which  exhibited  good  tuning  characteristics  with  little  ’mode  hopping’ 

(see  reference  33). 

Calibration  of  the  TDL  System 

Before  silane  density  measurements  could  be  made,  the  system  had 
to  be  calibrated.  The  basic  principle  behind  the  density  measurement  is 
straightforward;  the  density  is  determined  by  measuring  the  absorption  of 
the  laser  radiation  by  an  isolated  absorption  line  of  the  gas.  The  intensity 
of  the  laser  radiation  transmitted  through  an  absorption  cell  of  length  I 
which  contains  a  gas  of  density  n  is  given  by  (neglecting  window  losses) 
Beer's  law34  as 

l(v)  =  1 0(v)-Exp[-n-l-k(v)]  (eq.  2-2) 

where  lo(v)  is  the  intensity  of  the  radiation  on  the  cell,  v  is  the 
wavenumber  of  the  radiation,  and  k(v)  is  the  absorption  coefficient  of  the 
gas.  Figure  2-3  shows  a  transmittance  spectrum  of  the  v4  band  of  SiH4 
taken  at  1  cm-1  resolution*  with  a  Fourier  Transform  Infrared  (FT-IR) 
spectrometer  by  the  author.  The  band  is  overlapped  by  the  much  weaker  V2 


For  those  not  familiar  with  units  of  cm-1 ,  it  is  the  reciprocal  of  the 
wavelength  expressed  in  cm  (1  cm-1  =104/XmiCrons)- 
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Figure  2-3:  Silane  Absorption  Near  900  cm*1  (1  cm'1  Resolution) 


band  which  is  centered  approximately  59  cm-1  above  the  V4  band.3  The 
transmittance  is  defined  by 


T (v)  =  l(v)/l0(v)  (eq.  2-  3) 

Thus  the  transmittance  spectrum  is  directly  related  to  the  absorption 
coefficient.  As  can  be  seen,  the  absorption  coefficient  is  an  extremely 
complex  function  of  wavelength  (or  wavenumber).  If  the  resolution  is 
improved  to  0.06  cm-1  as  shown  in  figure  2-4  (FT-IR  in  this  laboratory), 
still  more  features  are  revealed.  Finally,  if  one  scans  the  region  at  edge  of 
the  Q-branch  near  910.7  cm-1  with  the  diode  laser  as  in  figure  2-5, 
individual  absorption  lines  can  be  seen.  The  lower  trace  in  figure  2-5  is  an 
etalon  fringe  pattern  made  by  inserting  a  1  inch  germanium  etalon  (see 
figure  2-2)  into  the  TDL  beam  and  repeating  the  scan.  The  etalon  scan 
gives  a  means  of  calibrating  the  wavelength  range  of  the  TDL  scan.  In 
figure  2-5,  the  spacing  of  the  etalon  peaks  is  0.04855cm-1.  It  should  be 
pointed  out  that  figure  2-5  is  nfli  a  transmittance  spectrum  but  is  simply 
a  plot  of  l(v)  vs.  v  from  the  diode  laser.  The  density  of  silane  was  deter¬ 
mined  by  measuring  the  absorption  of  one  of  these  isolated  lines.  For  an 
isolated  line  with  both  Doppler  and  pressure  broadening,  the  absorption 
coefficient  is  described  by  the  Voigt  function34: 

k(v)=S/yD-[ln(2)/jt]  /2-V(x,y)  (eq.  2-4) 

where  S  is  the  line  intensity,  Yd  is  the  Doppler  half-width  at  half  maxi¬ 
mum,  and  V(x,y)  is  the  Voigt  function  given  by 


V(x,y)=— J  Exp(-t2)/[y2+(x-t)2]dt 


(eq.  2-5) 


with  x  and  y  given  by 
x=[ln(2)]1/%c/YD 


(eq.  2-6) 
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Figure  2-4:  Silane  Absorption  Near  900  cm'1  (0.06  cm1  Resolution) 


(eq.  2-7) 


y=[ln(2)]1  /2(v-Vo)/yD 

where  YoandYc  are  the  Doppler  half-width  at  half  maximum  and  the 
pressure  broadened  half-width  at  half  maximum  for  the  line,  and  v0  is  the 
line  center.  Yc  is  given  by 

N 

Tc-X^-Pi 

i=i  (eq.  2-8) 

where  T\  is  the  broadening  coefficient  of  the  ith  species  in  the  gas  mix¬ 
ture  and  Pj  is  the  partial  pressure  of  the  ith  species.  Thus  from  this 
equation  it  is  clear  that  the  absorption  coefficient  k(v)  in  equation  2-2 
depends  on  the  pressure  and  composition  of  the  surrounding  gas.  This 
means  that  under  conditions  where  the  gas  composition  or  the  total 
pressure  vary  significantly,  a  measurement  of  the  peak  absorption  of  an 
isolated  line  is  not  sufficient  to  determine  the  density  of  the  absorbing 
species.  On  the  other  hand,  the  integral  of  k(v)  over  the  full  extent  of  the 
line,  given  by 


i.e.  the  line  intensity  is  a  function  of  the  gas  temperature  but  independent 
of  pressure.  Thus,  if  the  absorbance 

a(v)  =  -ln[T(v)J  (eq.  2-10) 

is  measured  over  the  full  extent  of  the  line,  its  integral  will  be  propor¬ 
tional  to  the  density  of  the  absorbing  species.  Thus,  if  the  line  intensity  is 
known,  or  can  be  measured  from  a  known  density  (and  pathlength)  of  the 
gas,  then  the  density  can  be  measured  in  an  unknown  sample.  This  was  the 
method  chosen  for  measuring  the  silane  density. 

A  calibration  was  performed  by  operating  the  system  shown  in 
figure  2-1  in  a  static  fill  mode  (no  gas  flowing)  without  a  discharge.  For 
the  calibration,  a  gas  mixture  of  0.95%  silane  in  argon  was  used.  This  gas 
mixture  was  supplied  by  Matheson  and  was  guaranteed  to  be  within  10%  of 
the  most  significant  figure  in  the  mixture  specification.  Hence  the  mix 
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was  guaranteed  to  be  0.95±0.005%  silane  in  argon.  Next,  a  line  was  chosen 
to  be  the  calibration  line  for  the  measurements.  The  strong  line  on  the 
right  of  figure  2-5  was  chosen  as  a  calibration  line  because  it  was  strong 
enough  to  give  good  sensitivity  and  isolated  from  other  interfering  lines. 
This  line  has  been  identified  by  Johns  and  Kreiner35  as  the  Q(3),A2,1  line 
at  a  position  of  910.759+0.003  cm-1.  The  system  was  filled  to  a  known 
pressure  with  the  mixture,  then  nitrogen  added  (to  simulate  conditions  in 
the  experiment)  and  a  scan  of  the  line  taken.  The  diode  output  was  ad¬ 
justed  so  that  the  zero-gas  background  intensity  was  flat  across  the 
width  of  the  line  and  on  either  side  for  at  least  four  full  widths  of  the 
line.  After  much  trial  and  error,  it  was  found  just  as  accurate  to  fit  a 
straight  line  across  the  top  of  the  line  and  to  use  this  as  the  zero-gas 
reference  as  to  scan  the  region  twice  (with  and  without  gas)  and  ratio 
these  files.  The  measured  l(v)  was  then  ratioed  against  this  constructed 
reference  to  obtain  a  transmittance  spectrum.  An  absorbance  spectrum 
was  then  calculated,  and  two  procedures  were  then  used  to  obtain  the  line 
intensity.  First,  the  absorbance  spectrum  was  integrated  numerically  over 
the  full  extent  of  the  scan.  Since  the  path  length  and  silane  partial  pres¬ 
sure  were  known,  the  line  intensity  S  could  be  measured.  The  second 
procedure  was  to  perform  a  least-squares  fit  of  a  Voigt  function  to  the 
measured  absorbance.  The  method  used  for  the  least-squares  fit  was  the 
Marquardt  algorithm  which  is  described  in  Bevington36.  Voigt  functions 
were  generated  using  an  algorithm  described  by  Drayson37.  This  later 
method  gave  much  more  consistent  results.  Five  measurements  were  made 
with  the  result 

S[Q(3),A2,1]  =  9.74±0.44  cm^atm"1. 

The  two  procedures  described  above  were  then  used  along  with  this  value 
of  the  intensity  to  measure  silane  density  downstream  of  the  afterglow. 

Visible  and  U.V.  Spectrometer 

The  spectrometer  used  to  acquire  visible  and  u.v.  emission  spectra 
of  the  reaction  flame  was  a  0.5  meter  Spex  monochromator  with  independ¬ 
ently  adjustable  entrance  and  exit  slits.  The  reaction  zone  was  imaged  1:1 
on  the  entrance  slit  using  a  50  mm  diameter  CaF2  lens  of  150  mm  focal 
length.  A  1200  lines/mm  grating  blazed  at  500nm  was  used  throughout 
the  measurements.  With  this  blaze  the  system  could  respond  down  to 
about  220  nm.  The  photomultiplier  used  was  an  RCA  C31034  which  is 
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sensitive  over  the  range  of  about  200  to  890  nm.  This  photomultiplier  has 
a  GaAs  photocathode  and  should  be  cooled  for  optimum  performance. 

Cooling  was  supplied  by  a  Products  for  Research  thermoelectric  cooler 
which  in  turn  was  water  cooled.  This  system  allowed  the  photocathode  to 
be  cooled  to  -20°C  which  resulted  in  a  reduction  in  dark  current  to  less 
than  2  nA.  Two  different  systems  were  used  to  measure  the  PMT  output 
during  the  course  of  the  experiment.  Early  in  the  experiment,  the  PMT 
anode  was  simply  connected  to  ground  through  a  20  Mohm  resistor  which 
had  a  0.01  uf  capacitor  in  parallel.  The  voltage  across  the  resistor  was 
measured  with  a  Hewlett  Packard  model  3478A  digital  voltmeter  operated 
in  the  auto  range  mode.  The  output  of  the  voltmeter  was  sent  to  the  9836 
computer  over  an  IEEE  488  line.  The  second  system  for  measuring  the  PMT 
output  consisted  of  a  Photochemical  Research  Associated  model  1770  dual 
channel  photon  counting  system  with  an  interface  module  for  transfer  of 
data  to  the  488  bus.  In  this  configuration,  the  PMT  dark  counts  were  on  the 
order  of  5  cps.  Both  systems  worked  very  well  and  gave  good  dynamic 
range  and  signal-to-noise  in  the  spectra.  The  latter  system  only  became 
available  near  the  end  of  the  experiment. 


Relative,  Calibratifln  ol  the .  Spectrometer 

Early  in  the  experiment,  it  was  found  that  intense  emission  was 
observed  in  the  wavelength  range  220-500  nm.  It  was  decided  to  calibrate 
the  spectrometer  for  relative  response  over  this  wavelength  range. 
Unfortunately  no  standard  sources  were  available  in  our  laboratory  that 
possessed  sufficient  intensity  below  300  nm  to  allow  calibration.  Hence 
the  system  was  calibrated  between  300  and  500  nm.  To  perform  the 
calibration,  a  tungsten  ribbon  standard  lamp  (Eppley  Laboratories)  was 
placed  at  the  focus  of  the  lens  in  figure  2-1  in  place  of  the  discharge  tube. 
The  difficulty  with  this  calibration  lies  in  the  fact  that  at  300  nm  the 
scattered  light  in  the  instrument  is  about  a  factor  of  3  larger  than  the 
radiation  from  the  standard  lamp.  To  determine  the  scattered  light 
contribution  at  a  given  wavelength,  two  sharp  cutoff  colored  glass  filters 
from  Corning  were  used.  Filter  #3-75  transmitted  down  to  380  nm  and 
filter  #3-72  cutoff  at  430  nm.  For  wavelengths  up  to  400  nm,  the 
procedure  was  to  record  the  signal  from  the  PMT  without  any  filters,  then 
with  filter  #3-72  and  then  with  filter  #3-75.  With  a  particular 
wavelength  blocked,  it  was  assumed  that  the  signal  was  due  to  scattered 
light.  But  this  was  the  scattered  light  present  with  the  filter  in  place.  To 
arrive  at  a  value  for  the  scattered  light  present  with  no  filter  in  place,  it 


was  assumed  that  up  to  400  nm  the  ratio  of  scattered  light  with  a  filter 
to  that  without  a  filter  was  a  constant.  Measurements  were  then  made  of 
the  scattered  light  at  200  and  220  nm  (where  the  lamp  output  is 
negligible)  with  and  without  the  filters  in  place.  The  ratios  of  these 
signals  were  then  used  to  scale  up  the  values  measured  for  scattered  light 
in  the  wavelength  range  of  interest.  An  attempt  was  made  to  use  a  third 
filter,  #3-67,  which  cutoff  at  550  nm  to  determine  the  scattered  light 
contribution  in  the  wavelength  range  400-500  nm.  Unfortunately,  the  data 
indicated  that  the  filter  cutoff  too  much  of  the  radiation  from  the  lamp  to 
accurately  determine  the  scattered  light  contribution  without  the  filter  in 
place.  Hence  it  was  decided  to  treat  the  scattered  light  contribution  as  a 
constant  over  the  range  400-500  nm  and  set  the  value  at  that  measured 
for  400nm  Fortunately,  the  lamp  signal  in  this  wavelength  range  was 
greater  than  a  factor  of  25  larger  than  the  scattered  light,  hence  the  error 
from  treating  the  scattered  light  value  incorrectly  should  be  less  than  4%. 
A  calibration  sheet  supplied  with  the  standard  lamp  gave  only  seven 
values  of  the  lamp  output  between  300  and  500  nm.  Signals  were 
measured  at  each  of  these  wavelengths,  scattered  light  contributions 
removed,  and  the  result  ratioed  to  the  radiance  of  the  standard  lamp.  The 
resulting  seven  'sensitivity  values'  for  the  system  were  then  least- 
squares  fit  to  a  third  order  polynomial  to  allow  the  relative  response  of 
the  system  to  be  calculated  at  any  wavelength  between  300  and  500  nm. 

Fourier  Transform  Spectrometer/Lono  Path  Cell 

During  part  of  this  work,  the  1.2  m  absorption  cell  shown  in  figure  2- 
1  was  replaced  by  a  20  m  absorption  cell  of  the  White-type38.  The  base 
length  of  the  cell  is  0.75  m,  but  by  reflecting  the  radiation  entering  the 
cell  back  and  forth,  a  maximum  path  length  of  20.25  m  can  be  achieved. 

This  cell  is  a  commercially  made  unit  from  Wilks.  Absorption  spectra  of 
species  in  the  cell  were  obtained  using  a  broadband  infrared  emitter 
(glowbar)  and  an  FT-IR  spectrometer.  The  cell,  spectrometer,  and  as¬ 
sociated  optics  are  shown  in  figure  2-6.  The  detector  used  was  a  liquid 
nitrogen  cooled  HgCdTe.  The  FT-IR  system  is  a  commercial  system  built  by 
EOCOM  and  interfaced  to  a  Nicolet  1180  minicomputer.  The  system  is 
capable  of  achieving  a  resolution  of  0.06  cm-1  and  with  the  KBr  beamsplit¬ 
ter  installed,  could  cover  the  wavelength  range  from  2.5  to  14  microns 
(4000-400  cm-1).  The  reason  for  using  the  FT-IR  system  instead  of  the 
diode  laser  is  for  its  great  scanning  range.  The  diode  laser  can  typically 
only  be  scanned  =0.5  cm-1  without  changing  the  heat-sink  temperature. 
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Figure  2-6:  FT-IR/Long-Path  Absorption  Cell  System 


Thus  the  diode  laser  is  useful  when  the  species  of  interest  is  well  iden¬ 
tified,  but  for  survey  work  the  FT  system  is  far  superior.  The  20  m  cell 
was  needed  to  give  the  FT  system  a  comparable  sensitivity  to  the  diode 
laser  for  low  pressure  gases.  Since  the  typical  linewidth  of  an  isolated 
line  at  1  torr  is  on  the  order  of  10-3  cm-1,  even  at  its  highest  resolution 
the  FT  system  is  a  factor  of  20  larger  in  FWHM.  Thus  peak  absorbances 
measured  with  the  FT  system  for  an  isolated  line  are  typically  a  factor  of 
20  smaller  than  those  measured  with  the  diode  laser. 
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Chapter  III 

Experimental  Results 


As  stated  in  chapter  I,  this  work  began  with  the  question  of  whether 
silane  would  be  dissociated  in  active  nitrogen.  After  the  apparatus 
described  in  chapter  II  was  assembled,  the  first  observations  made  were 
aimed  at  answering  that  question.  With  nitrogen  flowing,  no  discharge,  and 
a  small  amount  of  the  silane/helium  mixture  flowing,  the  diode  laser  was 
scanned  across  the  wavelength  region  shown  in  figure  2-5  to  observe  the 
silane  absorption.  With  the  diode  laser  signal  displayed  in  real  time  on  an 
oscilloscope,  the  discharge  was  initiated  and  the  silane  absorption  almost 
disappeared.  Indeed,  silane  was  dissociated  in  the  afterglow.  If  that  was 
all  that  had  happened,  the  experiment  probably  would  have  stopped  with  a 
measurement  of  the  rate  of  dissociation  of  the  silane.  But  before  any 
quantitative  measurements  were  made,  a  number  of  general  observations 
indicated  that  this  was  a  very  interesting  system  for  study. 

General  Observations  of  the  Reaction 

The  first  surprising  observation  was  that  at  the  point  where  the 
silane  was  introduced,  a  bright  lavender  reaction  flame  was  produced 
which  could  be  seen  with  the  room  lights  on  if  the  pressure  was  of  the 
order  of  few  torr  or  higher.  In  addition,  downstream  of  the  reaction  flame 
the  straw-yellow  afterglow  was  no  longer  visible.  This  reaction  flame  had 
also  been  described  by  Dewhurst  and  Cooper29  but  at  the  time  of  these 
first  observations  the  author  was  unaware  of  that  reference.  After 
observing  the  reaction  flame  for  only  a  few  minutes,  a  white  powder  could 
be  seen  coating  the  walls  of  the  reaction  vessel  and  also  the  silane 
injection  port.  However,  more  striking  than  the  flame  itself  was  its 
behavior  with  pressure,  silane  flow,  and  discharge  current.  For  example, 
at  a  few  torr  the  flame  could  be  made  to  move  far  downstream  of  the 
silane  injection  port  by  reducing  the  discharge  current  or  the  silane  flow, 
yet  the  physical  extent  of  the  flame  remained  about  the  same.  If  tfv  flame 
was  moved  far  enough  downstream  that  it  moved  out  of  the  absorption  cell 


toward  the  pump,  no  silane  dissociation  was  observed.  Under  slow 
nitrogen  flow  conditions,  the  flame  shape  changed  to  a  disk  which  moved 
upstream  of  the  silane  injection  port.  Under  some  conditions  of  pressure, 
discharge  current,  and  silane  flow  the  reaction  flame  could  be  made  to 
pulsate  at  about  1  Hz  or  less,  and  these  pulsations  could  be  seen  in  the 
silane  dissociation.  While  this  qualitative  behavior  was  very  repeatable, 
quantitatively  many  of  the  observations  appeared  to  be  influenced  by  the 
powder  on  the  tube  walls.  In  general,  as  the  tube  walls  became  coated,  the 
reaction  seemed  to  occur  more  slowly  (reaction  flame  moved  downstream) 
for  the  same  discharge  conditions  and  mass  flow  settings.  Because  of  the 
effects  of  powder  buildup  on  the  tube  walls,  frequent  cleaning  of  the  flow 
tube  was  required.  From  these  general  observations,  the  conclusion  was 
drawn  that  some  kind  of  chain  reaction  must  be  involved  in  both  the 
production  of  the  chemiluminescence  and  the  silane  dissociation.  Silane 
dissociation  by  direct  reaction  with  a  single  afterglow  species  was  ruled 
out  from  the  complex  behavior  of  the  reaction  flame.  Though  not  quan¬ 
titative,  these  general  observations  would  later  become  a  strong  guiding 
force  behind  the  description  of  the  overall  reaction  process. 

Reaction  Flame  Emission  Spectra 

Figure  3-1  shows  a  spectrum  of  the  reaction  flame  taken  at  a 
background  pressure  of  2  torr.  Resolution  (FWHM)  is  about  0.5  nm  and  the 
spectrum  is  not  corrected  for  system  response.  The  strongest  features  in 
the  spectrum  are  atomic  silicon  lines  at  250,  290,  and  300  nm,  SiN  emis¬ 
sion  between  390  and  430  nm,  and  an  unidentified  broadband  continuum 
emitter  between  450  and  870  nm.  Close  inspection  actually  shows  that 
the  continuum  extends  down  to  about  300  nm.  Figure  3-2  shows  a  more 
detailed  spectrum  of  the  Si  emission  at  a  resolution  of  0.1  nm.  In  figure  3- 
2,  the  log  of  the  intensity  is  plotted  so  the  weaker  lines  near  250  nm  can 
be  more  easily  seen.  Table  3-1  lists  the  Si  lines  observed.  The 
wavelengths  and  transitions  are  from  Striganov  and  Sventitskii39. 
Jevons40  has  observed  the  same  Si  lines  in  the  reaction  of  SiCU  in  active 
nitrogen.  Jevons  also  observed  an  additional  group  of  lines  between  220 
and  222  nm,  but  in  that  wavelength  range,  the  sensitivity  in  the  work 
reported  here  is  almost  zero.  The  only  other  lines  reported  by  Jevons 
besides  those  given  in  table  3-1  are  lines  at  243.8,  244.3,  245.2,  and 
253.2  nm.  Jevons  reports  that  for  each  of  these  lines,  the  intensity  is 
"weaker  than  in  a  discharge"  of  SiCU  vapor.  Since  these  lines  are  listed  as 
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REACTION  FLAME  EMISSION  AT  2  TORR 


Figure  3-1:  Emission  Spectrum  of  Reaction  Flame  at  2  torr 
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Figure  3-2: 


very  weak  in  Striganov  and  Sventitskii,  they  may  have  been  too  weak  to 
observe  in  this  work. 


Table  3-1 :  Si  Lines  Seen  in  Reaction  Flame 


Wavelength  (nml 

Transition1 

AK 

Relative  Intensity^ 

243.52 

3p2  1D-3d  1D° 

2-2 

67 

250.69 

3p2  ^s  3P° 

1-2 

656 

251.43 

3p2  3P-4s  3P° 

0-1 

536 

251.61 

3p2  3P-4s  3P° 

2-2 

2100 

251.92 

3p2  ’D^s  3P° 

1-1 

466 

252.41 

3p2  3P-4s  3P° 

1-0 

700 

252.85 

3p2  3P-4s  3P° 

2-1 

874 

263.13 

3p2  1P° 

0-1 

200 

288.16 

3p2  ’D^s  1P° 

2-1 

12300 

298.76 

3p2  ^s  3P° 

2-1 

334 

300.67 

3p2  3P-3p3  5S° 

1  -2 

10825 

302.00 

3p23P-3p3  sS° 

2-2 

32325 

390.55 

3p2  'S-4s  1P° 

0-1 

2548 

transition  listed  as  lower  state-upper  state 
total  J  vaules  are  upper  J-lower  J 

3Measured  from  data  at  10  torr;  not  corrected  for  spectrometer  response 

As  in  the  work  here,  Jevons  was  unable  to  identify  any  Si  lines  above 
390.5  nm  in  the  afterglow.  Jevons  does,  however,  point  out  that  lines 
above  390.5  nm  are  seen  very  clearly  in  a  discharge  of  SiCU.  Figure  3-3 
shows  a  more  detailed  spectrum  of  the  SiN  emission  which  originates 
from  the  B-X  transition  of  SiN41.  Vibrational  band  assignments  are  given 
as  (v'-v")  and  are  from  reference  43.  Besides  the  unidentified  broadband 
emitter,  only  Si  and  SiN  emission  could  be  identified.  No  NHX,  l-fe,  H,  or  SiH 
emission  could  be  observed  despite  careful  attempts  to  identify  emission 
from  these  species.  Dewhurst  and  Cooper29  also  reported  that  no  ammonia, 
hydrazine,  or  other  liquid  nitrogen  condensable  products  were  formed. 
Comparing  the  spectra  shown  here  with  the  recent  work  of  references  30 
and  31,  who  also  looked  at  silane  in  active  nitrogen,  the  spectra  in  this 
work  agree  qualitatively  with  those  of  Horie,  et  al.30,  but  differ 
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significantly  from  those  of  Piper  and  Caledonia31.  Piper  and  Caledonia 
observed  rather  strong  NH  emission  and  weak  SiH  emission  while  Horie,  et 
al.  could  not  detect  either  of  these  species.  The  SiN  emission  observed  by 
Piper  and  Caledonia  was  very  weak  while  in  the  work  of  Horie,  et  al.  the 
SiN  emission  was  a  dominant  feature  as  in  this  work.  In  addition,  Piper 
and  Caledonia  do  not  show  any  broadband  emission  above  450  nm,  but  show 
continuum  emission  between  about  280  and  430  nm.  On  the  other  hand, 
Horie,  et  al.  observe  strong  continuum  emission  between  450  and  810  nm 
but  show  no  broadband  emission  below  450  nm  The  differences  in  the 
experiments  are  primarily  in  the  construction  of  the  flowing  afterglow 
systems.  Horie,  et  al.  used  a  system  very  similar  to  the  one  used  here 
which  used  relatively  slow  gas  flow  while  Piper  and  Caledonia  used  a 
much  higher  flow  speed  system.  Specifically,  the  spectra  presented  in 
reference  30  were  taken  about  60  ms  after  the  silane  was  injected,  while 
the  spectra  presented  in  reference  31  were  acquired  at  about  10  ms  into 
the  reaction.  The  pressures  in  the  two  experiments  were  3  torr  and  1.5-2 
torr  respectively.  Thus,  Piper  and  Caledonia  are  observing  the 
chemiluminescence  from  very  early  times  in  the  reaction,  while  in  the 
work  of  both  Horie,  et  al.,  and  this  work,  the  observed  chemiluminescence 
is  from  much  later  in  the  reaction.  Scans  of  the  reaction  flame  given  in 
reference  30  indeed  show  that  for  early  times  in  the  reaction  both  the  Si 
and  SiN  emission  are  very  weak.  The  inability  to  characterize  the  flow 
conditions  in  this  work  prevented  quantitative  measurements  of  the  time 
dependence  of  the  chemiluminescence,  but  spectra  taken  at  different 
points  in  the  reaction  flame  also  indicated  a  buildup  of  the  Si  and  SiN 
emission.  These  data  support  the  idea  that  the  differences  in  the  observed 
spectra  are  a  result  of  the  measurements  being  made  at  different  times  in 
the  reaction  history.  On  the  other  hand,  it  was  found  in  this  work  that  the 
broadband  emission  recorded  at  800  nm  preceded  both  the  Si  and  SiN 
emission.  This  was  also  observed  by  Horie,  et  al.,  yet  as  mentioned  earlier 
Piper  and  Caledonia  show  no  long  wavelength  continuum  in  their  data.  In  a 
phone  conversation  with  Piper,  he  stated  that  looking  carefully  at  his  data 
it  appeared  "the  chart  recorder  pen  did  not  go  to  zero"  above  450  nm.  This 
statement  together  with  the  fact  that  they  used  a  grating  blazed  for  the 
u.v.,  may  have  lead  to  the  continuum  being  overlooked.  It  appears 
reasonable,  therefore,  that  the  qualitative  differences  between  the 
spectra  observed  in  this  work  and  reference  30,  and  those  observed  in 
reference  31  are  probably  due  to  differences  in  the  time  of  observation  of 
the  reaction  history.  The  failure  of  reference  30  to  see  the  origin  of  the 
continuum  closer  to  300  nm  may  have  been  the  result  of  the  pressure  used 
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in  their  experiment.  This  pressure  effect  will  be  discussed  in  the  follow¬ 
ing  section. 

Pressure  Effects  on  the  Emission 

Figures  3-4  through  3-7,  along  with  figure  3-1  show  emission 
spectra  obtained  between  0.5  torr  and  8  torr  total  pressure.  Resolution  is 
approximately  0.5nm.  In  each  case  the  nitrogen  flow  was  adjusted  to 
maintain  constant  flow  speed.  The  silane/helium  mixture  was  adjusted  so 
that  the  afterglow  was  fully  quenched  downstream  of  the  reaction  zone. 
The  spectra  show  how  the  continuum  emission  decreases  relative  to  the  Si 
and  SiN  emission  as  the  pressure  is  increased.  Because  of  this,  positive 
identification  of  the  SiN  emission  at  0.5  torr  was  difficult.  The  spectra 
also  show  that  the  SiN  emission  increases  faster  than  the  Si  emission  at 
the  higher  pressures.  Beyond  these  observations,  it  is  difficult  to  draw 
too  many  quantitative  conclusions  from  the  data  due  to  the  behavior  of  the 
reaction  flame.  At  higher  pressures  (8-10  torr)  the  flame  is  very  intense 
and  very  compact,  much  like  a  bunsen  burner  flame.  As  the  pressure  is 
decreased,  the  intensity  of  the  flame  decreases  and  the  flame  becomes 
more  diffuse.  For  example,  at  0.5  torr  the  'flame'  completely  fills  the  36 
mm  section  of  the  flowing  afterglow  system  (see  figure  2-1),  appearing 
as  a  very  diffuse  glow.  On  the  other  hand,  the  emitting  volume  sampled  by 
the  spectrometer  was  essentially  fixed  by  the  collection  optics  and  slit 
area.  This  volume  was  roughly  2  cm  x  0.3  mm  (slit  height  x  width)  x  3.6 
cm  (diameter  of  flow  section).  For  good  quantitative  data,  some  kind  of 
integration  over  the  full  extent  of  the  flame  would  be  needed,  and  this 
was  not  available.  As  can  be  seen  by  comparing  figures  3-4  and  3-5  with 
figure  3-6,  it  is  much  clearer  in  the  spectra  taken  at  less  than  2  torr  that 
the  continuum  emission  extends  down  to  =300nm.  This  may  explain  the 
failure  of  Horie,  et  al.30  to  observe  any  continuum  below  450nm,  as  all  of 
their  data  were  taken  at  pressures  very  close  to  3  torr.  Besides  pressure, 
another  possible  explanation  for  the  increased  intensity  would  be  the 
increase  in  N-atom  density  with  pressure  in  the  experiment.  If  the 
chemiluminescence  is  driven  by  N-atoms,  then  this  increased  density 
might  explain  the  increased  intensity  of  the  emission.  To  determine  the 
relative  increase  in  N-atom  density,  the  intensity  of  the  N2  first  positive 
emission  in  the  afterglow  was  used.  Figure  3-8  is  a  spectrum  of  a  typical 
afterglow  from  this  experiment  taken  at  2  torr.  The  intense  emission  of 
the  first  positive  system  gives  the  afterglow  its  straw-yellow  color. 
While  there  are  numerous  theories  that  attempt  to  explain  the  origin  of 
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Figure  3-5:  Emission  Spectrum  of  Reaction  Flame  at  1  torr 


REACTION  FLAME  EMISSION  AT  4  TORR 


Figure  3-6:  Emission  Spectrum  of  Reaction  Flame  at  4  torr 


NITROGEN  AFTERGLOW  EMISSION  AT  a  TORR 


Figure  3-8:  Emission  Spectrum  of  Nitrogen  Afterglow  at  2  torr 


the  first  positive  emission,  experimentally  it  is  observed  that  in 
the  pressure  range  of  1-10  torr,  the  intensity  of  the  first 
positive  afterglow  emission  is  proportional  to  the  square  of  the 
N-atom  density.42  Thus  measurements  of  the  first  positive  intensity 
before  silane  was  introduced  could  be  related  to  the  N-atom  density  in  the 
afterglow.  Before  each  of  the  spectra  shown  in  figures  3-4  through  3-7 
(and  3-1)  were  taken,  a  spectrum  of  the  first  positive  emission  was 
recorded  without  silane  flowing.  The  first  positive  spectra  were  then 
numerically  integrated  from  565nm  to  610nm  and  this  integral  used  as  a 
measure  of  the  first  positive  intensity.  For  example,  between  1  and  8  torr 
the  SiN  emission  increased  by  a  factor  of  360  while  the  first  positive 
intensity  only  increased  by  a  factor  of  68.  This  indicates  that  the  N-atom 
density  only  increased  by  a  factor  of  8.2.  The  Si  emission  likewise 
increased  more  than  the  first  positive  intensity.  This  indicates  that  there 
is  an  effect  of  total  pressure  on  both  the  Si  and  SiN  emission  which 
exceeds  the  rate  of  increase  one  might  expect  from  increased  N-atom 
density  at  higher  pressures.  While  it  is  difficult  to  be  very  precise  in  any 
quantitative  analysis  of  the  emission  data  for  the  reasons  stated  above,  it 
does  appear  that  both  the  Si  and  SiN  emission  increase  faster  than 
linear  with  pressure,  even  when  scaled  to  the  increased  N-atom 
density. 

SiN  B-State  Vibrational  Populations 

Figure  3-9  shows  a  set  of  Morse  potentials  for  SiN  constructed  from 
the  molecular  constants  given  by  Foster43.  Foster  does  not,  however,  give 
values  for  the  dissociation  limits.  The  dissociations  limits  were  taken 
from  Bredohl,  et  al.41  who  extrapolated  the  ground  state  potential  to  a 
value  of  D0=5.77  eV.  The  ground  state  dissociation  energy  given  by 
Herzberg44  is  Do=4.5  eV.  Since  the  dissociation  energy  of  SiN  is  important 
for  determining  the  energetics  of  various  reactions  in  the  afterglow,  a 
more  detailed  discussion  of  the  dissociation  energy  is  given  in  appendix  A. 
As  stated  earlier,  the  emission  shown  in  figure  3-3  is  from  the  B-X 
transition.  Bands  originating  from  the  first  "<0  vibrational  levels  of  the  B- 
state  could  be  identified  in  spectra  taken  during  the  experiment.  Using 
available  Frank-Condon  factors45,  the  relative  populations  of  the  B-state 
were  calculated  from  the  measured  intensities  of  the  B-X  emission.  These 
relative  populations  are  shown  in  figure  3-10.  These  data  were  obtained 
from  intensity  measurements  taken  at  10  torr  total  pressure.  The 
procedure  used  to  obtain  the  data  was  to  first  correct  the  spectra  for  the 
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MORSE  POTENTIALS  FOR  SiN 
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Figure  3-9:  Morse  Potentials  for  SiN 


MEASURED  RELATIVE  SiN  B-STATE  VIBRATIONAL 
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Figure  3-10:  Measured  Relative  SiN  B-State  Vibrational  Populations 


relative  response  of  the  spectrometer  system  using  the  calibration 
procedure  described  in  chapter  2.  Next,  intensities  of  the  R-branch 
bandheads  were  recorded  for  each  transition  that  could  be  observed.  Each 
resulting  intensity  value  was  then  divided  by  the  corresponding  Frank- 
Condon  factor  to  arrive  at  a  set  of  relative  populations.  Finally,  relative 
populations  having  a  common  upper  level  were  grouped  and  a  weighted 
average  taken.  The  weights  were  set  proportional  to  the  intensity  of  the 
transition.  These  average  values  of  the  relative  populations  of  the  upper 
levels  were  then  scaled  to  the  maximum  value  of  the  population  of  v=5.  No 
error  bars  are  given  in  figure  3-10,  but  reference  45  gives  the  Frank- 
Condon  factors  as  accurate  to  between  5%  and  10%.  Probably  an  error  of 
1 0%-20%  is  reasonable  for  the  data  of  figure  3-10.  As  can  be  seen,  the 
data  clearly  shows  a  peak  at  v=5  of  the  B-state.  This  distribution  is 
unusual  in  that  one  would  normally  expect  the  v=0  state  to  be  the  most 
populated.  What  is  perhaps  even  more  unusual  is  that  an  almost  identical 
vibrational  population  distribution  has  been  seen  by  Broida46  in  the  SiN 
emission  from  the  reactions  of  SiCU  and  SiBr4  in  active  nitrogen.  Thus, 
the  vibrational  populations  in  the  B-state  of  SiN  follow  the 
same  distribution  whether  the  starting  'fuel'  is  SiH4,  SiCU >  or 
SiBr4. 

Silane  Dissociation  Measurements 

As  discussed  at  the  beginning  of  this  chapter,  silane  was  indeed 
dissociated  in  the  reaction  with  active  nitrogen.  Since  a  dc  discharge  was 
used  to  produce  the  afterglow  in  this  experiment,  it  was  first  decided  to 
measure  the  amount  of  silane  dissociation  as  a  function  of  the  discharge 
parameters  such  as  current,  voltage,  and  pressure.  Early  on  in  the  inves¬ 
tigation,  it  was  felt  that  the  afterglow  produced  from  the  dc  discharge 
could  be  modeled  with  available  data  from  the  literature;  therefore, 
measurements  of  the  discharge  parameters  could  be  later  correlated  to 
specific  characteristics  of  the  afterglow.  Unfortunately,  the  modelling 
proved  to  be  unreliable  at  predicting  detailed  behavior  of  the  afterglow 
so  that  it  was  very  difficult  to  draw  conclusions  from  the  measurements 
of  the  dissociation  as  a  function  of  discharge  parameters.  For  example, 
model  predictions  of  the  N-atom  density  were  as  high  as  10%  while 
measurements  by  Caledonia,  et  al.47  give  typical  densities  of  less  than 
1%.  It  was  then  decided  to  try  and  correlate  the  dissociation  measure¬ 
ments  with  some  other  parameters  of  the  experiment  which  were  more 
easily  related  to  some  characteristic  of  the  afterglow. 
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As  discussed  earlier,  the  N2  first  positive  intensity  in  the  afterglow 
is  proportional  to  the  square  of  the  N-atom  density.  It  was  therefore 
decided  to  measure  the  silane  dissociation  as  a  function  of  the  first 
positive  emission  intensity.  To  do  this,  the  apparatus  shown  in  figure  2-1 
was  modified  slightly.  To  monitor  the  first  positive  intensity,  a  small 
side-looking  photomultiplier  with  a  slit  and  narrow  band  filter  was  placed 
just  upstream  of  the  silane  injection  tube.  Using  this  system,  the 
afterglow  intensity  could  be  measured  before  the  silane  was  introduced 
into  the  reactor.  The  PMT  was  a  Hamamatsu  R106  and  the  filter  was  a 
Baird  Atomic  #10-30-6  which  had  a  center  wavelength  of  578  nm  and  a 
bandwidth  of  approximately  13  nm.  It  was  found  necessary  to  locate  the 
PMT  slightly  away  from  the  silane  injection  port  because  first  positive 
emission  from  the  discharge  was  scattered  from  the  injection  tube  and 
overwhelmed  the  afterglow  signal.  Even  with  the  offset  of  the  PMT,  a 
small  amount  cf  scattered  light  from  the  discharge  still  entered  the  PMT 
and  had  to  be  corrected  for  in  the  measurements.  The  correction  procedure 
involved  quickly  reducing  the  nitrogen  mass  flow  while  allowing  the 
pressure  control  to  maintain  constant  pressure  in  the  system.  When  gas 
flow  reached  zero,  the  afterglow  intensity  effectively  reached  zero  and 
the  PMT  reading  was  taken  as  the  scattered  light  from  the  discharge.  This 
value  was  then  subtracted  from  the  reading  taken  with  nitrogen  flowing. 
The  PMT  signal  was  read  out  through  a  current-to-voltage  converter  con¬ 
nected  to  a  digital  voltmeter.  To  measure  the  silane  dissociation,  the 
following  procedure  was  used.  After  the  discharge  had  stabilized  and  the 
afterglow  intensity  measured,  a  small  amount  of  silane/helium  mixture 
was  introduced  and  the  downstream  silane  density  measured  with  the  TDL 
system.  The  density  was  determined  using  the  procedures  described  in 
chapter  2,  particularly  the  method  of  least-squares  fit  to  the  Voigt 
profile.  Depending  on  the  strength  of  the  absorption,  the  silane/helium 
mass  flow  was  either  increased  or  decreased  so  as  to  give  a  set  of 
downstream  silane  densities  that  could  be  accurately  measured  but  as 
close  to  zero  as  possible.  This  procedure  followed  the  approach  of  a 
titration  measurement  to  determine  what  silane/helium  mass  flow  would 
yield  zero  silane  density  downstream.  The  data  were  then  plotted  as  in 
figure  3-11  and  least-squares  fit  to  a  straight  line.  The  intercept  was 
then  taken  as  the  maximum  helium/silane  mass  flow  for  complete  silane 
dissociation.  This  raw  value  of  mixture  mass  flow  was  then  multiplied  by 
0.009  to  convert  to  silane  mass  flow.  This  value  was  taken  as  the  'dis¬ 
sociated  silane'.  Figure  3-12  shows  a  log-log  plot  of  the  first  positive 
intensity  vs.  the  dissociated  silane.  The  straight  line  is  a  least-squares 
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HELIUM/SILANE  MASS  FLOW  VS.  SILANE  DENSITY 
DOWNSTREAM  OF  REACTION  ZONE 
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DISSOCIATED  SILANE  vs.  FIRST  POSITIVE 
EMISSION  INTENSITY 


ssociated  Silane  vs.  First  Positive  Emission  Intensity 


fit  which  yields  a  slope  of  0.45±0.05.  This  indicates  that  the  silane 
dissociation  is  proportional  to  the  first  positive  intensity  to  the  0.45 
power.  Since  the  first  positive  intensity  is  proportional  to  the  square  of 
the  N-atom  density,  these  data  indicate  that  to  within  the  uncertainty  of 
the  measurements,  and  for  small  amounts  of  silane,  the  amount  of 
silane  dissociated  is  proportional  to  the  N-atom  density  in  the 
afterglow. 

FT-IR  Absorption  Measurements 

As  discussed  in  chapter  2,  during  this  part  of  the  experiment  the  1.2 
m  absorption  cell  in  figure  2-1  was  replaced  with  the  FT-IR/Long  Path 
Cell  system  shown  in  figure  2-6.  The  purpose  of  these  measurements  was 
to  survey  the  species  leaving  the  reaction  zone  with  absorption  spectros¬ 
copy  in  hopes  of  identifying  species  which  might  not  emit  in  the  visible. 
Initially,  measurements  were  made  at  4  cm-1  resolution  with  the  path 
length  at  at  20.25m.  Figure  3-13  is  a  transmittance  spectrum  of  the 
reaction  products  taken  at  a  total  pressure  of  10  torr  and  with  the  silane 
flow  adjusted  to  ensure  more  silane  than  could  be  dissociated  by  the 
afterglow.  The  sharp  features  at  *910  cm-1  and  *2180  cm-1  are  due  to 
unreacted  silane.  Likewise,  the  rapidly  varying  structures  on  each  side  of 
the  sharp  feature  at  2180  cm-1  are  also  due  to  silane.  The  structure 
centered  at  about  1500  cm-1  is  due  to  residual  water  vapor  in  the  lab  air 
which  did  not  "ratio  out"  when  the  transmittance  spectrum  was  cal¬ 
culated.  The  broad  features  at  900  and  2200  cm-1  which  overlap  the  silane 
are  due  to  new  species.  Figure  3-14  shows  that  if  the  silane/helium  flow 
is  reduced  so  that  all  of  the  silane  is  reacted,  these  broad  features 
remain.  Absorption  spectra  taken  at  0.06  cm-1  resolution,  though  much 
noisier,  show  the  same  broad  features.  No  structure  beyond  what  is  seen 
at  4  cm-1  resolution  could  be  identified.  Experiments  run  at  total  pres¬ 
sures  of  6  torr  and  5  torr  also  showed  these  features;  however,  at  5  torr 
the  absorption  was  approaching  the  detection  limit  of  the  system.  The 
lack  of  structure  in  the  observed  absorption  indicates  the  species  is 
probably  a  complex  molecule  with  many  degrees  of  freedom.  Positive 
identification  of  the  absorbing  species  has  not  been  made  at  this  time,  but 
some  suspected  species  have  been  ruled  out.  The  features  do  not  agree 
with  published  spectra  of  ammonia  or  hydrazine48.  This  is  consistent  with 
the  results  of  Dewhurst  and  Cooper29.  Comparison  with  absorption  spectra 
of  trisilylamine,  N(SiH3)3,  has  ruled  out  this  species  even  though  it 
possesses  absorption  bands  centered  in  exactly  the  same  locations  as  the 
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ABSORPTION  SPECTRUM  OF  REACTION  PRODUCTS  IN 
LONG-PATH  CELL  AT  10  TORR,  4  cm1  RESOLUTION 
(RESIDUAL  SIH4  SEEN  AT  910  AND  2180  cm*1) 
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Figure  3-13:  Long-Path  Absorption  of  Downstream  Reaction  Products 
(residual  silane  present) 
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ABSORPTION  SPECTRUM  OF  REACTION  PRODUCTS  IN 
LONG-PATH  CELL  AT  10  TORR,  4  cm*1  RESOLUTION 
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Figure  3-14:  Long-Path  Absorption  of  Downstream  Reaction  Products 
(reduced  silane  flow) 


observed  features.49  It  does  seem  reasonable  that  the  unidentified  species 
contains  Si,  N,  and  H  since  the  band  at  2180  cm-1  is  characteristic  of  an  Si- 
H  bond  stretching  vibration  and  the  Si-N  bond  stretching  is  at  about  1100 
cm-1.  It  is  very  tempting  to  assign  the  absorption  feature  to  SiaN4  which 
possesses  a  strong  band50  between  600  and  1250  cm-1,  however  the  band 
seen  here  at  2180  cm-1  would  indicate  that  some  hydrogen  is  also  present 
in  the  unidentified  species  or  perhaps  the  absorption  is  the  result  of 
multiple  species.  Therefore,  at  this  point  the  species  cannot  be  positively 
identified;  however,  it  appears  reasonable  to  assume  the  species  is  a 
higher  nitride  of  silicon  with  perhaps  some  hydrogenation.  This  indicates 
that  some  higher  nitrides  of  silicon  are  formed  in  the  reaction 
and  remain  in  the  gas  phase  at  pressures  as  low  as  5  torr. 

Disilane  Production  bv  Silane  Reaction  with  Nitrogen  Metastables 

While  the  FT-IR/Long  Path  Cell  system  was  still  in  place,  one  more 
set  of  measurements  was  made  with  a  modified  flow  reactor.  The 
modified  reactor  is  shown  in  figure  3-15.  In  this  system,  pure  argon  was 
flowed  through  the  discharge  to  produce  argon  metastables.  The  Ar  flow 
was  then  intersected  with  a  flow  of  pure  nitrogen.  The  gas  was  then 
reacted  with  the  silane/helium  mixture  further  downstream.  Stedman  and 
Setser51  have  shown  that  argon  metastables  will  readily  transfer  energy 
to  the  C3n  state  of  nitrogen  which  will  eventually  radiate  down  to  the 
metastable  A3X  state  of  N2.  This  reaction  is  fast  enough  that  conditions 
could  be  adjusted  so  that  all  of  the  argon  metastables  could  be  quenched 
by  the  time  the  flow  reached  the  helium/silane  port.  Thus  at  the  point  of 
silane  injection,  the  gas  flow  only  contained  groundstate  Ar  and  N2,  and 
some  metastable  N2(A3X).  All  measurements  made  with  this  configuration 
were  at  a  total  pressure  of  10  torr.  An  attempt  was  made  to  measure  any 
chemiluminescence  from  the  reaction,  but  to  the  sensitivity  of  the  ex¬ 
periment,  no  chemiluminescence  could  be  detected.  Piper  and  Caledonia31 
also  attempted  to  measure  chemiluminescence  from  the  reaction  of  N2 
metastables  with  silane,  but  likewise  were  unsuccessful.  From  the 
description  of  their  experiment,  it  is  felt  that  the  overall  sensitivity  of 
Piper  and  Caledonia's  experiment  was  better  than  this  work  below  about 
400  nm  but  above  that  wavelength,  the  experiment  described  here  was 
more  sensitive.  In  an  attempt  to  observe  whether  any  reaction  occurred, 
absorption  spectra  were  taken  of  the  gases  leading  the  reaction  zone  and 
entering  the  20.25  m  cell.  Initial  inspection  of  the  absorption  spectra 
showed  only  SiFU  present  in  the  cell.  Using  a  background  subtraction 
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ARGON  METASTABLE  TRANSFER  EXPERIMENT 
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Figure  3-15:  Schematic  of  Argon  Metastable  Transfer  Experiment 


scheme,  the  feature  shown  in  figure  3-16  (recorded  at  a  resolution  of  1 
cm-1)  was  uncovered.  In  this  subtraction  scheme,  an  absorbance  spectrum 
taken  without  a  discharge  was  subtracted  from  a  spectrum  with  a  dis¬ 
charge.  Before  the  subtraction  was  performed  the  'no  discharge'  spectrum 
was  multiplied  by  a  constant  less  than  1.  This  constant  could  be  varied 
about  1  while  the  difference  spectrum  was  observed  interactively.  When 
the  value  of  the  constant  was  such  that  a  flat  background  in  the  difference 
spectrum  was  obtained,  the  difference  spectrum  was  stored.  This  'absor¬ 
bance  difference  spectrum'  was  then  converted  back  to  the  transmittance 
spectrum  shown  in  figure  3-16.  While  the  feature  is  small  (note  the 
vertical  scale  in  figure  3-16),  the  absorbing  species  could  be  positively 
identified  as  disilane,  Si2H6-  Previous  work  in  this  lab  with  disilane  made 
possible  an  estimate  of  the  partial  pressure  of  the  species  from  the 
measured  absorbance.  Using  a  value  of  the  Q-branch  peak  absorption 
coefficient,  measured  in  this  lab  at  a  resolution  of  1  cm-1  by  Chaney,  et 
al.52,  of  ko=0.105±0.008  torr-1  cm-1  gives  a  value  of 

P Si2H6=  9.5±2.5  x  10-5  torr 

in  figure  3-16.  The  positive  identification  of  disilane  clearly  shows  that 
some  silane  is  dissociated  by  N2(A3£)  and  the  dissociated  products  react 
to  form  disilane.  The  most  likely  scheme  for  the  disilane  production  is 

SiH4  +  N2(A3I)  -  SiH2  +  H2  +  N2(X1Z) 
followed  by 

S1H4  +  SiH2  — ♦  Si2H6  — »  Si2H6- 
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Si2H6  ABSORPTION  SEEN  IN  LONG-PATH  CELL  AFTER 
SUBTRACTION  OF  RESIDUAL  SIH4  ABSORPTION 

O 

O 


Figure  3-16:  Disilane  Absorption  After  Substraction  of  Residual 
Silane  Absorption 


58 


Chapter  IV 

Development  of  the  Reaction  Model 


To  fully  understand  the  behavior  of  silane  in  active  nitrogen,  it  is 
necessary  to  know  what  chemical  reactions  take  place  and  the  rates  of 
those  reactions.  In  general,  this  is  a  very  difficult  problem.  In  light  of  the 
lack  of  detailed  understanding  of  the  nitrogen  afterglow  alone,  it  would  be 
foolhardy  to  suppose  that  such  a  detailed  picture  would  emerge  from  this 
investigation.  On  the  other  hand,  it  might  be  possible  to  propose  some 
likely  reactions  that  could  explain  the  observed  results.  Thus  an  effort  to 
"model"  the  reactions  was  begun.  Almost  by  definition,  this  part  of  the 
work  is  very  speculative.  Any  proposed  reaction  scheme  might  not  be  the 
only  scheme  that  would  describe  the  results,  however,  if  a  proposed 
scheme  is  plausible,  does  not  violate  any  laws  of  physics,  and  describes 
the  observed  results,  it  is  viable  until  some  new  result  or  other  data  rules 
it  out.  Therefore,  the  goal  of  this  part  of  the  investigation  was  to  arrive 
at  a  plausible  set  of  chemical  rerotions  which  would  explain  the  bulk  of 
the  observations.  Guidance  for  this  part  of  the  work  comes  from  a  number 
of  sources.  Wright  and  Winkler4  give  a  good  compilation  of  work  involving 
the  introduction  of  various  species  into  active  nitrogen.  Many  of  the 
references  therein  give  suggested  models  of  reaction  schemes  for  dif¬ 
ferent  species.  More  recently,  Horie  et  al.30  have  suggested  a  reaction 
scheme  based  on  their  measurements  of  silane  in  active  nitrogen  and  their 
work  has  guided  the  work  here.  Since  this  is  the  only  known  data  on  the 
reaction  mechanism  of  silane  in  active  nitrogen,  comparisons  will  be  fre¬ 
quently  made  between  their  work  and  the  work  reported  here.  When  the 
work  of  Horie  et  al.  is  cited,  it  should  be  understood  this  refers  to 
reference  30. 

General  Rate  Equation  Approach 

The  "model"  of  the  reaction  will  be  a  set  of  coupled  rate  equations 
which  govern  the  time  dependence  of  the  various  species  and  processes 
(such  as  chemiluminescence)  involved  in  the  reaction.  The  rate  equations 
are  essentially  a  set  of  coupled  differential  equations  which  can  be  solved 
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numerically  and  the  solutions  compared  with  observations.  For  example, 
suppose  the  following  reaction  occurs  in  a  system: 

A  +  B  -♦  C  (eq.  4-1) 

The  rate  equations  governing  the  species  would  be 

d[A]/dt  =  -K.[A].[BJ 

d[B]/dt  =  -K.[AJ.[B]  (eq.  4-2) 

d[C]/dt  =  +K.[A]-[B] 

where  [A],  [B],  and  [C]  are  number  densities  of  the  species  and  K  is  the  rate 
constant  for  that  reaction.  For  a  given  set  of  initial  conditions  one  can 
solve  this  set  of  coupled  equations  to  get  the  time  dependence  of  the 
species  involved.  From  this  simple  example,  one  can  see  that  the  number 
of  coupled  equations  is  equal  to  the  number  of  species,  and  the  number  of 
rate  constants  is  equal  to  the  number  of  reactions.  In  the  case  of  model¬ 
ling  a  completely  new  system  one  is  faced  with  putting  together  both  a 
set  of  proposed  reactions  and  a  set  of  proposed  rates  (if  measured  values 
do  not  exist).  What  makes  it  possible  to  bring  some  reality  to  such  a  model 
is  that  the  rates  one  can  choose  are  not  arbitrary.  The  rate  constants  basi¬ 
cally  are  a  measure  of  the  collision  rate  of  the  species  involved  in  a 
reaction.  This  means  the  rates  can  not  exceed  the  "gas  kinetic  rate",  i.e. 
the  rate  constant  calculated  by  assuming  the  particles  move  at  thermal 
velocities  and  collide  with  a  "hard  sphere"  collision  cross  section53. 
Specifically  for  two  particles  of  molecular  weight  m-i  and  m2  at  a 
temperature  T,  the  "gas  kinetic  rate"  is  given  by53 

r  1 2 

K(gas  kinetic)  »o12-[8jckT/m1  J  (eq.  4.3) 


where 


ai2-  0.5-fa1+a2)  (eq  4.4) 

with  ai  and  °2  the  hard  sphere  diameters,  mi  2  is  the  reduced  mass  and  k  is 
Boltzmann's  constant.  For  example,  if  one  assumes  for  each  colliding 
partner  a  hard  sphere  diameter  of  2.5  A,  a  molecular  weight  of  30amu,  and 
a  gas  temperature  of  300°K,  the  gas  kinetic  rate  is 
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K  =  1.3x10-1°  cm3/s. 


Therefore,  if  a  reaction  scheme  is  proposed  with  the  binary  collision  rate 
set  by  equation  4-3,  and  the  predicted  time  dependance  of  an  observed 
phenomenon  is  slower  than  measured,  the  scheme  is  in  trouble.  Another 
characteristic  of  this  system  that  makes  the  modeling  effort  tractable  is 
that  for  the  most  part  only  binary  collisions  need  be  considered  due  to  the 
low  pressures  involved  in  the  reactor. 


Characteristics  of  the  Reaction  Scheme 


Before  any  reaction  schemes  to  explain  the  behavior  of  silane  in 
active  nitrogen  could  be  developed,  some  knowledge  of  the  characteristics 
of  active  nitrogen  was  needed.  Caledonia,  et  al.47  have  made  measurements 
of  various  species  present  in  a  'typical'  afterglow  and  table  4-1  gives 
estimates  of  densities  of  various  species  one  might  expect  to  be  present 
at  a  pressure  of  3  torr. 

Table  4-1:  Catalogue  of  Abundant  Species  in  Active  Nitrogen 

(at  3  torr) 


Afterglow  Species  Approximate  Density  (cm=3.1 


N(4S) 

N(2D) 

N2(A3S) 

N(2P) 

N2(a’1I) 


1014 

1011 

1010 

1010 

<1010 


Table  4-1  must  be  taken  only  as  a  guide  because  conditions  such  as  gas 
purity,  cleanliness  of  the  flowing  afterglow  system,  type  of  discharge 
used  for  excitation,  and  a  host  of  other  factors  can  dramatically  affect 
the  number  densities  given  in  the  table.  Table  4-1  is  probably  a  reasonable 
guide  to  the  various  rankings  of  the  species  with  respect  to  number 


density.  One  species  no!  listed  in  table  4-1  is  vibrationally  excited, 
ground  electronic  state  N2.  The  number  densities  of  the  various  vibrational 
levels  depend  strongly  on  the  level  and  with  the  ground  state  potential 
capable  of  supporting  approximately  70  vibrational  levels,  there  are  too 
many  for  a  brief  table.  More  discussion  concerning  vibrationally  excited  N2 
will  be  given  later.  Thus,  table  4-1  indicates  that  in  most  nitrogen 
afterglows,  the  dominant  'active'  species  is  atomic  nitrogen.  Using  these 
general  characteristics  of  the  afterglow  and  the  results  presented  in 
chapter  III,  a  number  of  conclusions  about  the  reactions  occurring  in  the 
afterglow  were  formulated. 

The  following  is  a  list  of  the  results  and  the  conclusions  drawn: 

1)  From  the  general  observations  of  the  reaction  flame:  The  overall 
reaction  scheme  consisted  of  some  kind  of  chain  reaction  involving  more 
than  one  species  in  the  reaction  with  silane.  The  chemiluminescence  can 
take  place  in  time  short  compared  to  the  'induction  period',  i.e.  the  time 
before  the  chemiluminescence  begins. 

2)  From  the  observations  of  Si  and  SiN  emission  and  comparisons 
with  reports  using  other  silicon  compounds  as  'fuel':  The  Si  emission  and 
SiN  emission  resulted  from  reactions  only  involving  species  present  in  the 
afterglow,  atomic  silicon,  and  species  formed  from  atomic  silicon  and 
some  form  of  nitrogen. 

3)  From  the  general  observations,  the  dissociation  measurements 
and  the  N-atom  density  in  the  afterglow:  The  silane  dissociation  is  driven 
by  N-atoms;  however,  silane  does  not  readily  react  directly  with  N- 
atoms.* 

4)  From  the  FT-IR  absorption  measurements  and  the  results  of 
Dewhurst  and  Cooper29:  The  densities  of  any  NHX  species  must  be  very  low. 


The  latter  part  of  conclusion  3),  i.e.  silane  does  not  react  directly 
with  N-atoms  has  recently  been  confirmed  by  Piper  and  Caledonia  (ref.  31) 
who  measured  the  reaction  rate  of  N-atoms  with  silane  and  found  the  rate 
was  smaller  than  they  could  measure  with  their  apparatus: 

R<0. 006x1 0_1  ^rr^/s. 
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5)  From  the  results  of  emission  measurements  at  various  pressures: 
The  intensity  of  the  Si  and  SiN  chemiluminescence  increases  faster  than 
linear  with  pressure. 

These  general  conclusions  about  the  reaction  lead  to  a  very  general 
scheme  for  the  reaction  mechanism.  This  general  scheme  is  given  below: 

1)  Some  species,  [A],  present  in  the  afterglow  slowly  dissociates 

SiH4. 

2)  The  resulting  silicon-bearing  radicals  react  rapidly  with  N- 
atoms. 

3)  A  series  of  reactions  with  N-atoms  produces  the  observed 
chemiluminescence  and  produces  a  species,  [B],  which  can  dissociate  SiH4. 

4)  The  species  produced  in  step  3)  feeds  back  to  step  1)  to  dis¬ 
sociate  more  SiH4. 

In  the  above  scheme  the  two  unknown  species,  [A]  and  [B],  can  in  fact  be 
the  same:  however,  this  is  not  required.  The  problem  now  becomes  one  of 
filling  in  the  details  of  the  above  and  checking  if  indeed  the  results 
predicted  agree  with  observations. 

Initial  Silane  Dissociation  in  the  Afterglow 


To  determine  which  species  in  the  afterglow  might  be  responsible 
for  the  initial  dissociation  of  SiH4,  some  energetics  of  silane  dissociation 
are  needed.  The  energies  needed  to  dissociate  silane  thermodynamically 
into  various  radicals  are  shown  in  table  4-2.  The  dissociation  energies  are 
given  in  electron  volts.  The  dissociation  energies  are  calculated  from  heat 
of  formation  data  of  Ho,  et  al.54  Since  one  is  interested  in  dissociating 
silane  by  collision  with  some  energetic  species  in  the  afterglow,  table  4- 
2  is  somewhat  misleading  in  that  the  energies  apply  to  conditions  of 
thermodynamic  equilibrium.  In  a  collision  of  silane  with  an  energetic 
particle,  equilibrium  may  not  be  established  among  the  degrees  of  freedom 
of  the  particles.  As  discussed  in  chapter  I,  the  threshold  energy  for  dis¬ 
sociation  of  silane  by  electron  impact  is  believed  to  be  =7.8  eV,  even 
though  there  are  lower  energy  channels  as  shown  in  table  4-2. 


SiH4  — »SiH2  +  H2 


2.60 


SiH3  +  H 

3.98 

Si  +  H2  +  H2 

4.32 

SiH  +  H2  +  H 

5.89 

SiH2  +  H  +  H 

7.13 

Si  +  H2  +  H  +  H 

8.84 

SiH  +  H  +  H  +  H 

10.41 

Si  +  H  +  H  +  H  +  H 

13.37 

Data  to  support  the  hypothesis  that  dissociation  can  take  place  at  the 
lower  energies  shown  in  table  4-2  comes  from  two  sources.  In  chapter  III 
it  was  shown  that  Si2H6  was  formed  in  the  reaction  of  SiH4  with  ISfefA3!) 
metastables.  The  Si2H6  is  most  probably  produced  from  the  reaction  of 
SiH4  with  either  SiH2  or  S1H3  produced  by  the  dissociation  of  SiH4  by 
collision  with  N2(A3£).  Figure  4-1  is  an  energy  level  diagram  of  N2  based 
on  the  data  from  Huber  and  Herzberg55,  and  as  can  be  seen,  the  N2(A3£) 
state  lies  approximately  6.1  eV  above  the  ground  state  of  N2.  The  other 
indicator  of  lower  energy  dissociation  comes  from  experiments  with 
mercury  sensitization  of  silane24.  Silicon  deposition  by  mercury  sensitiza¬ 
tion  of  silane56  involves  deposition  by  irradiating  a  mixture  of  silane  and 
mercury  with  254nm  resonance  radiation  from  a  mercury  lamp.  It  is 
believed  that  the  deposition  is  driven  by  dissociation  of  SiH4  by  collisions 
with  the  Hg(6p3Pi°)  excited  state  which  lies  at  4.9  eV.  From  these  results, 
it  appears  that  dissociation  of  silane  by  heavy  particle  collisions  can 
occur  at  lower  energies  than  observed  for  electrons.  It  will  be  assumed 
then  that  the  energies  shown  in  table  4-2  are  valid  threshold  energies  for 
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plif ied) 


en  (simplified) 


dissociation  of  silane  by  heavy  particle  collisions.  From  this  assumption 
and  figure  4-1,  it  is  clear  that  a  number  of  energetic  species  of  N2  are 
available  to  dissociate  S1H4.  Both  the  metastable  A3I  and  the  a'1!  possess 
sufficient  energy  as  does  the  vibrationally  excited  X1!  state.  An  N^X1!) 
molecule  at  v=10  or  above  possesses  sufficient  energy  to  dissociate  SiH4 
by  the  low  energy  process  shown  in  table  4-2.  In  a  collision  of  this  type, 
the  N2  molecule  gives  up  a  large  number  of  vibrational  quanta  to  an 
electronic  state  of  SiH4  which  then  dissociates.  This  "v  to  E  transfer"  has 
been  observed  as  both  an  intramolecular  (within  the  same  molecule)  and 
as  an  intermolecular  (from  one  species  to  another)  process. 

Intramolecular  transfer  has  been  observed  by  Dunnwald,  et  al.57  in  nitric 
oxide,  and  by  DeLeon  and  Rich58  and  Farrenq  et  al.59  in  carbon  monoxide. 
Intermolecular  transfer  has  been  observed  by  Rich  and  Bergman60  in 
collisions  of  vibrationally  excited  CO  (and  N2)  with  C2  and  CN,  and  by 
Star61  in  collisions  of  vibrationally  excited  N2  with  sodium.  Thus,  there  is 
significant  evidence  to  suggest  that  vibrationally  excited  nitrogen  may  be 
capable  of  dissociating  silane  by  the  process  of  v  to  E  transfer.  Silane 
might  also  be  dissociated  by  v-v  pumping  in  collisions  with  vibrationally 
excited  N2.  In  this  case,  a  collision  is  expected  to  exchange  only  one 
vibrational  quanta  between  N2  and  SiH4.  By  repeated  collision  with  N2(v), 
SiH4  could  be  vibrationally  "pumped"  until  the  molecule  reached  the 
dissociation  limit.  As  mentioned  in  chapter  I,  SiFU  has  a  non-degenerate 
stretching  mode  at  2187  cm-1,  one  doubly  degenerate  bending  mode  v2  at 
978  cm-1,  and  two  triply  degenerate  bending  and  stretching  modes  V3  at 
2183  cm-1  and  V4  at  910  cm-1.  Exchange  of  vibrational  quanta  between 
N2(v)  and  SiH4  would  be  expected  to  occur  at  the  fastest  rate  between 
similarly  spaced  vibrational  levels.  Thus  one  might  expect  near-resonant 
transfer  between  either  vi  orv3  of  SiH4  and  v=6  of  N2(X1X).  The  lower 
energy  modes  of  SiFU  do  not  come  into  resonance  until  v=50  of  N2,  which 
is  too  high  to  be  considered  in  this  initial  dissociation  step.  To  model  this 
v-v  pumping  between  N2  and  SiH4  would  require  knowing  a  large  number  of 
rates  for  both  the  exchange  of  quanta  between  the  molecules  along  with 
rates  for  quenching  of  the  vibrational  energy  of  SiFU  by  various  gas  phase 
s  in  the  afterglow.  Due  to  the  complexity  of  this  calculation  and  the  lack 
of  rate  data,  it  was  simply  assumed  that  N2(X1I,  v>10)  would  have  suf¬ 
ficient  energy  to  dissociate  SiF^.  The  exact  mechanism  was  assumed  to  be 
v-to-E  transfer,  but  it  must  be  emphasized  that  v-v  pumping  may  also  be 
playing  a  role.  In  addition  to  the  energetic  species  of  N2,  the  metastable 
N(2P)  state  of  atomic  nitrogen  resides  at  3.57  eV  above  the  N(4S)  level, 
thus  it  may  also  be  capable  of  dissociating  SiFL}.  With  the  exception  of 
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vibrationally  excited  nitrogen,  the  estimated  densities  from  table  4-1  of 
the  remaining  species  are  on  the  same  order.  Calculations  which  will  be 
discussed  in  the  next  chapter  indicate  that  the  density  of  ^(X1!,  v=10)  is 
also  on  the  order  of  1010,  so  at  least  in  this  initial  model,  the  assumption 
is  made  that  any  of  the  four  species  just  discussed  may  contribute  to  the 
initial  dissociation  of  silane. 

Reactions  of  N-atoms  with  Silane  Radicals 

The  next  step  in  the  proposed  reaction  scheme  involves  direct  and 
rapid  reaction  of  N-atoms  with  radicals  produced  from  the  initial  dis¬ 
sociation  of  SiH4.  The  logic  behind  this  proposed  step  is  twofold.  First,  the 
observed  results  indicate  that  by  approximately  60ms  the  reaction  is 
complete  (at  3  torr).  Looking  at  table  4-1  and  assuming  a  gas  kinetic  rate, 
it  is  clear  that  N-atoms  must  be  involved  early  in  the  reaction  or  the 
reaction  will  evolve  too  slowly.  Second,  it  is  well  known  that  the  radicals 
produced  from  SiH4  and  N-atoms  are  all  (usually)  very  reactive.  Hence, 
what  is  contrary  to  logic  is  that  N-atoms  react  so  slowly  with  SiH4.  Table 
4-2  shows  that  the  three  silicon  bearing  radicals  and  atomic  silicon  may 
be  produced  in  the  dissociation  of  silane.  To  simplify  the  situation  some¬ 
what,  the  assumption  is  made  that  the  most  likely  radicals  produced  are 
those  originating  from  reactions  that  occur  at  the  lowest  energies  and 
require  the  fewest  hydrogen-silicon  bonds  to  be  broken.  For  example,  it  is 
difficult  to  visualize  a  collision  between  an  energetic  particle  and  SiH4  in 
which  all  four  hydrogen-silicon  bonds  are  broken  to  give  Si,  and  then  two 
hydrogen-hydrogen  bonds  reassembled  to  give  two  H2  molecules.  Using 
these  guidelines,  the  two  most  likely  radicals  to  be  produced  are  SiH3  (by 
line  2  of  table  4-2)  and  SiH2  (by  line  1  of  table  4-2).  The  possible  reac¬ 
tions  to  consider  in  this  step  will  therefore  involve  N-atoms  on  either 
SiH2  or  SiH3.  The  enthalpy  changes  for  each  suggested  reaction  will  be 
given  beside  each  equation  (Kcal/mole).  The  enthalpy  changes  are  cal¬ 
culated  from  the  heat  of  formation  data  given  in  appendix  B.  If  the  enthal¬ 
py  change  is  less  than  zero,  the  reaction  is  exothermic,  and  if  the  change 
is  greater  than  zero,  the  reaction  is  endothermic.  Since  this  reaction 
scheme  is  to  be  self  propagating,  reactions  which  are  endothermic  will  be 
ruled  out  as  possible  steps  in  the  chain.  For  S1H2,  some  possible  reactions 
are: 


SiH2  +  N  -SiNH  +  H  (-80) 


(eq.  4-5a) 
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HSiN  +  H 

-  SiH  +  NH 

-  SiN  +  H2 
-*  Si  +  NH2 

From  energy  considerations  reaction  4-5c  can  be  eliminated  from  the 
redaction  scheme.  Also  from  an  energetics  standpoint,  it  is  clear  that 
reactions  4-5a  and  4-5d  are  most  favored.  The  difficulty  with  reaction  4- 
5a  is  that  it  produces  H-atoms.  Since  no  H-atom  emission  was  seen  in  this 
experiment  one  is  tempted  to  conclude  that  the  H-atom  density  must  be 
low;  however,  no  H2  emission  is  seen  either,  and  Dewhurst  and  Cooper29 
found  approximately  stoichiometric  formation  of  H2  from  SiH4.  Still 
better  data  come  from  the  work  of  Horie,  et  al.  who  measured  the  time 
dependence  of  the  SiH4  dissociation  and  the  time  dependence  of  H2  density 
during  a  reaction.  Their  measurements  show  that  throughout  the 
time  history  of  the  reaction,  the  H2  formed  is  equal  to  twice  the 
SiH4  reacted.  So  that  even  early  on  in  the  reaction  scheme,  the  density 
of  atomic  hydrogen  must  be  low.  Therefore,  reactions  that  produce  H- 
atoms  are  considered  secondary  in  importance  to  reactions  which  produce 
H2.  Thus,  between  reactions  4-5a  and  4-5d,  reaction  4-5d  is  considered 
the  dominant  pathway.  Reaction  4-5e  is  considered  to  be  secondary  due  to 
the  lack  of  any  observed  NH2  emission  (not  a  positive  indicator),  no 
identifiable  NHX  absorption  in  the  long  path  data,  and  the  lack  of  any 
condensable  NHX  compounds  by  Dewhurst  and  Cooper29.  Therefore  the 
primary  reaction  step  for  N-atoms  on  SiH2  is  considered  reaction  4-5 d. 
More  realistically,  this  step  probably  proceeds  through  an  intermediate  so 
the  reaction  becomes 

SiH2  +  N  -  NSiH2*  -  SiN  +  H2  (eq.  4-6) 

It  should  be  pointed  out  that  eq.  4-6  does  not  energetically  allow  the 
formation  of  the  B-state  of  SiN  i.e., 

SiH2  +  N  -♦  SiN(B2X,  v=0)  +  H2  (+5)  (eq.  4-7a) 

-  SiN(B2I,  v=5)  +  H2  (+19)  (eq.  4-7b) 

Thus  the  reaction  given  in  eq.  4-6  cannot  account  for  the  SiN* 
chemiluminescence,  but  the  breakup  of  the  NSiH2*  complex  could  possibly 
give  rise  to  some  broadband  continuum  emission.  Horie  et  al.  have  also 
proposed  reaction  4-5d  as  a  key  reaction  in  the  chain,  however  no  mention 


(-12) 

(+18.4) 

(-64) 

(-27.8) 


(eq.  4-5b) 
(eq.  4-5c) 
(eq.  4-5d) 
(eq.  4-5e) 
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is  made  of  a  possible  complex  as  suggested  in  reaction  4-6.  Reactions 
involving  SiH3  are  probably  occurring  in  the  SiH4  reaction,  but  it  is 
difficult  to  produce  a  scheme  which  does  not  either  produce  H-atoms  or 
some  NHX  compound  which  have  been  previously  ruled  out.  In  addition,  the 
reaction  of  H-atoms  on  S1H3  is  considered  to  occur  at  a  near  gas-kinetic 
rate23,  so  that  any  SiH3  produced  by  the  reaction  given  in  line  two  of  table 
4-2  will  be  converted  to  SiH2  by  the  reaction 

SiH3  +  H  — *  SiH2  +  H2  (eq.  4-8) 

using  the  H-atoms  produced  in  the  initial  SiH4  dissociation.  Thus  one 
concludes  that  SiH3  molecules  do  not  play  a  primary  role  in  the  reaction 
kinetics.  This  does  not  mean  that  any  SiH3  produced  in  SiH4  dissociation 
does  not  react  with  N-atoms.  Reactions  of  the  type 

SiH3  +  N  -SiNH  +  H2  (-112)  (eq.  4-9a) 

-  SiNH3  (?)  (eq.  4-9b) 

may  very  well  take  place,  but  the  assumption  made  here  is  that  they  do 
not  participate  in  the  primary  chain  reaction  that  leads  to  SiH4  dis¬ 
sociation.  Due  to  the  large  exothermieity  of  reaction  4-9a,  there  is  again 
the  possibility  that  the  species  may  be  formed  in  an  excited  state  which 
could  possibly  contribute  to  the  unidentified  chemiluminescence. 

The  Si  and  SiN  Chemiluminescence 

As  stated  earlier,  the  emission  produced  from  Si  and  SiN  is  assumed 
to  originate  only  from  species  containing  Si  and  N-atoms.  The  SiN 
produced  in  reaction  4-6  is  probably  formed  in  a  lower  energy  state  which 
has  been  suggested  to  be  one  of  the  quartet  states  from  spin 
considerations62.  If  the  quartet  states  are  involved,  relaxation  to  the 
ground  state  can  only  occur  through  collisions  with  other  species  in  the 
afterglow.  The  assumption  is  made  that  most  of  the  SiN  is  relaxed  to  the 
ground  state  before  collisions  with  an  N-atom  can  occur.  The  most  likely 
reaction  for  the  formation  of  Si*  is 

SiN  +  N  -*  Si  +  N2  (-122)  (eq.  4-10) 

which  is  exothermic  by  =5.3  eV.  With  the  exception  of  the  lines  at  243.5 
nm  and  263.1  nm,  all  of  the  lines  listed  in  table  3-1  have  upper  levels  that 
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range  in  energy  between  4.13  and  5.08  eV  above  the  ground  state.  For  the 
lines  originating  from  levels  below  5.3  eV,  reaction  4-10  can  supply  the 
energy  for  Si  excitation.  The  Si  emission  is  produced  by 

SiN  +  N  -*  Si*  +  N2  -»  Si  +  hv  +  N2  (eq.  4-11) 

For  lines  originating  from  states  with  energies  above  5.3  eV  another 
mechanism  must  be  proposed.  Two  possible  reactions  may  be  involved  to 
supply  this  excess  energy.  The  first  reaction  involves  the  reaction  of  N- 
atoms  with  vibrationally  excited  SiN,  i.e., 

SiN(X2X,v>0)  +  N  -  Si*  +  N2  (eq.  4-1  la) 

The  vibrational  spacing  of  SiN  is  *0.1  eV  so  that  V=5  would  yield  an 
additional  0.5  eV  for  reaction  4-1  la.  The  line  at  263.1,  however, 

originates  from  a  6.62  eV  level  which  is  very  high  for  vibrationally 

excited  SiN.  The  second  proposed  reaction  to  explain  the  higher  lying  Si 
levels  involves  reactions  of  SiN  with  metastable  N-atoms.  These  are 

SiN  +  N(2D)  -  Si*  +  N2  (eq.  4-1 1b) 

SiN  +  N(2P)  -  Si*  +  N2  (eq.  4-1 1c) 

These  species  supply  an  additional  2.38  eV  and  3.57  eV  respectively  to  the 
reaction.  One  other  possible  reaction  to  explain  the  higher  lying  Si  energy 
levels  may  involve  the  reaction  of  N-atoms  with  an  electronic  metastable 
of  SiN,  but  at  this  time  the  metastable  levels  of  SiN  have  not  been  iden¬ 
tified,  and  thus  this  mechanism  is  considered  too  speculative  even  for  an 
admittedly  very  speculative  part  of  the  thesis.  Horie  et  al.  have  proposed 
collisions  of  Si  with  N2(A3X)  metastables  as  the  mechanism  for  the 
production  of  Si*;  however,  they  require  reaction  4-10  for  the  production 
of  Si  and  there  appears  to  be  no  reason  that  the  reaction  could  not  produce 
Si  in  an  excited  state.  From  the  energetics  of  reaction  4-10,  it  is  clear 
that  the  reverse  reaction  of  Si  on  N2  is  energetic  forbidden.  Horie,  et  al. 
have  proposed  that  Si-atoms  reassociate  with  N-atoms  by  the  process 

Si  +  N  +  M  -♦  SiN  +  M  (eq.  4-12) 

where  M  is  a  third  body.  However,  since  this  reaction  involves  a  three- 
body  process  among  two  species  whose  densities  are  comparatively  low  (N- 
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atoms  are  on  the  order  of  0.01  torr),  this  reaction  will  be  very  slow.  On 
the  other  hand,  a  reaction  of  the  type 

Si  +  N2  +  M  -♦  SiN2  +  M  (eq.  4-13) 

can  take  place  in  the  times  observed  because  it  involves  one  major 
species,  N2.  Horie,  et  al.  also  mention  this  reaction  but  consider  it  only 
secondary  to  reaction  4-12.  It  is  felt  that  a  three-body  reaction  of  this 
type  can  be  the  only  three-body  reaction  that  can  take  place  in  the  neces¬ 
sary  time  scales  needed  to  match  the  observations.  Thus  reaction  4-13  is 
proposed  as  the  next  step  in  the  reaction  chain.  The  next  step  in  the  chain 
is  the  reaction  of  N-atoms  with  the  SiN2  formed  in  reaction  4-13, 
specifically,  the  reaction  is 

SiN2  +  N  -*  SiN*  +  N2  — *  SiN  +  hv  +  N2  (eq.  4-14) 

Unfortunately,  the  heat  of  formation  of  SiN2  is  not  known  so  the  ener¬ 
getics  of  reaction  4-14  can  not  be  determined.  A  consistency  check  can  be 
made  by  comparing  the  energetics  of  reactions  4-13  and  4-14.  From  4-13, 
the  heat  of  formation  of  SiN2  must  be  less  than  that  of  Si,  i.e., 

H,(SiN2)  <  Hf(Si)  =  107.7  (eq.  4-15) 

and  from  equation  4-14, 

H,(SiN2)  >  H,(SiN[B2I,v=5])  -  Hf(N)  =  87.1  (eq.  4-16) 

Therefore,  at  least  the  two  reactions  do  not  violate  any  consistency 
arguments.  By  contrast,  Horie,  et  al.,  who  also  propose  a  reaction  similar 
to  4-14,  suggest  the  reaction 

SiN2  +  N  -♦  SiN  +  N2(A3Z)  (eq.  4-17) 

which  requires  that  the  heat  of  formation  of  SiN2  be 

Hf(SiN2)  >  Hf(SiN)  +  Hf(N2[A3I])  -  Hf(N)  =  146.2  (eq.  4-18) 

This  is  clearly  inconsistent  with  equation  4-13,  which  Horie  et  al.  also 
indicate  can  occur. 
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Reactions  4-11,  4-13,  and  4-14  represent  a  chain  reaction  which 
can  explain  the  observed  Si  and  SiN  chemiluminescence.  The  three-body 
reaction,  4-13,  is  responsible  for  the  observed  pressure  dependence  of  the 
Si  and  SiN  emission  and  the  involvement  of  only  Si-atom  and  N-atom 
bearing  compounds  satisfies  that  requirement.  Furthermore  the  chain 
described  by  these  three  equations  gives  an  almost  identical  time  depen¬ 
dence  to  the  Si  and  SiN  emission  which  is  observed  by  Horie,  et  al. 

Species  Generated  to  Dissociate  Silane 

In  step  3  of  the  general  reaction  scheme,  it  was  suggested  that 
during  some  reaction  with  N-atoms,  a  species  [B]  was  generated  which 
would  dissociate  SihU  and  speed  up  the  chain  reaction.  For  reactions  of 
this  type,  a  favorite  species  of  workers  who  have  studied  active  nitrogen 
is  the  N2(A3X)  metastable  (see  reference  4).  Following  this  line,  Horie,  et 
al.  have  proposed  that  N2(A3X)  metastables  created  in  reaction  4-17  are 
available  to  react  back  and  dissociate  SiH4.  The  difficulties  with  this 
scheme  have  been  discussed  above  but  will  be  reiterated  here.  If  N2(A3Z) 
are  to  be  generated  rapidly,  then  from  equation  4-17,  SiN2  must  also  be 
generated  rapidly.  From  the  energetics  discussed  above,  the  only 
mechanism  available  for  producing  SiN2  in  the  Horie,  et  al.  scheme  is 
reaction  4-12.  But  this  reaction  will  be  very  slow  unless  the  three-body 
rate  is  unusually  large.  Therefore,  it  is  doubtful  that  N2(A3X)  produced  in 
the  scheme  suggested  by  Horie,  et  al.  can  be  responsible  for  the  further 
dissociation  of  SiH4.  Even  if  some  other  method  for  producing  N2(A3Z) 
were  suggested,  it  is  still  doubtful  that  the  metastable  is  playing  a  role. 
As  discussed  by  Caledonia  et  al.47,  one  of  the  significant  loss  mechanisms 
for  N2(A3X)  is  by  the  "pooling"  reaction 

N2(A3I)  +  N2(A3I)  -  N^X1!)  +  N2(C3n)  (eq.  4-19) 

which  occurs  at  a  gas-kinetic  rate.  The  N2(C3ri)  state  quickly  (x=40  ns) 
radiates  down  to  the  N2(B3n),  which  in  turn  radiates  down  to  the  N2(A3X). 
Thus  two  A3E  molecules  react  to  give  back  one.  The  B-A  radiation  is  the 
familiar  first  positive  emission  from  the  afterglow.  Normally,  the  con¬ 
tribution  to  the  observed  first  positive  emission  from  the  pooling  reaction 
is  small  due  to  the  low  metastable  densities.  If,  however,  the  N2(A3I) 
density  is  raised  so  as  to  significantly  dissociate  the  SiH4,  then  this 
increased  metastable  density  would  lead  to  increased  first  positive 
intensity  during  the  reaction.  This  is  not  observed  in  this  work  or  the 
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work  of  Horie  et  al.  The  first  positive  intensity  is  seen  to  decrease 
monotonically  through  the  reaction  zone  after  the  introduction  of  SiH4 
into  the  afterglow.  To  further  clarify  this  point,  the  data  of  Horie,  et  al. 
can  be  used.  Figure  4-2  is  a  reproduction  of  figure  4  from  the  paper  of 
Horie  et  al.  and  shows  the  the  time  dependence  of  N-atoms  and  SiH4  (and 
H2  in  one  case)  during  the  reaction.  If  the  SiH4  is  dissociated  by  some 
species,  [B],  then  the  time  dependence  of  the  SiH4  is  given  by 


d[SiH4] 

dt 


-K-[SiH4][B] 


or 


1  t  d[SiH4] 
[SiH4]  dt 


(eq.  4-20) 


For  example,  measuring  the  slope  of  the  SiH4  curve  from  figure  4-2(c)  at  a 
point  where  the  density  has  gone  to  =0.5  gives  a  value  of  =150  s-1  and 
dividing  by  0.5  gives 

K.[B]  -  300  S"1  <g>  t=25  ms  (eq.  4-21 ) 


Now  using  the  N2(A3X)  quenching  rate  measured  by  Caledonia  and  Piper31 
of  5.9x1 0_1 2  cm3/s,  and  assuming  this  rate  is  also  the  rate  for  dissociat¬ 
ing  SiH4,  gives  a  value  for  [B]  of 


[B]  «  5x1 01 3  cm"3 


which  is  a  more  than  three  orders  of  magnitude  increase  over  the  density 
expected  at  3  torr.  This  is  a  very  high  N2(A3I)  density  and  would  certainly 
perturb  the  first  positive  intensity  in  the  afterglow.  Therefore,  the 
N2(A3I)  must  be  ruled  out  as  the  species  [B]  generated  in  the  afterglow 
which  dissociates  SiH4.  The  N2(aMI)  metastable  might  also  be  considered 
as  a  possible  species  for  dissociating  SiH4  in  the  latter  stages  of  the 
reaction.  The  difficulty  with  this  species  is  that  it  lies  at  8.2eV  above  the 
N2  ground  state,  hence  for  a  reaction  to  produce  this  state  it  must  be 
exothermic  by  >189  Kcal/mole.  No  reactions  could  be  postulated  that  were 
this  exothermic. 
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FIGURE  4-2:  FIGURE  4  OF  REFERENCE  30; 

[HOR1E,  ET  AL.,  CHEM.  PHYS.  LETT..  VOL.  4,  p.  231  (1988)] 


Fig.  4.  Concentration  profiles  of  N,  SiH4  and  H2  in  SiH4  + 
N/N$  ■  The  relative  concentrations  are  [N]/[N]  0,  [SiH4 ]  / 
[SiH4]0,  and  (H2]/2[SiH4) 0.  (a)  p  =  416  Pa.  (N]0=  2.38  x 
1014  cm'3 ,  [SiH4J o  x  4.21  X  1012  cm'3  (m  :  56).  (b)  p  = 
416  Pa,  (N]0  »  2.38  X  1014  cm'3 ,  (SiH4J 0  =  1.04  X  1013 
(m  -  23).  (c)  p  =  390  Pa,  (N)0  =  1.47  X  1014  cm'3 .  [SiH4)0 
-  8.8  X  1012  cm'3  (m  =  17). 


One  of  the  first  species  looked  at  in  this  work  to  fulfill  the  role  of  [B]  was 
atomic  silicon,  thought  to  react  with  SiH4  by  the  reaction, 


SiH4  +  Si  -  2SiH2  (eq.  4-22) 

This  reaction  has  a  rate  constant  of  5.33x1 0-1 3  cm3/s  calculated  by 
Coltrin  et  al.63.  While  the  value  of  this  rate  constant  may  have  been  too 
small  for  the  results  observed  here,  the  uncertainty  in  the  rate  was  large 
enough  to  accommodate  the  necessary  increase.  The  difficulty  with  this 
reaction  is  that  with  the  recent  calculations  by  Ho  et  al.54  and  measure¬ 
ments  by  Boo  and  Armentrout64  the  accepted  value  for  the  heat  of  forma¬ 
tion  of  SiH2  is  now  high  enough  so  as  to  make  reaction  4-22  endothermic 
by  *20  Kcal/mole.  In  light  of  this  recent  data,  it  is  very  doubtful  that 
reaction  4-21  is  exothermic  and  hence  cannot  be  used  to  drive  the  SiH4 
dissociation. 

Another  species  that  might  be  suspected  of  driving  the  SiH4  dis¬ 
sociation  is  atomic  hydrogen.  This  can  be  ruled  out  using  the  same  argu¬ 
ments  used  to  arrive  at  eq.  4-21.  The  rate  for  dissociation  of  SiH4  by  H- 
atoms  by  the  reaction 

SiH4  +  H  -»  SiH3  +  H2  (eq.  4-23) 

has  been  measured  at  2.7x1 0-1 2  cm3/s  (see  reference  23).  Using  eq.  4-21, 
this  would  give  an  H-atom  density  of  *1.1  xIO14  cm-3.  The  [SiH4]  density  at 
t=0  in  figure  4-2(c)  was  only  8.8x1 01 2  cm-3,  hence  the  H2  density 
measurements  made  by  Horie,  et  al.  would  certainly  have  differed  from  the 
2x[SiH4]  result  they  found. 

It  may  be  possible  that  some  other  silicon-bearing  species  formed  in 
the  afterglow  can  dissociate  silane;  however,  since  almost  nothing  is 
known  about  possible  reactions  between  many  of  these  species  (such  as 
SiN  or  SiN2)  and  SiH4,  the  search  for  the  species  responsible  for  SiH4 
dissociation  returned  to  the  active  nitrogen  itself.  Earlier,  the  exother- 
micity  of  reactions  4-10  and  4-14  was  used  to  excite  Si-atoms  and  SiN  to 
explain  the  chemiluminescence  from  these  species.  It  is  also  possible  that 
these  reactions  may  leave  the  N2  product  vibrationally  excited,  i.e., 

SiN  +  N  -  Si  +  N^X1!,  v>0)  (eq.  4-24) 


and 
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SiN2  +  N  -  SiN  +  N^X1!,  v>0) 


(eq.  4-25) 


From  the  energetics  of  reaction  4-24,  there  is  sufficient  excess  energy  to 
give  vmax=20,  and  since  it  is  assumed  that  reaction  4-14  is  capable  of 
exciting  the  B2X,  V=5  state  of  SiN,  this  yields  a  value  of  vmax-13.  As 
discussed  earlier,  vibrationally  excited  N2(X1X)  at  v>10  possesses  enough 
internal  energy  to  dissociate  SiFI4;  therefore,  vibrationally  excited 
nitrogen  produced  in  either  reaction  4-24  or  4-25  can  dissociate  SiFi4.  The 
difficulty  with  proposing  this  species  as  a  driver  for  SiFI4  dissociation 
lies  in  the  size  of  typical  rate  constants  for  v  to  E  transfer.  As  pointed 
out  by  Rich  and  Bergman61,  rates  for  v  to  E  transfer  typically  range  from 
10~3  to  10~4  of  gas  kinetic.  Using  a  rate  of  2.7x10'13  in  equation  4-21 
gives  a  value  for  [B]  of 

[B]  =  l.lxlO15  cm3 

at  a  total  pressure  of  3  torr.  This  would  represent  a  fractional  vibrational 
population  for  v>10  of  =1%.  This  is  a  very  high  fractional  population,  and 
it  is  very  doubtful  that  this  could  ever  be  realized  in  the  afterglow.  On  the 
other  hand,  there  have  not  been  any  measurements  of  v  to  E  transfer  to  a 
molecule  with  as  many  degrees  of  freedom  as  SiFt4,  thus  it  is  possible  that 
the  rate  in  this  case  might  be  larger  than  previously  measured  rates.  With 
a  larger  v  to  E  rate,  a  lower  vibrational  population  would  be  sufficient  to 
dissociate  SiH4  at  the  rate  observed  by  Fiorie,  et  al.  It  is  therefore 
proposed  that  the  species  [B],  generated  in  the  reaction,  which 
is  responsible  for  SiH4  dissociation  is  vibrationally  excited, 
ground  electronic  state  N2. 

Summary  of  the  Proposed  Reaction  Scheme 

The  proposed  reaction  scheme  is  summarized  below: 


N2(X1I,v>10)  +  SiFI4(X1A1)  ^  SiH2(X1A1)  +  H^X1!)  +  N^X1!)  (R-1 ) 

N2(A3I)  +  SiFI4(X1A1)  ^  SiH2(X1  A^  +  H^X1!)  +  N^X1!)  (R-2) 

SiFI2(X1  Ai)  +  N(4S)  -  SiN(4I)  +  HpCX1!)  (R-3) 

SiN(4I)  +  M  SiN(X2X)  +  M  (R-4) 
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(R-5) 


SiN(X2I)  +  N(4S)  -  Si(3p2  3P)  +  N2(X1I,v<20) 


SiN(X2I)  +  N(4S)  -  Si(4s  3P)  +  ^(X1!)  (R-6) 

Si(4s  3P)  Si(3p2  3P)  +  hv  (R-7) 

Si(3p2  3P)  +  N2(X1I)  +  M  -  SiN2(?)  +  M  (R-8) 

SiN2(?)  +  N(4S)  SiN(X2I)  +  N^I.v^lS)  (R-9) 

SiN2(?)  +  N(4S)  -»  SiN(B2I)  +  N^X1^)  (R-1 0) 

SiN(B2I)  -  SiNfX1!)  +  hv  (R-1 1 ) 


While  it  is  almost  certain  that  additional  reactions,  besides  the  ones 
listed  here,  are  also  occurring,  the  reaction  scheme  presented  here  seems 
capable  of  describing  qualitatively  the  bulk  of  the  observed  results.  For 
simplicity,  the  above  scheme  only  uses  N^X^.v^lO)  and  N2(A3Z)  as  the 
initial  drivers  for  SiFU  dissociation.  Since  it  is  assumed  this  initial 
process  is  slow,  there  is  nothing  to  be  gained  by  including  additional 
species.  Also  as  previously  discussed,  reactions  initially  involving  S1H3 
are  not  considered  as  primary  to  the  reaction  scheme.  Figure  4-3  is  a 
schematic  of  the  reaction  mechanism  showing  the  feedback  loops  that 
exist  in  the  scheme.  The  chemiluminescence  occurs  in  reactions  R-7  and  R- 
11  and  is  generated  in  a  loop  composed  of  reactions  R-6,  R-7,  R-8,  R-10, 
and  R-1 1 .  The  feedback  to  the  SiFU  dissociation  is  supplied  by  reactions  R- 
5  and  R-9.  All  of  these  reactions  are  driven  by  N-atoms  yet,  as  required, 
there  is  no  direct  reaction  of  N-atoms  on  SiF^.  Before  the  scheme  can  be 
quantitatively  assessed,  the  mathematical  model  described  at  the  begin¬ 
ning  of  this  chapter  must  be  constructed.  The  next  chapter  will  describe 
that  model. 
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Chapter  V 

The  Reaction  Model 

As  discussed  in  chapter  IV,  the  reaction  is  modeled  by  constructing 
a  set  of  coupled  rate  equations,  solving  these  equations,  comparing  the 
calculations  with  observations,  and  readjusting  the  model  if  necessary. 
Over  the  course  of  this  work,  a  large  number  of  models,  based  on  various 
proposed  reaction  schemes,  have  been  constructed.  While  these  models 
were  very  valuable  for  guidance  towards  a  "working"  reaction  mechanism, 
it  was  also  dominated  by  long  and,  for  the  most  part,  frustrating  excur¬ 
sions  down  dead  end  streets.  In  light  of  the  questionable  value  of  describ¬ 
ing  in  detail  the  many  failures,  this  chapter  will  concentrate  on  the  model 
developed  from  the  reaction  scheme  proposed  in  chapter  IV.  Some  back¬ 
ground  and  a  brief  summary  of  some  of  the  early  modelling  efforts  will  be 
described  below. 

Early  MQdgi? 

Initially,  the  modelling  effort  was  to  include  a  model  of  the  nitrogen 
discharge,  which  would  supply  input  to  an  afterglow  model,  which  would 
be  coupled  to  the  silane  reaction  scheme.  The  early  discharge  model  and 
afterglow  model  were  based  on  the  work  of  Slater,  et  al.65  and  modified 
for  the  low  pressure  dc  discharge  used  in  this  work.  It  was  quickly  dis¬ 
covered  that  the  discharge  model  did  not  yield  excited  state  densities 
consistent  with  measured  values  in  a  dc  discharge66  and  measurements  in 
the  early  afterglow47.  For  example,  model  predictions  of  the  N-atom 
density  were  as  high  as  10%  while  measurements  by  Caledonia,  et  al.47 
give  typical  densities  of  less  than  1%.  The  densities  were  also  a  strong 
function  of  the  residence  time  of  the  N2  in  the  discharge,  and  steady  state 
values  were  not  reached  for  very  long  times  (>10  s).  It  was  suspected  that 
the  electron  impact  dissociation  rate  of  N2  suggested  by  Slater  et.  al.65 
was  too  large  to  accurately  model  the  data.  During  this  discharge  modell¬ 
ing  effort,  comparisons  were  attempted  between  discharge  parameters 
and  quantities  measured  in  the  reactor.  This  proved  to  be  fruitless,  and 
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the  discharge  model  was  abandoned.  For  input  into  the  afterglow  model, 
excited  state  densities  measured  in  early  afterglows  were  taken  as  the 
initial  conditions  and  the  afterglow  model  was  allowed  to  "relax"  for 
approximately  30ms  before  introducing  the  silane.  The  afterglow  model 
contained  the  reactions  used  by  Slater  et  al.65  with  updated  rates  from 
Caledonia  et  al.47.  It  was  during  this  part  of  the  modelling  that  the  sig¬ 
nificance  of  the  pooling  reaction  to  the  first  positive  intensity  was 
clearly  identified  and  was  key  in  eliminating  the  N2(A3X)  metastable  from 
consideration  as  the  species  driving  the  silane  dissociation.  The  next 
modelling  effort  involved  adding  a  complex  model  of  the  neutral 
silane/hydrogen  chemistry  developed  by  Kushner23.  This  model  was  very 
complex  and  contained  as  many  possible  reactions  between  silane  and 
active  nitrogen  species  as  could  be  reasonably  imagined.  The  complexity 
of  this  model  and  the  many  unknown  rates  made  evaluation  of  the  various 
processes  in  the  model  very  difficult.  In  spite  of  this,  it  was  during  this 
modelling  effort  that  the  significance  of  the  atomic  silicon-silane 
reaction  discussed  in  chapter  IV  was  identified.  Inclusion  of  this  reaction 
in  the  model  gave  the  first  indications  of  the  chain  reaction  behavior  seen 
in  the  observations.  When  this  reaction  was  eliminated  from  consideration 
due  to  the  reevaluation  of  the  heat  of  formation  data  for  SiH2,  the  complex 
model  was  completely  scrapped.  The  reactions  that  appeared  "key"  in  the 
complex  model  were  pulled  out,  along  with  a  very  basic  afterglow  model, 
and  a  simpler  model  was  begun.  The  model  discussed  in  the  remainder  of 
the  chapter  represents  the  evolution  of  this  simple  model. 

The  Afterglow  Model 

The  afterglow  model  was  made  as  simple  as  possible  while  main¬ 
taining  the  most  important  reactions  among  the  dominate  species  present 
in  the  afterglow.  The  species  included  in  the  afterglow  model  were  N(4S), 
N(2D),  and  N2(A3X).  The  N(2P)  metastable  was  not  included  because  of  a 
lack  of  good  rate  data  on  this  species.  To  investigate  the  effects  of 
metastable  N(2P),  the  N(2D)  was,  at  times,  treated  as  "typical"  and  allowed 
to  dissociate  SiPU  as  described  in  chapter  IV.  However,  in  the  final  model 
this  was  not  the  case.  Vibrationally  excited  N^X1!)  was  included  as  a 
separate  set  of  species  and  will  be  discussed  later.  The  C3n  and  B3n  state 
were  left  out  of  the  model  once  the  dependence  of  the  first  positive 
intensity  had  been  established  in  the  more  complex  model.  Since  the  first 
positive  intensity  is  only  dependent  on  the  N-atom  density  through  recom¬ 
bination  and  the  A3Z  density  through  the  pooling  reaction,  the  densities  of 


80 


these  species  were  enough  to  predict  the  first  positive  intensity.  Since 
the  N-atom  recombination  is  assumed  to  proceed  into  the  B3n  state,  which 
radiates  down  to  the  A3I,  for  the  purposes  of  the  model  it  was  assumed 
that  the  N-atom  recombination  proceeds  directly  to  the  A3X  state.  Thus 
for  this  very  simple  afterglow,  the  following  reactions  were  used  (the 
equation  numbering  continues  from  chapter  IV) 


N(4S)  +  N(4S)  +  M  -*  N2(A3I)  +  M 

(R-12) 

N2(A3I)  +  N2(A3I)  -  N2(A3I)  +  N2(X1I) 

(R-1 3) 

N2(A3I)  +  N(4S)  -  N2(X1I)  +  N(2D) 

(R-14) 

N(2D)  +  N(4S)  -  N(4S)  +  N(4S) 

(R-i  5) 

N(2D)  +  N2(X1X)  -  N(4S)  +  N2(X1I) 

(R-1 6) 

Rates  for  these  reactions  were  taken  from  references  47  and  65.  The  rates 
used  for  the  above  reactions  were: 

R-12:  ki2  =  2.33x1 0"33  cm6/s 

R-1 3:  k13  =  I.OxlO-10  cm3/s 

R-14:  ki4  =  5.0x10-11  cm3/s 

R-1 5:  k15  =  I.OxlO-11  cm3/s 

R-1 6:  k16  =  2.0x1 0-1 6  cm3/s 

The  initial  conditions  used  at  the  point  of  injection  of  the  silane  were 
slightly  different  from  table  4-1.  The  N(4S)  density  could  be  entered  as  a 
parameter  in  the  model,  and  the  N2(A3Z)  and  N(2D)  densities  were  set  as 
4.3x1 0-6  and  2.4x1 0-5  times  the  N(4$)  density.  This  simple  model  supplied 
realistic  densities  of  N-atoms  and  N2(A3I)  metastables  which  were  the 
primary  drivers  (plus  N2(X1Z,v>0))  in  the  silane  reaction  model. 

Vibrationallv  Excited  N^fX1!) 

In  the  first  stages  of  this  latter  modelling  effort,  vibrationally 
excited,  ground  electronic  state  N2  was  simply  treated  as  a  single  species 
which  was  generated  by  reactions  R-5  and  R-9  of  chapter  IV  and  dis¬ 
sociated  SiH4  by  reaction  R-1.  While  this  simple  treatment  showed  the 
necessary  characteristics  of  a  chain  reaction,  a  more  realistic  treatment 
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was  needed  to  answer  questions  about  vibrational  relaxation  and/or 
vibrational  pumping.  The  method  chosen  to  treat  the  vibrationally  excited 
molecules  was  to  treat  each  vibrational  level  as  a  different  species  in  the 
model  and  supply  the  rates  for  vibrational  exchange  and  relaxation.  This 
method  is  often  referred  to  as  the  Master  Equation  Approach.  The  treat¬ 
ment  to  be  discussed  here  follows  very  closely  the  work  of  Bailey67  and 
also  the  work  of  Capitelli68.  The  equation  governing  the  density  of  a 
species  in  some  energy  level  I  is  given  by 

— -  =  X  Rj . iNjNk  -  X  Ri.iNiNi 

°l  i.j.k  i.j.k  (eq.  5-1) 


~i,k 

where  Ni  is  the  density,  and  the  Ri.j's  represent  generalized  rate  coeffi¬ 
cients  for  the  population  and  depopulation  of  state  I.  The  above  set  of 
equations  are  commonly  referred  to  as  the  Master  Equation.  This  set  of 
equations  is  very  difficult  to  solve  due  to  the  many  states  involved  so  the 
following  simplifying  assumptions  are  made  (see  Bailey,  reference  67): 

1.  The  translational  and  rotational  degrees  of  freedom  are  in  com¬ 
plete  Boltzmann  equilibrium. 

2.  The  molecules  in  each  vibrational  state  are  treated  as  separate 
species. 

3.  Only  single  quantum  vibrational  exchanges  occur. 

4.  The  molecule  is  represented  as  a  Morse  oscillator  rigid  rotator 
system. 

With  these  assumptions  the  Master  Equation  becomes 


dNv 

dt 


=  NWJ.iR 


v  +  1r*v  +  1  ,v 


+  Nw.iR 


v-1nv-1 


,v-  Nv[R 


v.v+1 


(eq.  5-2) 


where  the  Rg's  are  the  rates  for  change  of  vibrational  quanta  from  i  to  j. 
The  Rjj’s  can  be  written 


v  - 1  k  k+1 

Rv  +  1  ,v  =  N(VT)v„.v  +  2(VV)v;,i#Nk 


k=0 


(eq.  5-3) 
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(eq.  5-4) 


Rv.,.v  =  N-(VTB)v.1iV  +  ^(VV)ki'K.,,VNk 


k=1 


where  the  (VT)j,j*s  are  the  rates  for  the  loss  of  a  quantum  of  vibration  to 
translation  (vibration-to-transiation,  VT)  and  the  (VTB)j,j’s  are  the  rates 
for  the  gain  of  a  quantum  of  vibration  from  the  translational  energy  pool 
(VT  Back,  VTB).  The  VV's  are  the  rates  for  exchange  of  vibrational  quanta 
between  molecules  in  different  vibrational  states,  N  is  the  density  of 
colliding  species,  and  v*  is  the  maximum  vibrational  level  considered.  If 
the  VT  and  VV  rates  are  known,  the  set  of  equations  given  in  5-2  can  be 
integrated  from  some  initial  conditions  to  give  the  populations  in  each 
vibrational  level  as  a  function  of  time.  In  general,  the  VV  rates  only 
include  vibrational  exchange  among  like  molecules,  thus  in  this  work  only 
exchanges  between  two  N2  molecules  are  considered.  On  the  other  hand, 
the  VT  rates  are  strongly  affected  by  all  species  that  collide  with  the 
molecule.  In  this  experiment,  the  dominant  species  was  always  N2;  there¬ 
fore,  for  purposes  of  the  model  the  only  VT  rates  considered  were  those 
involving  N2-N2  collisions.  It  should  be  pointed  out  that  in  cases  where 
other  major  species  are  present,  the  VT  rates  for  collisions  between  that 
species  and  N2  must  also  be  included. 

For  most  molecules,  the  VV  and  VT  rates  are  only  known  for  the  very 
lowest  lying  vibrational  levels.  To  accurately  model  the  behavior  of  the 
ground  state  vibrational  manifold,  it  is  necessary  to  include  a  large 
number  of  levels,  thus  rates  for  high  lying  vibrational  states  are  needed. 
The  general  procedure  to  obtain  these  rates  is  to  extrapolate  the  known 
rates  from  the  lowest  levels  to  the  higher  states  using  a  theoretical 
expression  for  the  scaling  of  the  rates  with  vibrational  quantum  number. 
The  theory  usually  used  for  this  scaling  the  rates  is  due  to  Schwartz, 
Slawsky,  and  Herzfeld69  and  is  commonly  referred  to  as  the  SSH  theory.  In 
the  SSH  theory,  the  VT  transition  probability  (probability  per  collision) 
for  a  molecule  to  go  from  v+1  to  v  is  given  by 

Py+i  ,v  “  P(T)-(v+1)/[1-8-(v+1  )]  F(Yv+1  v)  (eq.  5-5) 

where  8  is  the  anharmonicity  of  the  molecule,  P(T)  is  a  parameter  related 
to  the  1-0  relaxation,  T  is  the  temperature,  and  Yv+i,v  is  given  by 
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Yv+i>v=(0.5)3/2Ve7T-[l-2-8-(v+1)] 


(eq.  5-6) 


where  6  is  a  parameter  related  to  the  molecule  (for  N2  Bray70  gives  a 
value  of  8=5.39x1 06  °K).  The  function  F  cannot  be  easily  evaluated  and  is 
usually  approximated  by  (see  reference  70) 

F(y)  -  0.5-[3-Exp(-2/3-y)]-Exp(-2/3-y)  (eq.  5-7) 

The  anharmonicity  of  the  molecule,  8,  is  defined  so  that  the  vibrational 
energy  levels  of  the  molecule  are  given  by 

Ev  =  E10-[(v+i)-5  (v  +  i)2]  (eq.  5-8) 

and  for  N2,  Ei0=2358.57  cm-1  and  8=6.1 27x1 0-3  (from  reference  55).  Using 
the  transition  probability,  the  VT  rate  is  given  by 

(VT)v+i,v  *  Z.Pv+i.v  (eq.  5-9) 

where  Z  is  the  gas  kinetic  collision  rate  calculated  using  equation  4-3 
from  chapter  IV  (for  N2,  the  molecular  diameter=  3.35A).  For  the  function 
P(T),  Bray70  gives  the  form 

P(T)  =  1 .07x1 0-4  *T  (eq.  6-10) 

however,  for  this  work  the  function  described  in  appendix  C  was  used  for 
P(T).  The  (VTB)  rates  were  calculated  from  detailed  balance67, 

(VTB)v,v+i  =  (VT)v+i,vExp[-(Ev+1-Ev)/(kT)]  (eq.  5-11) 

where  the  E's  are  calculated  from  equation  5-8. 

The  VV  rates  are  also  calculated  from  the  SSH  theory.  The  transition 
probability  (per  collision)  for  the  target  molecule  to  go  from  v  to  v+1 
while  the  colliding  molecule  goes  from  j+1  to  j  is  given  by67 

Pv,j  =  SViJ<T)-Exp{[(Ev+1-Ev)-(Ej+1-Ej)]/(2kT)}  (eq.  5-12) 


where  the  E's  are  from  equation  5-8  and  Sv,j(T)  is  given  by 
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(eq.  5-13) 


S  =  Q(TV 

V’J  [1-8*(v+1]  [1-5-(j+1)] 


F(Yv.i) 


where  F(Y)  is  given  by  equation  5-7  and  YVlj  is  given  by 


Yv.j.A.Vi7T.|v.j| 

where  0  is  given  above.  Q(T)  is  given  by  Bray70  for  N2  as 
Q(T)  =  3.7x1  (H-T 


(eq.  5-14) 


(eq.  5-15) 


As  in  equation  5-9  for  the  VT  rates,  the  VV  rate  is  given  by  the  product  of 
the  transition  probability  given  by  equation  5-12  and  the  gas  kinetic  rate 
given  in  equation  4-3.  Again  using  detailed  balance,  the  reverse  rates  are 
given  by67 


(VV)r;T.V=(VV)Vw+V1+i1w-Exp{-[(Ev 

+rEv)-(Ew+rEw)]/kT}  (eq.  5-16) 

Using  equations  5-5  through  5-16  and  appendix  C,  VT  and  VV  rates  were 
calculated  for  T=309°K.  Some  typical  rates  are  shown  in  figure  5-1.  One 
remaining  parameter  was  to  set  the  highest  vibrational  level  to  be 
included  in  the  model.  As  mentioned  in  chapter  IV,  the  N2  ground  electronic 
state  can  support  approximately  70  vibrational  levels,  but  due  to  the 
simple  expression  for  the  energy  of  the  levels  given  in  5-8,  the  number  of 
levels  should  be  set  lower.  In  his  work,  Bray70  used  a  value  of  vmax=47; 
therefore,  this  value  was  also  chosen  in  this  work.  Thus  a  total  of  48 
vibrational  levels  were  included  in  the  set  of  equations  given  by  eq.  5-2, 
along  with  the  calculated  VT  and  VV  rates  to  describe  the  interaction  of 
the  N2(X1X)  vibrational  manifold. 

Finally,  a  set  of  initial  vibrational  populations  was  needed  as  input 
into  the  model.  To  arrive  at  this  set,  the  following  procedure  was  used.  A 
model  was  constructed  which  included  only  the  vibrationally  excited 
N2(X1I)  species  and  the  VT  and  VV  rates.  The  initial  vibrational  popula¬ 
tions  were  set  to  a  Boltzmann  distribution  with  a  vibrational  temperature 
of  1500°K.  This  distribution  was  then  allowed  to  "relax"  for  1  second  at 
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Figure  5-1:  Selected  V-T  and  V-V  Rates  for  Nitrogen  at  300°K 


3  torr.  The  resulting  vibrational  populations  are  shown  in  figure  5-2, 
which  is  a  plot  of  the  fraction  of  molecules  in  a  given  level  vs.  the  vibra¬ 
tional  quantum  number.  As  can  be  seen,  the  distribution  is  "over  populated" 
with  respect  to  a  Boltzmann  distribution  at  thj  higher  levels,  but  the 
lower  levels  are  still  approximately  Boltzmann  at  the  initial  1500°K. 

Since  this  vibrational  temperature  is  typical  of  measurements  made  on 
the  lowest  levels  in  a  nitrogen  discharge71,  this  distribution  was  taken  as 
the  initial  distribution  for  all  runs.  It  should  be  pointed  out  that  one  would 
expect  the  distribution  to  be  affected  by  pressure,  but  since  only  the 
levels  above  v=10  interacted  with  silane,  details  of  the  distribution  were 
unimportant.  The  advantage  of  using  this  distribution  instead  of  a 
Boltzmann  distribution  was  that  this  distribution  was  relatively  stable 
for  times  on  the  order  of  1  s,  in  the  absence  of  any  energy  input,  while  the 
upper  levels  of  the  Boltzmann  distribution  tended  to  "come  up"  due  to  v-v 
pumping.  Therefore,  using  the  distribution  in  figure  5-2  allowed  changes 
due  to  chemical  kinetics  to  be  more  easily  seen. 


The  last  step  in  the  model  was  to  integrate  the  sila”°  reactions  R-1 
through  R-1 1  given  at  the  end  of  chapter  IV  into  the  model.  As  discussed 
in  chapter  IV,  the  rate  constants  can  not  exceed  the  gas  kinetic  rate  for  a 
given  reaction.  When  literature  values  of  molecular  diameters  were 
available,  equation  4-3  was  u^ed  to  set  the  limiting  rate  for  the  model.  In 
cases  where  molecular  diameters  were  not  available,  a  gas  kinetic  rate  of 
IxlO-10  cm3/s  was  chosen  as  a  typical  value.  It  should  be  pointed  out, 
however,  that  the  sum  of  rates  for  reactions  R-5  and  R-6  cannot  exceed 
the  gas  kinetic  rate,  since  the  two  only  yield  different  products  from  the 
same  reaction.  The  same  is  true  for  reactions  R-9  and  R-10.  The  v  to  E 
transfer  rate,  R-1,  was  never  allowed  to  reach  the  gas  kinetic  limit,  but 
was  limited  to  10~2  times  the  gas  kinetic.  The  rate  for  Sif-U  dissociation 
by  the  N2(A3X),  reaction  R-2,  was  taken  from  the  work  of  reference  31  as 
k2=5.9x10-12  cm3/s.  The  rates  R-7  and  R-1 1 ,  which  are  radiative  rates 
were  taken  from  the  literature  when  possible,  for  R-7,  an  average  rate  of 
IxlO7  s_1  was  taken  for  the  states  of  Si.  For  R-1 1 ,  the  recent  measure¬ 
ment  of  the  lifetime  of  the  E2I  of  SiN  of  210±10  ns  by  Walkup,  et  al.72 
was  used  to  derive  the  rate  of  4.8x1 06  s_1.  An  estimate  of  the  gas  kinetic 
three-body  rate  given  in  reaction  R-8  is  difficult.  If  the  process  is  a  true 
three-body  reaction,  requiring  the  three  bodies  to  be  present  within  some 
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Figure  5-2:  Calculated  Fractional  Vibrational  Populations  of  N2  in  the  Afterglow 


small  volume  for  the  reaction  to  occur,  then  a  calculation  based  on  the 
Thompson  theory53  can  be  used.  On  the  other  hand,  if  the  reaction  is 
typical  of  an  association  reaction,  an  intermediate  state  of  some  finite 
lifetime  is  probably  involved  which  must  be  stabilized  by  collision  with  a 
third  body.  Depending  on  the  lifetime  of  this  state,  the  three-body  rate  in 
this  latter  case  can  be  significantly  larger  than  that  given  by  the 
Thompson  theory.  In  the  Thompson  theory53,  the  two  reacting  species  and 
the  third  body  are  assumed  to  occupy  a  sphere  of  radius  r0  for  the  recom¬ 
bination  to  occur.  The  three-body  rate  scales  as  r03  so  that  uncertainties 
in  choosing  this  radius  can  affect  the  rate  dramatically.  Using  available 
data  for  the  diameters  of  N2  and  Si  and  assuming  the  interaction  radius  r0 
was  equal  to  the  average  diameter  of  the  two  gave  a  value  of  the  rate 
constant  of  reaction  R-8  as 

k8  =  6.6x1 0~32  cm6/s  (Thompson  theory) 

Because  of  the  nature  of  the  reaction,  the  above  rate  was  only  used  as  a 
guideline  in  the  model.  For  this  rate,  the  rate  constant  used  in  the  model 
was  allowed  to  exceed  the  above  rate  if  necessary.  For  the  vibrational 
pumping  reactions  R-5  and  R-9,  it  was  assumed  that  the  reactions  left  the 
molecules  in  v=20  and  v=13  respectively.  While  there  is  most  certainly 
some  distribution  of  vibrational  levels  that  result  from  these  reactions, 
with  no  knowledge  of  this  distribution  available,  this  simple  approach 
was  taken.  For  similar  reasons,  it  was  assumed  that  the  v-to-E  transfer 
reaction,  R-1,  could  occur  with  v>10  and  the  reaction  left  the  N2  molecule 
in  v=0.  Again,  there  is  probably  some  distribution  of  final  states  depend¬ 
ing  on  the  initial  vibrational  level  involved  in  the  reaction,  but  this 
simpler  approach  was  taken.  Diffusion  was  not  included  in  the  model 
although  it  almost  certainly  plays  a  role  in  the  reaction.  Changes  in  the 
behavior  of  the  reaction  as  the  tube  walls  became  coated,  as  discussed  in 
chapter  III,  indicated  that  diffusion  of  some  species  to  the  wall  affected 
the  reaction.  Since  the  proposed  reaction  scheme  involves  a  number  of 
feedback  loops,  there  are  a  number  of  steps  where  loss  of  a  species  to  the 
wall  (sticking)  could  change  the  behavior  of  the  reaction.  The  reason  for 
£LQi  including  diffusion  in  the  model  is  that  very  little  is  known  about  the 
sticking  coefficients  of  the  various  species  in  the  model,  hence  diffusion 
would  simply  add  more  uncertainty  to  the  modelling  effort. 

The  model  as  described  above  contained  1 1  "reactive"  species  plus 
48  vibrational  levels  of  N2  for  a  total  of  59  coupled  first  order  differen¬ 
tial  equations.  The  number  of  VV  and  VT  rates  needed  in  the  model  goes  as 
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2-(vtotaiMvtotai+1)  or  4704  rates  for  the  48  levels  used.  Finally,  the  16 
rates  used  to  describe  the  afterglow  and  the  silane  reaction  were  also 
included.  The  differential  equations  were  integrated  using  the  method  of 
Gear73  which  is  specifically  designed  to  handle  sets  of  equations  which 
exhibit  "stiff  behavior"74.  A  detailed  discussion  of  Gear’s  method  is  not 
appropriate  here,  but  it  should  be  pointed  out  that  for  stiff  sets  of 
equations,  use  of  this  method  is  almost  mandatory.  For  example,  it  was 
found  by  the  author,  on  a  much  simpler  set  of  stiff  equations,  that  Gear's 
method  was  a  factor  of  1000  times  faster  than  an  Adams-Bashforth 
method  for  integrating  the  equations.  This  is  a  significant  savings  in 
computer  time  and  as  a  result  of  this  savings,  the  calculations  presented 
here  could  be  performed  on  a  desktop  computer. 

General  Behavior  of  the  Model 


Detailed  quantitative  predictions  of  the  model  were  difficult  to  use 
in  assessing  its  validity  because  of  the  lack  of  specific  rate  data  and  the 
lack  of  detailed  quantitative  measurements  of  the  reaction.  However,  the 
general  behavior  of  the  model  could  be  used  to  evaluate  the  validity  of  the 
overall  scheme.  As  discussed  in  chapter  III,  the  Si  and  SiN  intensities 
increased  faster  than  linear  with  pressure  even  when  the  N-atom  density 
was  taken  into  account.  In  figures  3-1  and  3-4  to  3-7,  the  first  positive 
intensity  increased  roughly  as  the  square  of  the  pressure  while  the  Si  and 
SiN  intensities  increased  even  faster.  Since  the  first  positive  intensity 
goes  as  the  square  of  the  N-atom  density,  figures  3-1  and  3-4  through  3-7 
represent  data  where  the  N-atom  concentration  is  approximately 
constant.  To  model  these  data,  a  series  of  calculations  were  made  at  four 
different  pressures  with  both  the  N-atom  and  SiH4  concentrations  held 
constant.  The  rates  chosen  were  mostly  gas  kinetic  and  are  given  below: 


k1=2.7x10-12 
k4=1x10-10 
k7=1  xl 07 
k10=3.3x1 0_1 1 


k2=5.9x1 0-12 
k5=6.6x10-11 
k8=6.6x1 0"32 
kn=4.8x106 


k3=1x10-10 
k6=3.3x10-11 
k9=6.6x1 0_1 1 


Figure  5-3  shows  the  SiN(B2£)  density  as  a  function  of  the  reaction  time 
(time  after  the  silane  injection)  for  these  conditions.  The  Si  excited-state 
density  followed  a  similar  set  of  curves.  As  shown  in  the  figure,  the  SiH4 
density  was  set  to  0.1%  of  the  total  density  while  the  N-atom  density  was 
set  to  0.2%  of  the  total  density.  Since  the  intensity  of  the  emission  is 
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proportional  to  the  excited  state  density,  these  curves  are  proportional  to 
the  chemiluminescence  intensity.  In  figure  5-3,  the  ratio  of  the  peak 
density  at  10  torr  to  that  at  1  torr  is  351  while  the  calculated  first 
positive  intensity  increases  by  a  factor  of  100.  This  is  similar  to  the 
increases  observed  in  the  emission  spectra.  In  addition  to  faster-than- 
linear  increase  in  the  SiN*  emission  predicted  by  the  model,  the  model 
also  predicts  a  longer  time  for  the  occurrence  of  the  peak  emission  and 
the  spreading  in  time  of  the  chemiluminescence.  Both  of  these  effects 
were  observed  in  the  reaction  flame.  The  exact  position  in  time  of  the 
peak  and  the  ratio  of  the  peak  intensities  at  different  pressures  is  a 
function  of  the  specific  rates  put  into  the  model,  but  the  general  behavior 
remains  the  same.  It  was  also  observed  that  the  peak  intensity  could  be 
made  to  move  back  and  forth  in  the  reaction  zone  by  varying  the  SiH4  flow 
(and  to  a  lesser  extent  the  discharge  current).  This  is  also  predicted  by 
the  model  as  shown  in  figure  5-4.  In  this  case,  both  the  pressure  and  N- 
atom  density  were  held  constant  while  the  SiH4  density  was  varied.  As  can 
be  seen  in  figure  5-4,  as  the  SiH4  density  is  decreased  the  peak  in  the 
emission  occurs  at  later  times.  A  similar  (but  less  sensitive)  effect  is 
also  observed  for  varying  the  N-atom  density.  Since  the  occurrence  of  the 
emission  at  a  later  time  corresponds  to  the  chemiluminescence  occurring 
farther  downstream  of  the  injection  port,  this  agrees  with  the  observed 
movement  of  the  flame  in  the  experiments.  As  before,  the  specific  loca¬ 
tion  of  the  peak  depends  on  the  rates  chosen  in  the  model,  but  the  general 
behavior  remains  the  same. 

The  overall  design  of  the  model  insures  that  the  SiH4  dissociated 
will  be  proportional  to  the  N-atom  density  if  a  "titration"  experiment  is 
performed  as  described  in  chapter  III.  On  the  other  hand,  the  resulting 
proportionality  constant  depends  on  the  rates  put  into  the  model.  To  better 
understand  this,  consider  the  four  reactions  R-5,  R-6,  R-9,  and  R-10. 
Reactions  R-5  and  R-9  each  result  in  the  loss  of  an  N-atom  but  produce  a 
vibrationally  excited  N2  molecule  which  can  dissociate  an  SiH4  molecule 
by  reaction  R-1.  On  the  other  hand,  reactions  R-6  and  R-10  also  result  in 
the  loss  of  an  N-atom,  but  in  this  case  no  SiH4  dissociation  results.  Thus 
it  can  be  argued  that  more  SiH4  will  be  dissociated  per  N-atom  if  the  ratio 
(k5+k9)/(k6+kio)  is  increased.  This  is,  in  fact,  observed  in  the  model, 
but  the  functional  dependence  of  the  SiH4  dissociation  is  much  more 
complex  than  the  previous  expression  would  lead  one  to  believe. 
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CALCULATED  TIME  DEPENDENCE  OF  SiN  EXCITED  STATE 
DENSITY  AS  A  FUNCTION  OF  SiH.  DENSITY 
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Figure  5-4:  Calculated  Time  Dependence  of  SiN  Excited  State  Density  as  a  Function  of  SiFU  Density 


Since  the  time  dependence  of  the  various  species  in  the  model  are 
strongly  affected  by  the  magnitudes  of  the  rates  used,  time  dependent 
measurements  of  the  reactive  species  could  further  validate  the  reaction 
model.  As  discussed  earlier,  and  as  shown  in  figure  4-2,  Horie,  et  al. 
(reference  30)  have  made  measurements  of  the  time  dependence  of  the 
SiH4,  the  N-atom,  and  (in  one  case)  the  H2  density  in  the  reactor  at  ap¬ 
proximately  3  torr.  Since  these  measurements  could  not  be  made  in  this 
experiment  (Horie,  et  al.  used  a  mass  spectrometer),  it  was  decided  to  try 
and  model  their  results.  To  allow  comparison  of  the  model  results  with 
their  data,  the  data  points  from  figure  4-2  were  digitized  using  an  HP 
plotter.  Since  the  data  in  figure  4b  of  reference  30  were  the  only  group  to 
contain  H2  data,  this  group  was  chosen  for  modelling.  To  arrive  at  the 
rates  needed  to  model  their  data,  a  least-squares  approach  was  attempted. 
Because  of  the  2:1  correspondence  between  H2  and  SiH4  in  the  model,  only 
the  N-atom  data  and  the  SiH4  data  were  included  in  the  fitting  effort.  The 
approach  followed  the  Marquardt  algorithm  as  given  by  Bevington36,  but  in 
this  case  the  method  did  not  work  well.  It  was  found  very  early  on  in  the 
fitting  that  the  iterative  method  of  Marquardt  often  "over  shot"  a  par¬ 
ticular  rate  by  an  amount  large  enough  to  drive  the  rate  negative.  This 
condition  played  havoc  with  the  integrator  and,  of  course,  gave  unreal 
results.  To  suppress  these  excursions,  the  iterations  were  "damped"  so 
that  a  rate  could  only  be  increased  or  decreased  by  a  factor  of  ±0.8  (times 
the  rate)  during  an  iteration.  Another  problem,  which  was  never  solved, 
appeared  to  have  its  origins  in  the  shape  of  the  chi-squared  hypersur¬ 
face36.  Many  times  the  algorithm  would  begin  searching  for  the  chi- 
squared  minimum  along  a  path  that  pushed  one  particular  rate  to  a  very 
large  value  (much  higher  than  gas  kinetic).  An  attempt  was  made  to  solve 
this  problem  by  "capping"  the  rates  to  some  value  less  than  gas  kinetic, 
but  even  this  method  sometimes  failed  to  yield  a  minimum  in  chi-squared. 
The  worst  problem  of  all  was  that  the  chi-squared  minimum  found  by  the 
procedure  was  a  strong  function  of  the  starting  values  used  for  the 
iteration  of  the  rates.  Even  with  only  8  "unknown"  rates,  this  made  finding 
a  true  minimum  very  difficult.  Exacerbating  this  problem  even  further  was 
the  fact  that  the  calculations  were  by  nature  very  slow.  A  typical  "ad¬ 
justment"  of  the  rates  might  require  10  hours  of  desktop  computer  time, 
only  to  find  the  procedure  had  failed  to  find  a  minimum.  As  difficult  as 
this  procedure  was,  it  was  still  more  manageable  than  random  shuffling  of 
the  rates  to  fit  the  data.  In  the  end,  a  combination  of  "cut  and  try"  and 
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least-squares  was  used  to  fit  rates  in  the  model  to  the  data  of  Horie  et  al. 

Figure  5-5  shows  the  results  of  this  rate  adjusting  exercise  along 
with  data  from  reference  30  as  digitized  from  their  figure  4b.  The  con¬ 
tinuous  curves  are  calculations  from  the  model,  and  the  symbols  are  the 
data  of  reference  30.  The  "best  fit"  rates  are  given  below: 


ki=1 .74x1 0~12 
k4=1x10_1° 
k7=1x107 
k10=2.0x1 0-1 1 


k2=5.9x10-12 
k5=7.8x1 0-1 1 
k8=5.25x10-31 
kn=4.8x106 


k3=  1x1 0-1 0 
k6=2.2x10-11 
kg=8.0x1 0_1 1 


It  must  be  stressed,  however,  that  the  above  rates  may  not  yield  an 
absolute  minimum  in  chi-squared,  but  represent  the  rates  which  yielded 
the  smallest  chi-squared  observed  during  the  work.  Since  no  uncertainties 
were  given  in  reference  30,  the  actual  value  of  chi-squared  could  not  be 
calculated,  however  for  fitting  purposes  the  data  points  were  all  assigned 
equal  weights.  One  of  the  problems  in  trying  to  model  these  data  is  that 
Horie  et  al.  used  a  gas  mixture  of  N2  and  He  in  their  discharge,  but  they  do 
not  state  what  the  mixture  was  for  the  data  given  in  figure  4-2.  The 
model,  of  course,  assumes  the  gas  mixture  entering  the  discharge  is  pure 
N2.  The  effects  of  He  quenching  of  the  vibrationally  excited  N2  and  perhaps 
SiH4  dissociation  (or  ionization)  produced  by  collisions  with  He  metas¬ 
tables  would  need  to  be  included  in  the  model  if  one  were  to  try  to  ac¬ 
curately  model  the  data  of  reference  30.  The  vibrational  distribution  used 
as  input  to  the  model  might  also  be  considerably  different  if  a  large 
amount  of  He  were  used.  In  addition  to  these  factors,  the  measurements 
themselves  are  very  difficult.  One  of  the  difficulties  which  Horie  et  al.  do 
not  address  is  the  problem  of  the  introduction  of  the  SiH4  into  the 
afterglow.  As  can  be  seen  from  figure  5-4,  the  reaction  rate  is  a  strong 
function  of  the  SiH4  density;  therefore,  if  the  SiH4  is  introduced  through  a 
small  tube  (introduced  at  a  point)  the  higher  density  of  the  SiH4  at  the 
point  of  introduction  will  tend  to  "speed  up"  the  reaction  beyond  what  one 
might  expect  if  the  SiH4  were  more  evenly  distributed  upon  its  introduc¬ 
tion.  Considering  all  of  these  problems,  the  fit  is  reasonable. 

If  one  now  takes  this  rate  set  and  tries  to  fit  either  the  data  from 
figure  4-2(a)  or  4-2(c),  the  fit  is  very  poor.  Since  the  gas  mixture  used  in 
reference  30  was  not  known,  this  lack  of  consistency  could  be  explained 
by  assuming  that  any  of  the  reactions  involving  N2  would  be  changed  by  the 
introduction  of  He  into  the  gas  mixture.  Since  these  are  basically  reac- 
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COMPARISON  OF  MODEL  RESULTS  WITH  MEASUREMENTS 
OF  REFERENCE  30  (Fig.  4b) 


Comparison  of  Model  Results  with  Measurements  of  Reference  30 


tions  R-1  and  R-8,  if  the  data  from  reference  30  could  be  fit  by  only 
adjusting  these  two  rates,  then  the  model  would  at  least  remain  consis¬ 
tent.  This  procedure  proved  possible  and  the  "best  fit"  rates  for  the  other 
two  data  sets  are  given  below: 

Figure  4-2(a):  k1=3.07x1 0“12  k8=2.32x10-31 

Figure  4-2(c):  ki=5.34x1 CT12  k8=4.05x10'31 

One  point  that  arises  from  the  fitting  procedure  is  that  both  ki  and  k8, 
derived  from  the  fit,  are  somewhat  larger  than  was  expected.  In  all  cases, 
k8  exceeds  the  rate  calculated  from  the  Thompson  theory  and  ki  exceeds 
the  10-3  times  gas  kinetic  rate  for  the  v-to-E  process.  While  it  is  dif¬ 
ficult  to  assess  the  true  rates  from  the  rates  arrived  at  using  the  un¬ 
known  Fle-N2  gas  mixture,  at  least  the  rates  calculated  above  do  not 
exceed  the  basic  gas  kinetic  rate.  Both  the  general  behavior  of  the  model 
and  the  reasonable  fit  to  the  data  of  reference  30  give  credibility  to  the 
model  but  in  no  way  "prove"  the  model  is  the  correct  one. 

Behavior  of  the  Vibrational  Distribution 

During  the  reaction,  the  N2  vibrational  distribution  changes  sig¬ 
nificantly.  The  vibrational  pumping  supplied  by  reactions  R-5  and  R-9  is 
fast  enough  to  overcome  any  v-T  losses  for  the  levels  v=20  and  v=13,  and 
subsequent  v-v  exchange  results  in  all  the  vibrational  levels,  up  to  the  48 
considered,  being  excited.  Figure  5-6  shows  the  N2(X1X)  vibrational 
distribution  at  different  times  under  the  same  conditions  as  in  figure  5-5. 
Initially,  only  levels  13  and  20  are  pumped  but  as  the  reaction  proceeds, 
and  v-v  exchange  occurs,  all  of  the  levels  become  significantly  populated. 
The  behavior  of  the  vibrational  levels  above  v=9  is  especially  important 
because  these  levels  are  available  for  dissociating  SiFU  in  the  reaction. 
Figure  5-7  shows  the  time  dependence  of  this  group  of  levels  during  the 
reaction.  The  reason  the  group  "hangs  up"  after  the  reaction  is  complete  is 
due  to  the  small  v-T  rates  in  N2.  Certainly  if  diffusion  were  included  in 
the  model,  these  levels  would  be  relaxed  by  diffusion  to  the  walls  of  the 
tube.  The  peak  fractional  population  above  v=9  occurs  at  t=37.2  ms  and 
reaches  a  value  of  7.6x1 0-4.  This  number  is  more  than  three  orders  of 
magnitude  above  the  fractional  population  at  t=0. 
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CALCULATED  TIME  DEPENDENCE  OF  FRACTIONAL 
VIBRATIONAL  POPULATIONS  IN  N2 
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Figure  5-6:  Calculated  Time  Dependence  of  Fractional  Vibrational  Populations  in  Nitrogen 


FRACTIONAL  VIBRATIONAL  POPULATION  AVAILABLE  FOR 
SiH.  DISSOCIATION  DURING  THE  REACTION 


Final  Products.  Higher  Nitrides 

The  final  products  in  the  reaction  are  SiN,  SiN2,  H2,  SiH2  and  SiH4  (if 
complete  dissociation  does  not  occur).  At  lower  pressures  very  little  SiN2 
is  formed.  It  is  doubtful  that  in  a  real  reactor  SiH2  remains  for  very  long 
in  the  gas  phase.  It  may  be  dissociated  by  collisions  with  vibrationally 
excited  N2  (as  SiH4  is)  or  it  may  diffuse  to  the  walls  and  stick.  Dewhurst 
and  Cooper29  analyzed  the  white  powder  which  coated  the  tube  walls  and 
found  it  to  contain  mostly  higher  nitrides  of  silicon  of  the  form  SixNy.  The 
question  of  whether  the  powder  is  formed  in  the  gas  phase  and  "deposits" 
on  the  walls  or  whether  the  powder  is  actually  formed  on  the  walls  is  an 
interesting  one.  During  the  course  of  this  investigation  an  attempt  was 
made  to  observe  particulates  in  the  gas  phase  by  directing  a  HeNe  laser 
down  the  reaction  vessel  and  looking  for  scattered  light.  No  scattering 
could  be  observed,  but  the  sensitivity  of  this  method  is  unknown.  At  the 
higher  pressures  in  the  experiment,  it  might  be  possible  for  three-body 
reactions  to  occur  among  some  of  the  nitrides  to  form  higher  nitrides  such 
as 


SiN  +  SiN  +  M  -  Si2N2  +  M 
or 

SiN  +  SiN2  +  M  — ♦  Si2N3  +  M 

but  it  is  doubtful  that  complex  nitrides  such  as  ShN4  are  actually  formed 
in  the  gas  phase  during  such  short  reaction  times.  Reactions  between  some 
of  the  simpler  nitrides,  such  as  SiN,  and  SiH4  may  also  occur  in  the  gas 
phase  resulting  in  some  hydrogenated  nitride.  A  species  of  this  type  may 
have  been  the  cause  of  the  broadband  absorption  seen  in  the  FT-IR  ex¬ 
periment.  Because  of  the  low  densities  of  the  products  and  the  low 
pressures  in  the  experiment  the  most  likely  place  for  formation  of  more 
complex  nitrides  of  silicon  must  certainly  be  the  walls  of  the  tube.  It  is 
therefore  suggested  that  any  higher  nitrides  produced  in  the  reaction  are 
produced  on  the  reactor  tube  walls  and  not  in  the  gas  phase. 
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Chapter  VI 

Summary,  Conclusions,  and  Recommendations 

Visible  and  ultraviolet  spectra  of  the  chemiluminescence  from  the 
reaction  of  silane  with  active  nitrogen  have  been  obtained  over  the 
wavelength  range  230-850  nm  and  a  pressure  range  of  0.5  to  10  torr.  The 
spectra  at  higher  pressures  are  dominated  by  emission  from  Si  and  SiN  and 
at  low  pressure  by  an  unidentified  broad-band  emitter.  The  measurements 
indicate  that  both  the  Si  and  SiN  emission  intensities  increase  with 
pressure  more  rapidly  than  the  first  positive  intensity  of  the  N2 
afterglow.  The  measurements  also  show  that  the  vibrational  levels  of  the 
upper  electronic  state  of  SiN  are  populated  to  the  same  extent  as  those 
observed  in  the  reactions  of  both  SiCI4  and  SiBr4  with  active  nitrogen. 
From  the  general  behavior  of  the  chemiluminescence,  it  is  clear  that  the 
reaction  involves  some  kind  of  chain  reaction,  and  that  direct  reaction  of 
SiH4  with  N-atoms  can  be  ruled  out.  Silane  dissociation,  resulting  from 
the  reaction  with  active  nitrogen,  has  been  measured  downstream  of  the 
reaction  zone  and  correlated  with  the  relative  N-atom  density.  In  spite  of 
the  evidence  ruling  out  direct  reaction  of  SiH4  with  N-atoms,  these 
measurements  show  that  for  small  amounts  of  SiH4  the  dissociation  is 
proportional  to  the  N-atom  density  in  the  afterglow.  To  explain  this 
apparent  paradox,  a  reaction  model  has  been  developed  which  can  describe 
the  general  behavior  of  the  Si  and  SiN  chemiluminescence  and  the  silane 
dissociation  observed  in  the  reaction.  The  rate-equation  model  involves 
only  neutral  species  and  uses  as  initial  conditions  the  characteristics  of  a 
"typical"  nitrogen  afterglow.  A  key  reaction  in  the  model  is  the  dissocia¬ 
tion  of  SiH4  by  vibrationally  excited,  ground  electronic  state  N2.  This 
proposed  reaction  involves  a  transfer  of  vibrational  energy  to  a  dissociat¬ 
ing  electronic  state  of  SiH4.  While  reactions  of  this  type  have  been 
observed  in  other  systems,  no  direct  evidence  of  SiH4  dissociation  by 
vibrationally  excited  N2  has  ever  been  seen.  The  choice  of  this  reaction  as 
the  driver  of  SiH4  dissociation  was  arrived  at  by  eliminating  other  pos¬ 
sible  reaction  channels  which  could  dissociate  SiH4  at  the  rates  observed 
in  time  dependent  measurements30  and  without  perturbing  the  observed  N2 
afterglow  emission.  Dissociation  of  silane  by  vibration-to-vibration  (v-v) 
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exchange  with  N2(v)  must  also  be  considered  as  a  possible  mechanism,  but 
no  modelling  of  this  scheme  was  attempted.  If  v-v  pumping  of  SiH4  is  a 
significant  driver  of  dissociation,  the  pumping  must  be  enhanced  during 
the  reaction  since  the  observed  dissociation  increases  dramatically  during 
the  reaction.  This  implies  the  density  of  N2(v)  levels  responsible  for 
pumping  SiH4  must  also  increase  dramatically  during  the  reaction.  Because 
of  rapid  v-v  up-pumping  in  N2,  the  densities  of  the  higher  lying  vibrational 
levels  are  greatly  altered  by  any  mechanism  of  vibrational  excitation; 
therefore,  a  further  implication  is  that  if  v-v  pumping  of  SiH4  occurs  it 
must  come  from  the  higher  lying  levels  of  N2(v).  The  vibrationally  excited 
N2  is  produced  by  exothermic  reactions  of  two  nitrides  of  silicon  with 
atomic  nitrogen.  The  model  demonstrates  that  the  vibrational  pumping 
available  from  these  reactions  is  sufficient  to  drive  the  SiH4  dissociation 
at  the  rates  observed.  In  addition,  the  model  gives  the  SiH4  dissociation 
proportional  to  the  N-atom  density  in  the  afterglow.  The  model  includes  a 
complete  description  of  the  exchange  and  loss  of  N2  vibrational  quanta  and 
therefore  allows  accurate  modelling  of  the  time  dependence  of  the  SiH4 
dissociation.  The  model  results  have  been  compared  with  time  dependent 
mass  spectroscopy  measurements  of  SiH4  dissociation  made  by  Horie,  et 
al.30  in  a  similar  experiment,  and  the  model  can  reasonably  predict  their 
observed  results  for  a  given  set  of  conditions.  The  model  could  not 
describe  all  of  the  data  of  reference  30  under  all  conditions  because  the 
experiments  in  reference  30  were  performed  with  a  gas  mixture  of  N2  and 
He  and  no  He  was  included  in  the  model.  Specifically,  the  effect  of  He  on 
the  vibrational  levels  of  N2  and  SiH4  dissociation  induced  by  He  metas¬ 
tables  were  not  included  in  the  model.  In  general,  atomic  species  such  as 
helium  are  known  to  effectively  relax  molecular  vibrational  levels  so  it 
might  be  expected  that  significant  concentrations  of  helium  in  the  gas 
mixture  would  alter  the  behavior  of  the  N2  vibrational  populations  and 
therefore  affect  SiH4  dissociation.  Helium  metastables  should  not  directly 
perturb  the  amount  of  SiH4  dissociation,  since  their  density  is  probably  on 
the  order  of  1  part  in  106  of  the  He  density,  but  the  time  dependence  of  the 
dissociation  can  be  affected.  It  was  found  in  the  modelling  that  the  time 
of  occurrence  of  the  peak  in  the  chemiluminescence  and  the  rapid  drop  in 
SiH4  density  can  be  dramatically  affected  by  the  initial  rate  of  SiH4 
dissociation  caused  by  various  species  in  the  afterglow.  Thus  the  presence 
of  He  metastables  could  affect  the  time  dependence  of  the  SiH4  dissocia¬ 
tion  while  not  affecting  the  total  dissociated  density.  Horie,  et  al.  did  not 
state  the  specific  gas  mixture  used  in  their  experiments;  therefore,  it 
was  impossible  to  include  He  in  the  modelled  results.  While  there  is  no 
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guarantee  that  the  model  developed  here  is  the  only  model  that  will 
describe  the  observed  results,  the  model  is  plausible,  does  not  violate  any 
laws  of  physics,  and  is  generally  successful  at  describing  the  observed 
results. 

Application  of  the  model  reaction  scheme  to  commercial  reactors 
producing  silicon  nitride  (S^N^  cannot  be  strictly  made  due  to  differen¬ 
ces  between  those  systems  and  the  system  studied  here.  Since  most 
silicon  nitride  reactors  do  not  operate  in  an  afterglow  mode  as  described 
in  this  work,  but  instead  operate  as  a  low  pressure  discharge,  the  reaction 
scheme  under  reactor  conditions  is  expected  to  be  more  complex.  On  the 
other  hand,  this  work  was  begun  for  the  reason  that  very  little  SiH4 
dissociation  occurred  in  a  SiH4-N2  discharge  due  to  electron  impact,  thus 
many  of  the  reactions  described  here  would  be  expected  to  still  occur  in  a 
low  pressure  reactor.  In  addition,  many  commercial  reactors  use  gas 
mixtures  of  silane  and  ammonia  instead  of  silane  and  nitrogen,  thus  the 
reaction  scheme  in  these  systems  should  be  considerably  different.  The 
implications  of  the  model  to  the  formation  of  ShN4  are  that  the  reactions 
in  the  model  rapidly  lead  to  the  formation  of  atomic  silicon  and  SiN.  At 
higher  pressures  or  for  very  long  reaction  times,  SiN2  is  also  formed  at 
densities  comparable  to  the  dissociated  silane  density.  These  densities 
are  typically  a  factor  of  between  2  and  5  less  than  the  density  of  N-atoms 
in  the  afterglow  which  is  typically  less  than  1%  of  the  N2  density.  It  is 
clear  then  that  these  silicon-bearing  radicals  (along  with  smaller  quan¬ 
tities  of  S1H2  and  SiHa)  are  at  very  low  densities  for  the  formation  of 
higher  nitrides  by  three-body  processes  in  the  gas  phase.  Since  reactors 
for  ShN4  production  usually  operate  at  pressures  below  1  torr,  it  is 
doubtful  that  the  nitride  is  formed  in  the  gas  phase  and  "rains  down"  on 
the  substrate  surface.  It  is  far  more  likely  that  the  reactions  leading  to 
Si3N4  formation  occur  on  the  substrate  surface  and  rely  on  the  species 
generated  by  the  gas  phase  reactions  described  in  this  work  to  act  as  the 
"starters"  for  the  surface  reactions.  The  apparent  rapid  reaction  of  N- 
atoms  with  SiH2  (and  perhaps  SiH3)  means  that  even  under  discharge 
conditions  and  even  in  reactors  utilizing  NH3  the  simpler  nitrides  (espe¬ 
cially  SiN)  will  be  formed  and  be  available  for  surface  reactions.  The 
preference  of  most  commercial  producers  for  gas  mixtures  containing  NH3 
and  the  low  pressures  involved  may  indicate  that  species  such  as  SiN2, 
formed  in  reactions  with  N2  at  higher  pressures,  may  be  deleterious  to 
good  quality  film  production.  The  fact  that  surface  reactions  are  impor¬ 
tant  for  the  formation  of  the  higher  nitrides  opens  a  much  larger  arena  for 
affecting  the  final  nitride  product.  Specifically,  the  preparations  of  the 
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walls  of  the  reactor,  materials  involved,  and  even  the  geometry  may  alter 
the  film  quality  or  growth  rate. 

Recommendations  for  Further  Work 

As  is  often  the  case,  this  work  has  created  almost  as  many  ques¬ 
tions  as  it  has  answered.  Of  importance  in  any  future  work  is  accurate 
measurement  of  the  time  dependence  of  any  of  the  observed  species  and 
chemiluminescence.  Good  measurements  of  the  time  dependence  of  the 
chemiluminescence  would  be  particularly  useful  for  checking  the  validity 
of  the  model.  Particular  attention  to  good  flow  conditions  in  the  afterglow 
and  to  the  method  of  silane  introduction  are  imperative  to  achieving  good 
time  resolved  measurements.  Better  construction  of  the  afterglow  system 
could  give  more  uniform  gas  flow  and  perhaps  a  "shower  head"  system, 
similar  to  those  used  in  commercial  reactors,  could  be  used  to  introduce 
the  silane.  Mass  spectroscopy  of  the  reaction  products  might  help  shed 
some  light  on  the  higher  nitrides  which  are  formed  at  higher  pressures  as 
evidenced  by  the  infrared  absorption  discussed  in  chapter  III.  While  Horie, 
et  al.30  used  mass  spectroscopy  in  their  work,  it  is  not  clear  over  what 
mass  range  they  surveyed  the  reaction  products.  The  question  of  SiH4 
dissociation  by  vibrationally  excited  N2  is  another  area  that  might  be 
further  addressed  in  another  experiment.  Unfortunately,  it  is  difficult  to 
prepare  vibrationally  excited  N2  with  significant  populations  of  v>10  in  a 
manner  that  is  reproducible  and  measurable.  One  possible  experiment 
would  be  to  heat  N2  to  approximately  2500°K  in  an  oven  and  introduce  SiH4 
into  the  gas  stream  leaving  the  oven.  The  reaction  would  be  slow,  but  a 
slow  flow  system  with  diagnostics  downstream  of  the  reaction  could 
measure  either  the  SiH4  dissociation  or  perhaps  a  higher  silane,  such  as 
Si2H6,  formed  from  the  radicals  produced.  One  final  suggestion  for  further 
work  would  be  to  bring  a  technique  such  as  laser  induced  fluorescence 
(LIF)  to  the  diagnosis  of  species  formed  in  the  reaction.  It  would  be 
particularly  useful  if  the  species  SiN2,  which  has  been  proposed  as  one  of 
the  major  species  produced  at  higher  pressures  could  be  identified  in  the 
gas  phase.  Spectroscopy  of  SiN2  is  scarce,  however  an  ultraviolet  absorp¬ 
tion  system  has  been  identified  in  the  wavelength  range  300-360  nm  in 
matrix  isolation  at  4°K75.  The  estimated  densities  of  SiN2  in  these  ex¬ 
periments  range  from  6.5x1 01 2  cm-3  at  3  torr  to  2x1 01 3  cm-3  at  10  torr.  A 
technique  such  as  LIF  should  be  capable  of  measuring  densities  of  this 
order. 
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Appendix  A 


Heat  of  Formation  of  SiN 

The  heat  of  formation  of  a  diatomic  molecule  can  be  calculated  if 
the  dissociation  energy  of  the  molecule  is  known  and  the  heats  of  forma¬ 
tion  of  the  dissociated  products  are  known.  If  Hf(AB)  is  the  heat  of  forma¬ 
tion  of  species  AB  then, 

H,(AB)  =  H|(A)  +  Hf(B)  -  D0(AB)  (eq.  A-1) 

where  D0(AB)  is  the  dissociation  energy  of  AB  measured  with  respect  to 
the  zero  point  energy  of  the  molecule.  For  SiN  the  dissociation  energy  has 
not  been  clearly  established.  Herzberg44  gives  the  dissociation  energy  as 
D0=4.5  eV  from  the  work  of  Gaydon.  However,  in  his  later  book,55  Herzberg 
does  not  list  a  value  for  the  dissociation  energy.  Bredohl,  et  al.41  used  an 
extrapolation  of  their  calculated  ground  electronic  state  potential  which 
gave  a  value  of  D0=5.84  eV.  Since  it  has  only  been  since  the  work  of  Foster 
in  198443  that  the  electronic  energy  of  the  A2n  has  been  accurately  deter¬ 
mined,  another  look  at  the  dissociation  energy  was  performed  during  this 
work.  The  extrapolation  technique  given  in  Herzberg44  was  used  for  each 
state  for  which  Foster  gave  at  least  three  terms  describing  the  vibration¬ 
al  energy  levels.  The  vibrational  energy  levels  of  each  electronic  state 
were  calculated  until  the  vibrational  energy  reached  a  peak  value  at  some 
V max-  The  value  of  the  energy  at  this  point  was  considered  the  dissociation 
energy  De.  For  the  X2£  and  A2n  states,  the  calculated  energies  were  5.84 
and  5.37  eV,  and  for  the  D2n  state  the  calculated  energy  was  7.93  eV. 

Since  the  D2n  state  dissociates  into  Si(3P)+N(2D),  its  limit  should  lie  2.38 
eV  (the  N(4S)-N(2D)  energy  spacing26)  above  the  two  lower  states  used  in 
the  calculation.  Using  this  difference  gives  a  third  value  for  the 
Si(3P)+N(4S)  limit  as  5.55  eV.  Since  the  zero  point  energy  is  =0.07  eV, 
these  extrapolated  values  give  values  for  Do  between  5.30  and  5.77  eV.  On 
the  other  hand,  extrapolation  of  the  B2X  state  gave  a  value  of  De=4.69  eV. 
Since  the  B-state  also  dissociates  into  Si(3P)  +  N(2D),  this  would  give  a 
value  of  the  dissociation  energy  of  the  ground  state  as  4.69-2.38  or  2.31 
eV.  This  is  an  extremely  low  value.  One  might  argue  that  the  B-state  value 
is  anomalous  and  should  be  thrown  out,  but  in  light  of  the  uncertainty  over 
the  years  in  placing  the  A2n  state  it  is  not  clear  that  this  would  be  wise. 


Averaging  the  three  "good"  values  gives  a  dissociation  energy  of 
De=5.59±.24  eV.  Subtracting  the  zero  point  energy  gives  D0=5.52±.24  eV. 
Since  the  extrapolation  method  is  prone  to  give  high  values  for  the  dis¬ 
sociation  energy,  and  with  the  anomalous  behavior  of  the  B-state  ex¬ 
trapolation,  the  value  originally  given  by  Herzberg  will  be  used  for  the 
heat  of  formation  calculation.  Using  the  Hf(@298°K)  values  for  Si  and  N 
from  appendix  B, 

Hf(SiN)  =  107.7  +  112.9  -  4.5-23.03  =  117  (Kcal/mole). 

Considering  the  uncertainty  in  the  dissociation  energy,  this  value  must  be 
given  an  equal  uncertainty.  It  is  doubtful  the  true  value  exceeds  this  value 
by  very  much,  but  the  true  value  could  be  very  much  smaller  than  the  value 
given  above.  The  heat  of  formation  of  the  SiN  B2X(V=5)  state,  discussed  in 
chapter  III,  can  be  calculated  from  the  molecular  constants  of  Foster43 
and  the  above  value  for  the  heat  of  formation  of  the  ground  state.  This 
gives 


Hf[SiN(B2I(V=5))J  =  117  +  82.6  =  200  (Kcal/mole). 
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Appendix  B 


Heats  of  Formation  Used  In  This  Work 


Species 


Heat  of  Formation 


Referenc< 


H 

52.1 

76 

h2 

0.0 

76 

HSiN 

117 

77 

HNSi 

49 

77 

N(4S) 

112.9 

76 

N(2D) 

167.8 

76,26 

N(2P) 

195.2 

76,26 

N2(X1X) 

0.0 

76 

N2(X1Ii  V-10) 

62.8 

76,55 

N^X1!,  V-17) 

101.9 

76,55 

N^X1!,  V=20) 

117.3 

76,55 

N2(A3I) 

142.1 

76,55 

l\H 

85.2 

76 

nh2 

45.5 

76 

Si 

107.7 

54 

SiH 

91.7 

54 

SiH2 

68.1 

54 

SiH3 

47.8 

54 

SiH4 

8.2 

54 

SiN(X1Z) 

SiN(B2I,  V=5) 

117 

200 

(appendix  A) 
(appendix  A) 

‘Units  are  Kcal/mole  at  296°K 


Appendix  C 


The  Function  Pm  in  Equation  5-10 

For  the  function  P(T)  in  equation  5-10,  the  following  procedure  was 
used.  Capitelli68  gives  an  expression  for  the  time  constant  of  the  v=1 
relaxation  at  one  atmosphere  as 


t(T)  =5.7x1  O'1 2-Exp(234.9T'1/3) 


(eq.  C-1) 


where  x  is  in  sec.  and  T  is  in  °K.  In  terms  of  the  VT  rate,  t  is  given  by 


-,-1 


*(T)  =  Pn2(VT)i  ,0  (at  pN  =1  atm) 


(eq.  C-2) 


Using  equations  5-10,  11,  12,  and  14,  P(T)  was  calculated  in  terms  of  t(T). 
Two  hundred  values  of  P(T)  were  then  calculated  over  the  temperature 
range  150-1000°K  and  a  rational  function  of  the  form, 


2 

p^T)  _  a+bT+cT 
d+eT+fT2 


(eq.  C-3) 


was  least-squares  fit  to  the  calculated  values.  The  values  of  the  coef¬ 
ficients  were  found  to  be: 


a=1. 140616 
b=-0. 002328813 
c=4. 009294x1  O'6 


d=0. 7228253 
e=-0. 01295966 
f=5. 635459x1  O'5 


Using  this  function,  the  x  rates  could  be  calculated  to  within  1%  of  the 
values  given  by  eq.  C-1  over  the  temperature  range  150-1500°K. 
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